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SYNOPSIS 


The thesis presents the statistical analysis of field data 
on conpressive strengths of concrete cubes for characteristic 
concretes of M150, M200, M250 and M350. The statistical analysis 
indicates that the most of the groups of concrete follow normal 
distribution with one percent level of significance. The most 
probable value of the failure of concrete cubes, computed by the 
method of least squares, is found to be 5,5 percent. This value 
checks closely with that recommended by CEB-FIP Committee for 
international code of practice. Statistical analysis is performed 
for the data on strengths of high tensile steel and cold twisted 
deformed steel. Load survey carried out on floors in office rooms 
and statistical analysis of the same is presented. Statistical 
analysis of the data collected on geometric parameters of a typical 
concrete section is also included. The strengths of steel and 
concrete are always treated as random variables; whereas loads 
and geometric parameters are treated as deterministic in some 
examples and probabilistic in some other examples. 

The reliability analysis of simply supported and continuous 
prestressed concrete beams designed as per Indian Standard 
specifications is presented using statistical analysis of random 
variables and Monte Carlo technique. The reliability analysis of 
the above beams is done at limit states of (i) strength, (ii) initial 
crack, (iii) deflection and (iv) strength at transfer of prestress. 
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The failure of a prestressed concrete section at limit state 
of strength is divided intd two cases - (i) under-reinforced and 
(ii) over- reinforced. The probability of failure of a section is 
evaluated based on the above two conditional failure events. 

A method for the reliability based design of sinply si;5)ported 
and continuous prestressed concrete beams at limit state of strength 
for a given probability of failure is formulated. This involves 
solving of an integral equation and finding the mean value of 
ultimate resisting moment of the section. 

A rational study on load factors is made taking into account 
of the coefficient of variation of the resistance and probability 
of failure of the beams. Load factors for dead load and live 
load are suggested for strength and initial crack limit states. 

Limit state design of a prestressed concrete beam, based on send- 
probabilistic approach, is illustrated. 

The probability of failure at transfer condition of prestressed 
concrete beam designed by Indian Standard Code is fPund to be higher 
than that at limit states of strength and deflection during working 
load condition. The probability of cracking of the beams is generally 
hi^er than the probability of failure at limit states of strength 


and deflection 



CiyPTER 1 


INmjllJCTION 


1 « 1 General 

Ihe probabilistic approach to structural safely in civil 
engineering has been one of the subjects of interest in the last 15 
years* Great Strides have been made to accurately compute stresses, 
deflection, buckling loads etc* If loads and resistances aa^ known 
exactly, Ifeen Significant reductions in the factor of safety and 
is^roved economical benefits can be achieved. In fact, for maty 
problans, loads and resistances can be described as statistical 
variables and therefore, the safety of the structure can also be 
statistical variable* It would appear that improved structtiral 
analysis and design techniques should be accompanied by improved 
methods of predicting structural safety. It has been proposed for 
sometime that a rational criterion for the safely of the structure 
is its reliabilily or probabili'^ of survival* Statistical frequency 
distelbutions for loads and strengths must bo considered in determi- 
ning Ihe reliability* iTeudenthal said (4)* : ’Because the design 
of a structure embodies uncertain paredictions of the structural 
materials as well as of the expected load patterns and intenslMes, 
the concept of jacobabiliiy must form an integral part of aiy rational 

*Numbers in paranthesis refer to corresponding items in References. 
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design; any conceivable condition is necessarily associated with a 
numerical measure of the probability of its occurrence# It is by 
this measure alone that the structural signifioatxse of a specified 
condition can be evaluated*. 

Much woifc was done on the basic concept of stiuctural 
safety, in 1950s. In 1960s, most of the contribution was towards 
the probabilistic analysis and design of structural i^stems# In 
1970, CEB-PIP Committee (45) reccmmended semiprobabilistic limit 
state design method in which basic parameters are considered as 
random and the coefficient related to safely is introduced by 
using limit state method where increased loads are compared with 
the relevant resistance of the structure and where effects of 
service loads eire compaored with specified values. Currently Ang, 

Amin, Cornell and Und (l8, 24, 30, 34-, 37) have developed reliability 
design formats where uncertainty is expressed in terms of Ihe 
coefficient of variation without frequency distributions of individual 
variates. literature dealing with reliability analysis and design of 
structures is reviewed in the present thesis, 

1 .2 EEVIEW OP HIERAOUEE OK PRCBABIHSTIC AHAIYSIS AHD DESIGN 

Most of the literature available in Engli^ language is 
reviewed in chronological order* 


A comprehensive treatise on the safety of structures and 
on the factor of safety in particular was presented by Preudenthal 
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in 1947. The basic concept of structural safely analysis was first 
outlined by him (l )• (Ehis problem was later discussed by Asplund ( 2 ) 
and Bigsley ( 3 ). 

It was shown ( 4 ) that the probability of failure, p^, can 
be computed as follows. 


r(H < s) » p^ = / 

0 

[ 1 - EgCr) 3fg(r) dr 

(t.la) 


ys) fg(s) ds 

(l*1b) 


o 


in which fg(i*) and fg(s) are equal to the derivatives of disribution 
functions ^ 5 ( 1 *) jrespectively, E and S represent resistance 

and action (load) 3:espectively. Ihe p^ can be ea^jressed in tessns of 
a design factor of safety and thus is used to provide for a piobabi- 
listic interpretation of the fee tor of safety. 

Julian ( 5 ) presented statistical data on the yield strength 
and ultimate strength of A-7 steel, the yield strength of intermediate 
grade steel and ultimate compiressive strength of concrete. 

The paper (6) deals with the methods of providing engineering 
safety and additional safety necessary for social purposes. The 
methods considered fbr providing safety were by specifying (i) a 
minimum ratio of resistance to maT-iTniiin design load (ii) a maximum 
value for probability of failure and (iii) a minimum factor of 
safety coupled with maximum probability of failure. The concept of 
structural life in design was discussed* 
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Milik Ticky and Milos Yorliecek (?) foamulated the problen of 
safety analysis of reinforced concrete framed structures subjected 
to load from one source and with several collapse mechanisms and 
load from several sources. Pearson's type III curve was suggested 
for Idle ultimate strength of the structure. Ultimate strength of 
statically determinate and indeterminate structures were ccsupared to 
find the relation between statistical parameters of their ultimate 
strengths produced under identical conditions and demonstrated that 
variability in ultimate strength of redundat^t etruoture is lower than 
that of statically determinate one in all cases. It was also shown 
how the defcarmation ability of critical sections can also be taken 
into account in the study of safety of structures. Aa approximate 
method was developed for evaluating the multiple integrals for 
reliability expression. 

inhe work of the Task Committee on factors of safely since 
1965 were summerised by Ereudenthal, Garretts and Shinoeuka (s). 
Reliability functions of structures consisting multiple members were 
derived for the cases when the loads are applied at equal intervals 
or at prescribed instants euad the number of occurrences of the load 
is governed by a Poisson law. Hmerical examples were also given. 

In 1967 j Cornell (9) derived expressions for bounds on 
the reliability of a structure, having number of failure modes and 
subjected to number of successive loads, fbr independent, perfectly 
and imperfectly correlated loads and modal resistances. It wsus shown 
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that althou^ the precise determination of the reliability of structural 
system is a most complicated effort, upper and lower bounds can be 
readily evaluated, 

Moses and Kinser (lO,il) presented a paper on optimum sizing 
of elements for a multielement and multiload conditions of stiuoture 
with safety defined in tenns of allowable probability of failure. 
Significant wei^t savings were obtained by computing the statistical 
correlation between failure modes and developing an ordinary melixod 
to find probabilities of failure of a mode conditional upon the 
suiTrival of the other modes. Th^ also presented a method for 
evaluating the reliability of a structure (pinned -connected truss) 
with possible fail-ure modes subject to two load conditions taking 
strength and load as probabilitic, 

Tuifestra (12) formulated the choice of safety level in 
struotural design on the basis of minimum expected loss criterion. 
Bounds on sensitivily to failure probabilities had been established 
based on the insignificance of initial costs, Ihe collapse or limit 
failure of a system was examined by Shinozuka and Hanai (15) and it 
was shown that the exact reliabilily ejpression requires the evaluation 
of the multiple integrals. 

One of the important tools in reliabilily analysis is the 
hazard function or risk function. In 1968, Ang and Amin (l4) 
established a monotonic property for the hazard function of structures, 
33ie paper (l5) describes some of the advantages of rational probabil 
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analysis and design concepts with existing deterministic procedures. 

Bie design decision process for seiemdLc loading was illustrated using 
statistical decision theory. A Monte Carlo siinulation technique was 
described by Warner and Kabaila (16) which can be used to determine 
the cumulative distribution functions of stochastic variables such 
as ultimate strength and factor of safely. Using this the variabilily 
in strength of an exially loaded short reinforced concrete column was 
investigated and the results were compared with closed farm. solutLon. 

Procediires were presented by Sexsmilh (l7) for determining 
the probabilily distribution of the safely margin and related quantities 
for several 1ypes of reinforced concrete members. Dae procedure 
examines liie structure and its components at a factored load (near 
ultimate) and considers the probability of failure under this loading, 

A fonmilation of structural safely was proposed by Ang 
and Amin (is) in which the lack of knowledge and information is 
handled through a judgement factor whereas the statistical variables 
of information are treated probabilistically. 

Ihe paper (l9) deals with the selection of probabilily 
distributions for wind load, snow load and plastic moment of a section. 
She plastic collapse failure of a i^stem waus investigated, Ihe method 
of determining probabilily of failure of a frame by individual failure 
mechardsias was presented. It was shown that 'toe probabilily of 
failure of the frame involves the evaluation of multiple integrals. 



7 


Approximate methods were suggested for evaluating the probability of 
failure of a fraae. 

lUie advantages of probabilistic ooneepts were discussed by 
Shah (20), He explained a format for the calculation of probability 
of failure of the structure, using only mean values and standard 
deviations of resistance and load and assuming probability distribution 
of (R~s) or (R/I). 

Costello and Chu (21 ) assumed Beta distributions for strengths 
of steel and concrete and based on these, probabilities of failures of 
reinforced concrete rectangular beams were presented. 

Ihe paper (22) illustrated various difficulties, such as 
(i) choice of the appropriate probabilily model, (ii) introduction of 
subjective elements into the probabilistic structure and (iii) inter- 
pretation of probabilistic information in a form vdiich leads to 
rational decisions, in the application of probabilistic concepts. 
Benjamin and land (23) investigated the approximate probabilily measures 
associated with the existing code provisions. Reliabiliiy based 
procedures were illustrated using a set of assigned probabilities. 

The concept is to relate probabilily measures to existing deterministic 
code provisions thereby allowing improvement of design procedures 
within the deteministic existing code format, 

Cornell (24) presented a consistent first order reliability 


code format based on first and second moments of all stochastic variables. 
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Yao and Ydi fonmlated the safely analysis of structural 
systems wMch yields a systaaatic method of counting the failure paths. 

The reliability of a parallel and redundant ^stem was foimilated in 
this manner (25). Some examples (trusses and frames) of the optimal 
reliability based design of both weakest lirk structures and frames 
for minimum wei^t were presented by Stevenson and Moses (26). Shey 
developed a method of analysis (27 ) tor determining the overall probability’' 
of failure of fraaed structures (suitable for plastic mechanism analysis) 
vhose failure mechanisms is expressed as a linear combination of several 
structural resistances and loadings. The correlation between failui*e 
modes was taken into account. To determine the reliability of a failure 
mode, recursive integration method and Pearson distribution method were 
developed. 

Allen (28) made a detailed study on es-fcLmated probability 
distributions for the ultimate bending strength and ductility ratio 
(curvature at ultimate to curvature at yield) using Monte Carlo 
technique for rectangular concrete sections reinforced in tension only. 

He concluded that a section under-r’einf or ced according to ACI 318-63 (29) 
will undergo a brittle compressive failure. The -variabili-ty of ductility 
ratio is much hi^er than "the -variability of ul-timate moment. 

land (30) presented reliability code formats of arbitrary 
higher order than first and second moments of stochastic -variables 
and demonstrated the calibra-feLon of a partial safety factor format to 
Cornell's format as well as to the Ang-Amin format. 
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{Qie paper (51 ) proposed a distribution function composed of 
a central gaussian portion and an exponential tail with a sin^e point 
of discontinuily as a convenioit representation of statistical demand 
and capacity data to evaluate reliability, regardless of any conjecture 
about the parent distributions. !I!he convolution integral of the 
probabilily of failure is evaluated in closed form in terms of simple 
statistics of the data. 

In 1972, !Pa^ committee on structixral safely reviewed the 
available literature till December 1971 and compiled a bibliograply 
on structural safety (32). The feasibility of designing concrete beams 
for mi nimuTn cost with a reliability based constraint was demonstrated 
by S.S. Eao (33 )• Instead of speci::^ing small values for probability 
of failure, the statistical constraint was specified as the design is 
safe when the expected moment capacity of the beam exceeds the design 
moment by a certain number of standard deviations. 

The importance of risk and uncertainty conditions in 
structural evaliiation and design was emphasized and practical methods 
for ride assessment in terms of failure probability and for development 
of reliability based design criteria were described and developed 
by Ang (34). He suggested the systematic analysis of uncertainty 
through first order statistical analysis and the explicit use of 
failure probability as a measrore of risk and safety. 


Ered Webster (35), in 1973, analysed the structural 
response of a reinforced concrete portal made up of one beam and two 
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spirally reinforced columns using probability tbeoiy. ilhe portal was 
subjected to an increasing sequence of static test loads and at each 
load the probabilities of plastic hinge formation were calculated. 
Possible paths of failure mode were represented in a Markov chain. 

A probabilistic based structural design method was introduced by 
Paloheimo Eero and Hannus Matti (36). Tae reliability is replaced 
by a deterministic design equation in sdiich all random variables are 
substituted by appropriate design values which depend upon the mean 
values and standard deviation of the corresponding variables, the 
reliability factors that reflect the acquiried levels of reliability 
and the sensitivity factors that measure the effects to individual 
random variables upon the total system. Methods of application were 
also considered. 

Models of risk evaluation were developed by Ang and 
Cornell (37). ihe distribution function of (E-S) or An (e/S) was 
considered without the distributions of individual variates. It was 
emphaisized that the distribution of the safely margin (e-s) or An (e/s) 
is important in the calculation of probability of failure. Design 
formats for reliability based design were paresented. Uncertainty 
(including uncertainty associated with errors in estimation and 
perfection in mathematical model) was expressed in terms of the 
coefficient of variation. This is amenable for systematic evaluation 
Finfl analysis of uncertainty in structures. Eoimulation for multiple 
load cases was also presented, Ihe quantitative analysis of design 
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uncertainties was illustrated by ELlingward and Ang (58) and was shown, 
how these uncertainties affect tiie level of risk, ttie risk associated 
with the existing design procedures were evaluated with reference to 
reinforced concrete. It was found 1iiat the percentage of steel and 
effective depth are the only variables contributing significantly 
uncertainty to resisting moment, Daese risks serve as initial basis 
for the formulation of rii^ based design. Using the reliability based 
design method proposed by Ang and Cornell, design of several beams, 
colimins and cold formed steel beam elements were illustrated by 
Ravindra, Lind and Si (39 ) and was shown that code parameters such as 
the safety indices, can be selected to match the safety level of current 
designs. Only mean values and standard deviations were used to represent 
the design random variables, Glie final design criteria were developed 
on a probability distribution free basis. 

The design formats suggested by Ang and Cornell were ■used 
by Moses for reliability analysis of structural ^sterns, Moses (40 ) 
introduced partial safety factors which are quantities by which the 
element safety factor should be changed to accoiuat for its presence in 
a structural assemblage so that the overall probability remains roughly 
equal to what is desired for the element. Examples of solutions for 
reliability analysis and design of structural systems were presented. 

Probabilistic dynamic approach and problems that are 
encountered in dealing with the random process were discussed by 
Shinozuka (41 ). Bfodels that have been used in the structural response 
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analysis involving winds and waves were presented. Yao James ( 42 ) 
presented a simple reliabiliiy based procedure for the fatigue design 
of structures, 

Putcha C3aandraseikar and Dsyaratnam (45^44) piresented the 
reliabili-ty analysis of prestressed concrete simply supported beams 
designed on the codes of Indian Standards, jSmerican Concrete Institute 
and Keeamnended Code of Practice by CjEB-PIP, Probable strengths of 
basic materials used in the beam were generated as random variables 
subject to code specifications, Ihe probabilily of failure was calculated 
based on the governing equation of the ultimate strength of the section 
fixed by deterministic analysis and normal distribution was assumed for 
the strength and resistance of the section. !I!he design of prestressed 
concrete beam for a given reliabilily was farmulated as a polynomial 
equation. Optimization of prestressed concrete beams subject to the 
constraint of probabilily of failure besides the usual constraints was 
also done, 

1,3 STATMMT OP IHE HIOBLM 

®ie object of the present investigation is as follows : 

(a) Statistical analysis of the available data on strexjgths of 
materials, loads and geometric properties of the section; 

(b) Presentation of a method of reliability analysis of prestressed 
concrete sections at limit state of strength; 

(c) Probability failure analysis of prestressed concrete simply 
supported and continuous beaas designed as per the existing 
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I.S. Code (46) at limit states of (i) strength at design 
load, (ii) crack initiation, (iii) deflection and (iv) strength 
at transfer of prestress} 

(d) Design of simply supported and continuous prestressed concrete 
beams for given probability of failure at limit state of 
stiangth; 

(e) Study of the variation of load factors with coefficient 

of variation of resistance of the section and probability of 
failure and to suggest the values of load factors. 

Die thesis presents the statistical analysis of field 
data on compfessive strength of concrete cubes for characteristic 
concretes of Ml 50, M200, M250 and M550 and the most probable value 
of the failure of concrete cubes is computed using the method of least 
squares. Statistical analysis has been performed for the data on 
strengths of hi^ tensile steel and cold twisted defortnei steel, lead 
surv^ carried out on floors in office rooms and statistical analysis 
of the same is presented. Statistical analysis of the data collected 
on geometric properties of a typical concrete section is also included, 

A general method of analysis of probability of failure of 
a prestressed concrete section at limit state of strength is presented, 
Monte Carlo technique is used to stu<^ the reliability of sections at 
limit states of (i) strength, (ii) initial crack, (iii) deflection 
and (iv) strength at transfer of prestress. Reliability analysis of 
typical simply supported and continuous prestregsed concrete beams 
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for deterministic and probabilistic loads at different limit states 
are presented considering probabilistic variations of (i) strengths 
of concrete and steel and (ii) strengths of concrete and steel and 
geometric properties of the section, A method of design of simply 
supported and continuous prestressed concrete beams for given reliaibilily 
at limit state of strength is formulated and examples are given. 

Variation of load factors with the coefficient of variation of resis- 
tance and probability of failure are studied and the values of load 
factors for dead load and live load are suggested based on the results 
of the reliability analysis of prestressed concrete beams. Semi- 
probabilistic limit state design of prestressed concrete beam is 
ill;astrated. Conclusions drawn based on the discxjssions of the results 
are presented at the end, 

Uie main purpose of the thesis has been the demonstration of 
possibilities furnished 'j„ he i'ioate Carlo simulation method. Since the 
data on random behaviour of various quantities are currently not available 
in sufficient amount certain simplifications has to be accepted in the soluti 
Daey mainly concern the mutual dependence of some random quantities and the 
deterministic models of limit states. The general approach presented in the 
thesis can easily be developed after sufficient data are available. 
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STATISTICAL AEAHSIS OP SIERMOTBS, OP CONCEETE 
MD SQSEL, lOAD AMD DBIEMSIONS OP 1 SECMOM 

2.1 IMERODaCHON 

Acceptable strengths of materials and tolerances and charact- 
eristic loads are specified by the codes of practice. A rational 
approach to the specifications on structural safely needs to be 
associated with the field data and statistical analysis. A proba- 
bility analysis and design requires a knowledge of the statistical 
nature of strengths of materials, loads and the geometry of 
invidual sections, 

lihe collection and statistical analysis of field data on 
strengths of concretes used indifferent projects is first 
presented*. 3he freouency distribution of the samples and the 
probability of failure (pp) of a specified concrete are computed 
and presented. A characteristic coefficiait, which ensures a 
preassigned reliability of the test results less than the 
characteristic strength, is obtained by the meihod of least 
squares. The work also investigates the relation between 
characteristic and mean value of the strength of concrete used 
in some parts of India and compares it with that suggested in 

* ilhis work is reported in refer^ence nrmber (62). 
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CEB-EEP recommendations. The statistical analysis of data on 
strengths of high tensile steel (hts) and cold twisted defoimed 
(cTD) steel is presented with the probabilily distributions of 
the same, The survey carried out on floor loads in office 
building at Indian Institute of Technology, Kanpur, is also given 
and the frequency distribution of the load is specified. Lastly, 
a study of the variations in geometric properties of a typical 
prestressed concrete section is included and the statistical 
analysis of the data is presented. 

2.2 STATISTICAL ANAITSIS OP STRENGTH OP CONCRETES 
2,2,1 Collection of Pield Lata 

Structural Engineering Laboratory at the laadian Institute 
of Tfechnology, Kanpur,, has been associated with testing of variety 
of materials used in the construction industry. Liring the last 
ten years a large number of samples of cement, concrete and steel 
have been tested in the laboratory, 150 mm cubes were prepared as 
per Indian Standard Specifications (47 ) and supplied either by the 
builders or by the project supervisor® to the laboratory. The 
projects include small and large scale office, residential and 
workshop type buildings, bridges and foundation structures. This 
thesis presents the strength of the cubes after 28 days of water 
curing.Only those tests with a mininum number of 35 specimens have 
been considered for the analysis, Ohe specimens of a given project 
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were supplied in sets of three to six at rajadom. Bierefore, all 
the cubes of a given project were tested as and when received and 
not in a continuous series. Ihe cubes were tested as per the 
Indian Standard Specifications (47 ). Olie characteristic strength 
of the concretes vdiich are also called as the specified strengths 
of the concretes consist of Ml 50 and M200 based on the nominal 
mix specifications and M250 and M350 proportioned by a design method. 

2.2.2 Histograms and Statistical Analysis 

Die total number of samples of a particular concrete mix 
(say M150) are divided into project groups. Each group is 
separated into a number of sets as supplied to the laboratory. 

Each set consists of usually three samples and sometimes four 
or sdx. The mean value and the standard deviation of the strength 
of concrete of each set is confuted and typical variations are 
given in Jigs, 2.1 to 2.4, The figures indicate the random 
variation of the strength of the concrete. Histograms for three 
lypical groups a 2 *e shown in ligs, 2.5 to 2.7 indicating the 
frequency distribution of the strengths of the samples. The 
mean value and the standard deviation of the strengths belonging 
to the three groups are given in these fibres. 

All the samples belonging to a class of concrete are 
combined and the histogram for such a concrete is drawn. 
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ilga, 2.8 to 2*11 give the histograms and camzlative distributions 
of the strengths of the concretes. Ihe approximate nmber of 
intervals for drawing the histograms is selected (48) as 

z » 1 + 3*3 (2.1) 

in which z is the number of intervals and n is Ihe number of 
samples. A convenient value of 1he interval is chosen &om the 
appropriate value of z . 

The suitability of a mathematical model to fit the data is 
arrived at after applying the chi square test with one percent 
level of significance. The value of the chi square is computed 
from the formula (49) 

where 

X*" * the value of chi square 

= the observed frequency in the ith interval 
= the expected frequency corresponding to the model 
in the ith interval 

a = the number of categories considered. 

The number of degrees of freedom is given by 


2 ? 

C = I Y 

i=1 ^i 


( 2 , 2 ) 


M = a - r - 1 


(2.3) 
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wheire M = the number of degrees of freedom 

r = the number of parameters estimated from the data 

The chi sq^uare values computed fc^om Bq, 2.2 for Ihe assumed 
distributions are checked with the standard chi square values gLven 
in the standard tables (49). The normal distribution is found to 
satiaJ^ the chi square test for M200, M250 and M350 concretes at 
one percent level of significance, Hbwever, the normal distribution 
is found unsatisfactory for "toe complete class of M150 concrete 
samples, even though three out of the six groups satisfied the test. 
3he lognormal distribution is found satisfactory for M150 concrete 
samples. The computed values of chi square and the values of 
chi square obtained frcan the standard tables (49 ) for four classes 
of concretes are given in table 2.1. 

2# 2, 3 Analysis of Probability of Pailure 

The probability of failure of concrete is given by 

®cuk 

Pfe* ‘ “cuk) “ I ®cu 

fmOO 

where p^^ =» the probability of failure of concrete cube 
= the strength of concrete cube at 28 days 
‘^cuik” characteristic strength of concrete cube 
f(cf^u)" the density function of strength of concrete. 

Por a noamal density flinction, the Bi, 2.4 reduces to 
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P 


cu 



a 

cum 


s 

0 


in which 


(2,5) 


= the cimilative distribution function of the 

standairdized normal variable (50 ) 

0 ,s = parameters of the concrete naneiy mean value and 
cum ^ G ^ 

standard deviation 
%tan expressed as (49) 

n 

Ocum " 

8 * y Co . a ) /(n— 1 ) 

e cToi cum' ' ' 

1=1 

itor lognoimal distribution (48), the Bq,. 2*4 simplifies to 
log 


(a.6) 

(2.7) 


^cu ^ s 


] 


cl 


( 2 . 8 ) 


in which and s^^ are the parameters of the lognormal distri- 
bution, Ihey are expressed as 


°cul “ ‘^cum ^“^*5 s^^) 


(2.9) 


®cl ~ ( 2 . 10 ) 

where 6 is the coefficient of variation of a„„« Table 2.2 
c cu 

shows the probability of failure of the concrete for different 
classes of concrete, computed from Bjs, 2,5 and 2.8. The percentage 
of the tests of concrete cubes satisfying the Indian Code specifications 
on quality control (51 ) is also listed in the same table. 
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2*2#4 Relation Between Characteristic and Mean Strengths 
of Concrete 

Dae aim of a builder is to produce a concrete in the field 
whose strength is as close as possible to that of characteristic 
strength. Hbwever, there will be certain variation of the mean 
strength of the test specimens from the characteristic. One of 
the aims of this thesis is to find the relation between characteristic 
and mean value of tiie strength of concrete. (this thesis establishes 
■the expected probabilily of failure of the concrete for a given 
characteristic strength, Assuming the noiuial distribution, Bq., 2.5 
gives the probability of failuia of the concrete. This equation 
can be rearranged as 

( ’’ . cuk - °cum ) ^ ^-1 ( ) = ^ (2.,1) 

' s '^cu 

c 

in #iich k is the coefficioat depending on the probability of 
failure p . The Bi. 2,11 can be rewritten as 

0,»o -ks=- (l- k6) (2.12) 

cuk cum c °cum 

or 

B = 1 - k6 (2.13) 

vhere 

^ = the coefficient of variation 
® ~ cm 

The Bq, 2.12 yields the relationship between the mean strength 
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and characteristic strength for the preassigned reliabilily. Using 
normal distribution and substituting the corresponding values of 
and of each group of samples in Bi- 2,12, the values 

of 6 and k have been computed and given in table 2*3 • 

It can be seen from table 2,3 that the values of k vary from. 
1.21 to 2.01. Iherefore, it is desirable to establish a particular 
value of k which can be used in the design. Ihis value of k is 
obtained by using the method of least squares, based on the field 
data. The CEB-PIJ Committee has recommended a five percent proba- 
bility of failure based on which the value of k works out to be 
1,64 assuming normal distribution for the strength of the concrete. 

2,2,5 least Square lit 

Bor each group of samples, the values of B and 6 are computed 
and shown in ELg. 2,12. The value of 3 is considered as a random 
variable fbr an assigned value of 5 , A linear regression model 
connecting 3 and 6 is selected based on the assumptions that the 
values of $ corresponding to a given value of 6 follow normal 
distribution and -the variance of the values of B for a given 6 is 
independent of the magnitude of 5 , 33ie equation of the model, 

0 = b^ - k« (2.14) 

in which and k are parameters of the model, is selectetf 
in oonfoimiiy with that of the eqmtion recommended by CEB-Pip(45 ). 
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Table 2.3 • Values of 6 and k for various groups of concrete 


Sroup Hb. 

8 

k 

& 

1 

0.571 

1.7463 

0.2450 

2 

0.628 

1.4386 

0.2580 

3 

0.635 

2.0100 

0.1820 

4 

0.740 

1.5173 

0.1720 

5 

0.674 

1.4289 

0.2290 

6 

0.520 

1.7561 

0.2740 

7 

0.687 

1.6650 

0.1880 

8 

0.842 

1.1967 

0.1331 

9 

0.817 

1.5704 

0.1177 

10 

0.783 

1.5243 

0.1360 

11 

0.860 

1.2100 

0.1131 
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However when the variance is zero, it is known that the mean value 
of the teat results should coincide with the characteristic value. 
GMs leads to B = 1, when 6=0* Substituting the above values 
of 6 and 6 in Bi* 2.14, one gets b^= 1, Hence Hi, 2,14 reduces to 
the Bi* 2.13 with one parameter k* Using the method of least 
squares, k is chosen to minimize the sum of the squared differences 
between the observed values 0^ and the estimated expected values 
given by Bq, 2.13, that is k is chosen to minimize 

H 2 

I -(l - k6^) 1 (2.15) 

i=1 

where N is the nimber of the groups of the data. Taking 
derivative of Bq. 2.15 with respect to k and equating to zero, 
the value of k is given the following equation. 

H - N 

I «i - Z q 6i 

k = — (2.16) 

I «i 

i=1 

The computed value of k from Bq. 2.16 is equal to 1.6. Rearrange- 
ment of tile Bq, 2.11 gives 

The value of p^^ is equal to 5.5 percent for the established 
value of k equal to 1.6. !Ihe line of regression of $ on 6 
for the Bq. 2.13 is shown in the ELg. 2.12. 
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2.2,6 Regression Analysis 

!I3ie regression analysis is clone to indicate the level of 
confidence and I4ie acceptability of the most probable value of the 
failure of concrete cubes established earlier. 


The standard deviation of g , Sg , is computed from the 
following eqiAation (49) 


2 



^ 7 

I [Bi-Cl-ted,)! 

i=^1 


( 2 . 18 ) 


and it is fo\ind to be 0,0412. The variance of the mean value 
2 

of 6 , Sg^ , is given by 


2 

2 _ 

®8m ~ N 


The confidence inteoTval for the estimate of the mean value of 
is given by 


3-ts- <0 <0 + tSo 

0m — "^m — m 0m 


in which 


(2.19) 

6 


(2.20) 


0 = the mean value of B 

m 

t =» the ’t' statistic 

0 = the possible range of 0 . 

m m 

The number of degrees of freedom of the problem is ecLual to (f-i ). 
For one percent level of significance, the confidence interval for 
the estimate of 0^ is found from the standard table (49) and it 
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varies from 0.6659 to 0,7499. Die correspording confidejace limit 
lines are shown in ilg» 2.12. 

The variance of k is con^juted using the value of Sg in 
the following equation 


2 

H 


®k U 


( 2 . 21 ) 


in which 




I 

x=1 

= the variance of k 


6 = the mean value of fi 

m 


The confidence interval for the estimate of k is given by 


k - tSj^ i h ^ k + tSj^ (2, 22 ) 

in which k is the possible range of k. Ttor one percent level 
of significance and ten degrees of freedom, the confidence interval 
for the estimate of k is from 0,880 to 2.320. It can be seen from 
the table 2.3 that all the values of k for different groups of the 
data lie within the upper and lower confidence limits. The result 
of the significance test establishes the credibility of the 
estimated value of k. 


The correlation coefficient P is calculated from the 
formula (49) 


N 

I 

= i=1 


I 


/-S — ^ — s — - 

^ I < I 

i=1 ^ i=1 


(2.23) 
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and its value is equal to 0,9l« 5br ten. degrees of freedom the 
value of p f obtained from standard tables (49) for one percent 
level of significance, is less than the calculated value. Hence 
there is correlation between 3 and S , 

2.5 SIAHSHCAl MAIZSIS OH SDREHOIH OH STEEL 
2.3*1 Histogram and Statistical Analysis 

The structural esagineering laboratory at III Kanpur has 
conduojed number of tests on strengths of different steels. This 
thesis presents statistical analysis of the data on strengths of 
hi^ tensile steel (hts) and cold twisted deformed (cTD) steel. 

Ihere are seven groups of data on 7mm ^ HTS wire, Humber 
of specimens in each group varies from 2 to 22, All the samples 
belonging to each group are combined and the histogram and 
cumulative distribution of the ultimate strength of 7mm <P HTS wire 
are given in Jig, 2.13. The normal distribution is found to 
satis:^ the chi square test at one percent level of significance 
(satisfies for 5^- significance level also). 

Horty tvro specimens of 5mm <{• HTS wire belonging to one coil 
have been tested, Qhe histogram and cumulative distribution of 
the ultimate strength of 5fflm <l> HTS wire are given in Hig, 2.14. 
normal distribution is found to satis:fy the chi square test at one 
percent level of significance (satisfies for 5^ significance level also). 



f* of Observations 



relative Frequency 







176.5 178.5 




180.5 182.5 

(kg/mm^) 

(b) CUMULATIVE DISTRIBUT 
FIG. 2. 14 HIGH TENSILE 5TEE 
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There are 8 sets of data on strength of CTD steel* IJumber 
of i^edmens in each set is 6. Nominal diameter of Bjiecimens is 
equal to 25mm for four sets, 28tnm for 2 sets and 52mm for one 
set* As the characteristic strength of CTD steel upto 32mm i{> 
is same (52), all sets have been combined to fbrm a group. Ultimate 
strength and proof strength have been recorded. Ih.stogram and 
cumulative distribution of the data on proof strength and ultimate 
strength of torsteel are given in figs, 2, 15 and 2. 16 respectively. 
Normal distribution is found to satis:fy the chi square test for 
the data on proof strength at one percent level of significance. 

It has not been possible to fit ary standard distribution for the 
data of fig. 2.16, 

2.3*2 Analysis of Trobabilily of failure 

for a normal density function of ultimate strength of steel, 
the probability of failure of steel, p^ , is given by 

Ps= (2.24) 

s 

where o cr and s are the characteristic value, mean 
SK! SUL s 

value and standard deviation of the ultimate strength of steel 
respectively. Table 2,4 shows the probability of failure of HTS 
wires and CTD steel computed from Eq. 2,24. It is observed that 
the characteristic strength specified by I, S. Code (46 ) for 5mm <ti 
HTS wire is very much on the conservative side. The probabilily 
of failure for proof strength of CTD steel is observed to be hi^. 
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2.4 STAO?isan;cjffi ahahsis op loads 

A survey of floor loads in office roans has been carried 
out. Hhe live loads in the surv^ include the occupants, the 
movable partitions, all furniture and their contents and other 
miscellaneous items viz typewriters, table fans etc, A volume 
estimate has been made of the contents of cabinets, desks and book 
cases and a note has been made of the nature of the contents. 

After observing liie arrangement of files, stationeries and other 
items in a few almairahs, the maximum volume of an almairah that is 
usually used has been found to be 67 percent. All the rooms have 
been surveyed during peak working hours to record the number of 
persons woikxng and visitors and the furniture and storage in each 
room. The wei^t of individual items have been assessed as follows. 
Chiairs ard tables have been weighed. Weight of empty steel almairah 
has been obtained from the manufacturer, IHie almairah has been 
assumed to he occupied by paper to an estimated volume of 67 percent 
of its total capacity , An average weight of a person has been 
estimated after knowing the weight of a few persons. Other items 
either have been weired or assessed. The total load has been 
calculated after adding the weight of individual items. 

The magnitude of the response of floor changes with the 
position of the load and its distribution. Concentrated load 
effect should be t^en into account in determining the static 
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equivalexib floor load* An analytical study has been made for 
different magnitudes and occupancies of load on the floor to deter- 
mine the static equivalent load. llreating floor as a plate, 
sinply supported on all four sides as shown in lig^ 2.17, the 
deflection, u, of the floor at ary point (x,y) can be represented by 

u = sin ~ sin ^ (2,25) 

I si b 

where is the coefficient of the series defining the deflection 
surface and 1 and b are length and breadth of floor. IXiiing load 
survey it has been observed that all heavy loads are arranged near 
the walls of the room. Assuming ’A ‘area of the floor is occupied 
as shown in lig. 2.17a, consider two ^sterns of loads as shown in 
Jigs. 2.17a and 2.17b. The hatched portions shown in ligs. 2,17a 
and 2.17b represent the area of the floor occupied by loads. The 
total load is same in both cases, the method of virtual work 

concept the effect of both loads can be con^ared and the static 
equivalent load, Q^, distributed over the whole area of the floor 
can be found out, Hie virtual work done, V, by the load Q, 
distributed over the area 'A 'as shown in lig. 2.17a is given by 

V=-^ 11 sin ^ dx dy (2,26) 

The hatched area in lig. 2.17a represents the area of integraHon, 
Tbr fifty percenb of the total area of the floor occupied by the 
load as shown in Pig, 2.17a, the value of T, computed from Bq. 2,26, 
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is ). The virtual woife done, V , by Q (Pig. 2.17b) 

= 4r I I •? ^ dx dy 

_ 

TT^ 

Equating V = V , the value of Q , the ratio of Q to Qj is obtais®*^ 

as 0*412. Similarly, the values of have been calculated and 

given in table 2.5. During load surv^, as stated earlier, heavy 

loads are close to walls and it has been observed that load distri" 

bution in most of the rooms belong to cases 6 and 7 of table 2.5. 

The value of Q for these cases is less than 1. However, to be 

on the safer side, the value of Q is taken as equal to Q i. e. the 

total load is considered to be unifoimly distributed over the whol® 

area of the room (This calculation of Q is valid when floors are 

supported by walls. In the caseof large halls where partition 

walls are built and floors are designed as continuous slabs and 

flat slabs, the value of will be greater than one). Based on this, 

equivalent uniformly distributed floor load has been calculated by 

dividing the total load in each rocmi by its area. The details 

of room usuage, total load, area of the room and the equivalent 

viniformly distributed load in each room, are given in table 1 in 

2 2 

Appendix , The areas of rooms vary from 14m to 147, 8n • The 

2 2 

equivalent load intensily varies from 56 kg/m to 293 kg/m . Slie 
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Table 2,5. Values of Q for different distribution of loads on 
floor ^ 


Case 

Load Distributed 
Close to Walls 

Load Distributed 
at the Middle Portion 

<^r 

Area 

Magnitude 

Area 

Magnitude 

1 

bOfo 

Q 

KIL 

- 

0.412 

2 

73% 

Q 

ITEL 


0.557 

3 

50^ 

0.75Q , 

. 30% 

0,S5Q 

0.706 

4 

50^ 

0.50Q 

3<Ofo 

0.50Q 

1.0 

5 

60^ 

0.50Q 

AOfo 

0.50Q 

1.125 

6 

60^ 

0.75Q 

^o% 

0,25Q 

0.812 

7 

75^ 

0.90Q 

23% 

0.10Q 

0.954 
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influence of room area on load is not considered. Ilg. 2.18 presents 
histogran and cvimulattve distribution Ibr equivalent unifomOy 
distributed load in office rooms such as accounts, registrars 
office, departmental office etc. at IIT, Kanpur. lognoimal distri- 
bution is foimd to satis:fy the chi square test at one percent level 
of significance. 


Indian Standard (53) specifies a characteristic live load 
2 

of 250 to 400 kg/m for office buildings. If Qj^. is the characteristic 
live load and and s^ are the parameters of the live load following 
lognormal distribution 


P(Q >01;)= t 


108 ( 1 ^^,) 

s 


(2.27) 


in which p^ is the probability of the load greater than Qj^ . 
Parameters and s^ are given by (48) 

Sq = log + l) (2.28) 

^ exp (-0.5 s^) (2.29) 

vdiere 6 and Q are the coefficient of variation and mean value 
q m 

of the live load respectively, Oieir values are 0,381 and 
2 

130,3 kg/m . Using "Eas. 2,28 and 2.29, Qjj^ and are calculated 
2 

as 121.7 kg/m and 0.368 respectively. Substituting the above 

2 

values in Bq, 2,27 and assuming Qj^. is equal to 400 kg/m , the 


value of Pg^ is given by 
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log (400/121.7) ^ 

psl-4' [ ] = 6 X 10 

^ 0.368 

2.5 STATISnCAl MAOJSIS OP GEOMETEIC IROPERTIES OP A SECTION 

A number of prestressed concrete beams of the section shown 
in Pig. 2.19 have been cast and tested in the stmictural engineering 
laboratDiy at IIT, Kanpur for a research project. Measurements 
have been made to the accuracy of 0, 02jnmon the geometric properties 
viz, width of flange at top (b^), thickness of flange (t), thickness 
of web (b’)j total depth (h), depths of prestressing steels 
^*^1*^2*%^ of the section. IHie histogram and cumulative distribution 
of the data of bf, b‘, d^ and h are given in Pigs, 2,20 to 2,24. 
The normal distribution has been found to satisfy the chi square test 
at one percent level of significance in all the cases. It has been 
found that the coefficient of variation of b^, b’ and h are negligible 
and significant variation has been observed only in the thickness of 
flange and the depth of cable d^. Specified value, mean value, 
standard deviation, coefficient of variataon and the ratio of ihe 
mean value to its specified value of an paremeters of the section 
are given in table 2.6. Por the purpose of reliability analysis 
and design of beams, the variation of thickness of flange at top, 
t^, is taken as the average variation of t^ aind t^ and the variation 
of effective depth, d, as the average variation of and d^. 

Por considering the probabilistic variation of area of steel 
in the reliability analysis and design of prestressed concrete beams, 




1.8 204.2 206.6 209.0 

Width of Top Flange (m 
(b) CUMULATIVE DISTRIBU' 
'ARIATION OF WIDTH OF TO 
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robabtlity 









ative Freq 










204.0 


205.6 207.2 


208.8 210.4 


Total Depth (mm) 

(b) CUMULATIVE DISTRIBUTION 
FIG.2.24 VARIATION OF TOTAL DEPTH OF THE 
SECTION SHOWN IN FIG. 2.19 



59 


Table 2.6. Values of specified, mean, standard deviation, and the 
ratio of mean value to its specified value of geometric 
properties of the section of Pig, 2.19. 


SI. 

No. 

para- 

meter 

Specified 

Value 

cm 

Mean 

Value 

cm 

Standard 

Deviation 

cm 

Coeff. - 
of vari- 
ation 

Mean Value 
Specified Value 

1. 


20 

20.63 

0.164 

0.00794 

1.031 

2. 

b‘ 

14 

13.93 

0.128 

0.00918 

0.995 

3. 


2.5 

2.783 

0.152 

0.0546 

1.113 

4. 

^2 

4 

4.126 

0.193 

0.0468 

1.031 

5. 


2 

2.11 

0.181 

0.0858 

1.055 

6. 

^2 

13 

15.3 

0.244 

0.0159 

1.176 

7. 


16.3 

16.59 

0.226 

0.0136 

1.017 

8. 

h 

20 

20.62 

0.135 

0.0065 

1.031 
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the probability distribu'jtion of diameter of steel, is required, 
lig, 2»25 represents the histograa and cumulative distribution of 
the diameter of Tmm «|» HIS wire, Noimal distribution is found to 
satisfy i±ie chi square test at one percent level of significance, 
Bhe coefficient of variation is observed to be very anall. 
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LIABILITY ‘iJJALYSIS OF PEESTRESSED CONCSETE BEAMS 
AT STRENGTH LIMIT STATE UNDER DESIGN LOAD 

3.1 GEMERAL 

The probability distribations which have been determined 
in Chapter 2 for the field data on the strengths of M550 concrete 
and high tensile steel (7inm(fr), load and geometric properties of 
the section, are used in this chapter for the analysis of probability 
of failure of prestressed concrete (PSC) beams. A method of 
formulation for the reliability analysis of PSC beams at limit of 
state of strength under design load, using Monte Carlo techniq^ue, 
is presented for deterainistic and probabilistic loads considering 
probabilistic variations (PV) of 

(i) strengths of concrete and steel and 

(ii) strengths of concrete and steel and geometric properties 
of the section* 

3.2 INTRODUCTION 

Generally PSC structures have symmetrical or unsymmetnoal I 
or T sections. Putoha Chandrasekar (44) analysed the motions 
based on the governing equation for the ultimate resisting moment, 
of a section fixed by the deterministic analysis and assumed 
normal distribution for the resistance of the section. However, 
in I 8-nd T sections, because of the random variations of the 
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parameters of the resistance of a section, the ikilure can take 
place with anyone of the following eventss 

- the section is under-reinforced with the neutral axis in 
the flange 

- the section is under-reinforced with the neutral axis in 
the weh 

Y^ - the section is over-reinforced with the neutral axis in 
the flange 

Y. - the section is over-reinforced with the neutral axis in the weh. 
4 

The above events are assumed as mutually exclusive. The occurrence 
of each event has a certain probability. The probability tree for 
failure of a section at limit state of etrenglh under design load is 
gi'ven in Pig. 3 ‘I* If can be seen that the probability of failure, 

p^, of a section is the sum of the conditional probabilities of 
failures of the section under each given event and the same can be 
exnressed as 

4 

p^= J p(p!y.) P(Y.) (5.1) 

i=1 

where P(p|y^) denotes the conditional probability of P for given 
event Y^. P denotes the event ’failure* . P(Y^) is read as the 

probability of the event Y^. The conditional probability of 
failure of a section for any given event (say Y^) is given by 

P(PjY^) « P [(R-S) < 0 |y^] (3.2a) 

or, P(p1y^) » P[(| < 1)1 Yj^l 


(5.2b) 





FIG.3,2 TYPICAL PRCSTRESSED CONCRETE 
SECTION (IDEALISED) 
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where E is the resistance of the section and S is the action 
(load or bending moment) on the section. The resistance of the 

section is a function of various material and geometric properties 
of the section, 

R = f(x^j Xgj,*,,, Xj^) (3*3) 

Because parameters X. are usually random variables, the resistance 

3 

is also a random variable with density function fj^ and cumulative 

distribution Assuming the X. in E>i. 3.3 are statistically 

E j 

independent, their joint density function is 

n 

f(x^, Xg, X^) = n 

tl 3 

and its cumulatire piobabilily is 

n 

P„(r) = p(r < r) = / .♦* / n fy (x. ) dx (3.5) 

E - Q 3t=1 3 3 

The restriction R ^r define# the region of integration G in 
Bq^, 3»5» The integral contained in Bq, 3,5 cannot be evaluated in 
closed form. Added to this, the evaluation of p^ requires the 
evaluation of the probabilily of the occurrence of each given 
event A^. Defining 

= the event that the section is under-reinforced 
Bg = the event that the section is over-reinforced 

B_ = the event that the neutral axis lies in the flange 

3 

B, = the event that the neutral axis lies in the web 

4 
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the probabilily of the occurrence of the event is 

P(Y^) = p(B^n B^) (5.6) 

The events and. B^ are dependent on each other and the densily 
function of each is again the function of mteilal and geometrio 
properties of the section. Hence the evaluation of Sj, 3# 6 is 
again difficult. Finally, to calculate the conditional probability 
of failure for the given event, Bi.3.2a or 3.2b is to be used which 
involves numeiical integration. fhe evaluation of structural 
safety can thus become a fonnidable ta^, even when adequate statis- 
tical data are available. 

3.3 MONIE CAEOiO EECffiSlIQUE 

One of the usual obj ectives in using Ihe Monte Carlo 
simulation technique is to estimate certain parameters and proba- 
bility distributions of random variables whose values depend on the 
interactions with random variables whose probability distributions 
are specified. As it is known that the ultimate resisting moment 
of a section is a function of several random variables, the 
probability distribution of depends on the equation connecting 
these random variables. As explained in the previous section, 
as closed form solution to calculate cumulative probability of 
IS not possible, Monte Carlo method has been used in this thesis. 
Secondly, as seen in the previous section, the failure of a flanged 
section can take place under different events. Hence to study 
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and siuulate the oomplete random behaviour of the section at limit state 
of strength^ Ifonte Carlo teGhni(iufi is the beat suited method and 
used in this thesis to obtain the probability distribution and 
paraneters of IMs method consists of the following three 

steps (l6): 


1 ) Generating a set of values for the material properties 

and geometric parameters X. in accordanoe with the empiri- 

3 

cally deteiaoined or assumed density functions f^ f 

3 

2 ) Calculating the value of r^ corresponding to Ihe set of 
values x., obtained in step 1, by means of the appropriate 
resistance equation 3 •3* 

3 ) Repeating steps 1 and 2 to obtain a large sample values of R 
and therefore estimating 


Die parameters for drawing kth set of input values x^j^. from 
the corresponding density functions f^ is to generate first a 
sequence of n random numbers r with unifoim density in the range 
0 i.r ^1.0« These generated random numbers are then transformed to 
their corresponding values of Ihe particular distribution with its 
parameters, Humber of methods are available to generate uniformly 
distributed random numbers (54), 


2 

Assuming a and a are distributed as h( 422,8,56) kg/cm 
0 s 

and N( 15680, 488) kg/cm respectively, samples have been generated for 
M^, using Bi,3,f^, of 1he section shown in Fig, 3,2, by Monto Carlo 
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meliiod. denotes noimal distribution. vd.th parameters X.j 

and X^* rigs, 3*3a, 3«3b and 3*5c sOaow the plot of generated 
random samples of and 

3,4 SMILE SIZE 

lEhe generated data is used for estimating mean and standard 
deviation of the resistance of a section, Aa larger and larger 
samples are used, 1he estimates are closer to the population -values, 
The miniiEum size of the sample depends on the desired accuracy of 
estimates. 


Itor the es-timate of the population mean of a randan -variabl 
X, the mi ninum sample size is specified (54) such -that the probabil 
of the true mean falling within the confidence interval 

is (l-o) percent where X^ and s^ are sample mean and standard 
de-viation of X . If the length of the confidence interval, i, 
acceptable error in the es-feimate of mean value of X , xs specifiec 
as 


s 



then mininum size of sample for the estimate of the population 
of X is gi-ven by 



(^f 

m 


n 


(5,9) 
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5br large aaniple size (say n, > 120 ), the standard deviation of 

e 

standard deviation of x is equal to — Hence for the estimte of 

vSn’ 

standard deviation of Xy the minijiu'm size is specified such that the 
prohahilily of true etandard deviation falling within the -confidence 
interval (55) 

"x i Va 

is (l-o) percent. Speci:^ingi 

'5.11) 

the mininum sample size for the estimate of the population standard 
deviation of X is given by 

^ = 1 ^a/2 ( 5 . 12 ) 

s 

where e^ Is the acceptable error in the estimate of standard 
deviation of X. 

In all examples in this thesis more samples, than required by 
Eqs. 5.9 and 3.12 for a = 5 ^ and acceptable error of 5 percent in 
the estimates of mean and standard deviations, are generated. 

3.5 iPUmOHS IDE BEMMINATION OF ULUMATE STREEGOH OE A PSC SECTION 

It is assumed that the beams are flally bonded and the stress- 
strain cftirve fbr steel is as shown in ELg. 3.4. Equations, which 
are necessaiy fbr the Monte Carlo method, to determine for each 


event and Y^ are given b^ow. 
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Event Y, 


Ilie ultimate' resistxng moment of the section, '1'^® 

event in which the section is under-reinforced and the neutral 
axLS (EA) lies in the flange (ELg, 5»5a) is given hy (46) 


M ,► = A. d (l - ^ 

ruf s ts b^d 


(3.13) 


where A is the area of prestress steel. The section is under- 

uS 

reinforced if 


t cu 


(3.14) 


Event Y2(Eig. 3.5b) 

The ultimate resisting moment of the section, M_, , for the 

^ lUW' 

event Y^ in which Ihe section is under-reinfbiced and the neutral 
axis lies in the web, is given (46 ) 

(to^-b') t.(d-0.5 t. ) 
riOT ou t T X 


4. A rl ^tsw^s V 

°s \sw ” b'd a ^ 

CU 


(3.15) 


where 


A. = A . A . ^ 

tsw ts tsf 


(3.16) 


Q 


(3.17) 


The Eq,. 3. 15 is valid when 


A. a 
tsw 3 

b ‘d a 


<0.24 


(3.18) 
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Event Y_ 

3 

Eor the event Y^, i. e« the neutral axis lies in the flange and 

the section is over-roinforced, the stress in steel, a , at the 

' ' su' 

time of failure of tho section is given by 




( 3 . 19 : 


where 


V _ 

1 0„015 


(0,006 + e - £ - £ ) 

' c sp ce"^ 


^ “ o7oin: " s) 

ts 


in which is the limiting compressive strain in concrete and 

is the strain in prestressing steel and is the strain in 
concrete dixo to effective prestress. Assuming the d^^pth of stress 
block is equal to 0,8 times the depth of neutral axis (45), the 
depth of stress block, a, is calculated from the formula (56 ) 


a = 0,8d (■ 


^0 Sp” "^ce 


( 3 . 20 ) 


where e is the strain in steel at the time of failure of the 
s 

section. Assuming the stress-strain curve for steel shovi in 


Jig. 3e4, the value of e is calculated from the equation 


t = 0.006 (o - 0.8 (3.21) 

S 

Bie ultimate resisting moment of the section, for the event 

Y_ is given by 
3 
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M - = b . a Cd - 0.5a) (3^2a) 

Event 

Ibr- tbe event 'tiae n6ut]>al axis- lias in the ilaisge and 

the seotion ia oveav-xeinibKsed. tha value of o ia given by 

Oli 


whare 


= T^ [0.68 0^^ (b^- b') t^. ] 


2 A 


ts 

o 


® (0,006 + e -e -e^) 

'' ' o sp ce' 


0.015 

0.68 cf o 

- oToiT^ V°-°***o-Sp-'ce) J 

TS 


After evaluating the determination of and ’a' are same as 

explained for event Yy Values of ’a» and are given by B(is.5.20 
and 3.21 respectively. !I3ie ultimate resisting moment of the section 
M . for the event Y . , is given hy 

TOW 4 ^ 


M =0,68 0 [(b.-b*)t4.(d-0,5t. )tb*a(d-0.5a)] (3*24) 

row ou t t V 


In eveiy casej it can be seen that the resistance is a function 
of strengths of conocbete and steel and geometric properties of the 


section. 


5.6 HIOBABILITT ECSTRIHJELON OP M OP A SEOHON EOR PROBABIUSTIC 
VARIATIOJrS OP SiffiESrOTHS OP MATERIALS 

Even thou^ all parameters in the equations for deteimination 
of of a section are random variables, 


and a are 
s 
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first considered and the reliability analysis presented. Ihe typical 
prestressed concrete section shown in Eig, 5.2 has been used for the 

purpose of study of tbs probabilily distribution of Mth 

2 2 
strengths of concrete isr(422.8, 56) kg/cm and steel 11(1 5680, 488) kg/cm. 

llonte Carlo technique b^^s been used and samples have been generated 

■''<5ing Bjo, 5. 13 to 3.21- for of the section. Different sample 

^izes have been considered from 500 to 10,000. Olie effect of sample 

"ize and minifum sample size required for the determination of 

parameters of a random T,-anable have already been given in article 3,4^ 

iig. 3.6 shows the histogram of the generated data for of the 

section shown jn iig, 3.2 for a sample size of 10,000. It is seen 

that the distribution is negatively ^ewed, Bbimal, beta, lype I 

^cctremal (smallest) and Q^pe III extremal (smallest) distributions 

have mt been found to satisfy the chi square test for the generated 

data. The samples, fer which the section is under-reinforced and 

cver-reinforcod have been arranged separately and histograms and * 

cumulative distributions of the data are shown in ilgs. 3,7 and 3,8. 

With the deterministic analysis the section Is under-roinforced 


A, a 
f tsw s 

b "do 


O Q 

= 0,2) for o = 350 kg/cm and cr = 15000 kg/cm . It is 


cbseived from Ihg, 3,8a that about 10 percent of samples lie in the 
ovei^relnforced case. Hence there is a probability of the section 
>'ecoming over-reinforced even thou^ the section is under-reinfo iced 


'^“tPiminist1''‘'Tlye It is found that the nomal distribution. 


satisfies the generated data (tinder-reinforced case) shown in 
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at one percent level of significance^ However, rioimal 
distribution does not satisfy the generated data (over-reinforced 
ease) shown in ilg, 5.8a and it is found that 3ype III esjctremal 
(smallest) distribution satisfies the same data at one percent level 
of signxficanoe. Hence normal and !]ype III ertemial (smallest) 
distnhutions are adopted for ultimate resisting moments of the 
section for under-reinforced and over-reinfo3?ced cases respectively, 

OIhe probability density function of the normal distribution of 
the ultimate resisting moment of the section for under-reinforced 
case, is given by 

% (V=7=^ (5.25) 

ru /2ir s ru 

ru 

vdiere M and s ar« parameters mean and standard deviation 
rum TO 

of the resisting moment of the section for under-reinforced case. 

The cumulative distribution of M is written as 

TO 

M - M 

The probability density function of the l^pe III extremal 
(smallest) distribution of the ultimate resisting moment of 
the section for the over-reinforced case, is given by (48), 


- (tjr \ ^ ^ro / ^ro . r / ^ro \^ro 

ro rom rom rom U > o 

TO “ 


( 3 - 27 ) 


vdiere ymd are paraneters of the distribution, of 
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Hie parameters are estimated (48) as follows. 


M 


M 


im 


rom r(l + s^J 


(3#2S)' 


"ro 


^(1 

TO 

ro 


- 1 


(3.29) 


where M and 6 are the mean and coefficient of variation of 
rom ro 

M^. Dae relationship of 5.29 is graphed in Hg, 3.9* 
permitting one to find the parameter s^ if 5^ is known &r the 
data. The cumulative distribution of is given by 


(M„). 1 -«cp 

ro rom 


(3.30) 


M >0 
ro — 


3.7 EELLABIIITr MifflCSIS BOEMULATIOSr OP A SECTIOir 


Having arrived at the distributions for the ultimate resisting 
moment of a section for under-reinforced and over-reinfojoed cases, 
the reliability analysis is formulated as given below. 

The probability of failure of a section is given by, 

P^ = p(p|tj) p(tj) + p(p|o) t(o) (3.51) 

where TJ and 0 represenb the events iJiat the section is under- 
reinforced and over-reinforced respectively, Par any given 
section, using Monte Carlo method, number of samples are generated 
and the sanple means and standard deviations of the generated 
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data for the resistance of the section for under-reinforced and over- 
reinforced cases are found out. iix)in sample means and standard 
deviations, the parameters of th§ distributions of and 

can be estimated. 3!he value of l(u ) has been taken as 

P(U) = ■— (3.32) 

where n^ is the number of sanples of under-reinforced case and 
n is the total rwmber of generated saaples. Hence, 

p(o) = 1 - p(u) (3.33) 

3.8 CXDMHJTATIOir OP POR DETERMINI SHC 10 AD 

Assuming the load is deterministic, the conditional probability 
of failure of the section, p , for the given event D, is given by 

J. I t4. 

Pf|u = ‘ “el'') 

where M is the external bending momenb. The hatched area in 
Pig. 5.10 is the area vmder integration for the evaluation of 
probabilities of failure for deteiministio loads. ibr normally 
distributed M , 

I'll' 

P(P|u)-Pf,^= U— - g -— ) (3.35) 

Similarly, the conditional probability of failure of the section, 
for the gi-ven event 0, using Hi, 3.30, is given by 

M s 

p(p1o) = p = 1 - exp l-i (3.36) 

* rcm 
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The eztertal moment can be TOrking moment, or ultimate moment, 


M . 

eu 


M = F M + F M 
eu q q g g 


(3.37) 


where, 


F 


P 


g 


M 


g 

M 

q 


= live load factor 
= dead load factor 
= working moment due to dead load 
= working moment due to live load 


The of the section is calculated using Hi» 3»31. 
The procedure is illustrated with an example. 


Example 3-1 

Siii 5 )!ly supported PSC beams, supporting a floor, are spaced 

at 5m» live load on the floor is 400 kg/m^. Total dead load (including 

self wei^t of the beam) on any intermediate beam is 1500 kg/m. 

2 

The beams are designed with concrete cube strength of 350 kg/cm 

2 

and steel strength of 15000 kg/cm . The cross section of the beam 

arid qable profile are shown in Fig. 3-11- The effective span and 

2 

area of steel are 10 mand 11,83cm respectively. The beams are 
designed assuming th^ are cast separately and the composite action 
of slab and beam is not considered. One of the intermediate beams 
only is considered and Ihe reliabilily analysis of the same at limit 
state of strength is illustrated, 

live load on any intermediate beam = 5 ^ 400 * 2000 kg/m. 

From Chapter 2, for M350 concrete, 
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°cum ~ ^22.8 kg/cm^ ; = 56 kg/cm^ 

- 15680 kg/cm^ s s = 4-88 kg/cm^ 
sm ^ s 

using Monte Carlo method, 10000 samples are generated it is found 
that number of samples for under-reinforced and over-reiiaforced 
cases are equal to 8570 and I63O respectively. Using Bi* 3*^32, 

p(u) = 0.837 ; P(o) = 0.163 


Mean and standard deviation of liie generated samples of and 
are calculated (this is considered as a part of 1iie Monte Carlo 
method). Their values are given below, 

Por the under-reinforced case, 

M = 90,57 tm ; s = 2.84- tm 
rum ru 

Por the over-reinferced case, 


M = 82,58 tan; s^ = 4.85tm 
rom ’ ro 


6 


4.85 


ro 82.58 


S= 0.0586 


Using I^rpe III extromal (smallest) distribution for the parameter 
s^ is obteiined as 21.6 from the PLg. 3,9 for the calculated value 
of equal to 0,0586, Using Bq, 3.28, the value of M ^^ is 
calculated and found to be equal to 84.8 tm. The value of M^ 
at the section is 43,75 ‘tai. The conditional probability of failure 
for the given event U, using Eq. 3.55» is evaluated as 




■f.’a 


.75 - 90.57 
2.84 


) < 10 ‘ 


r24 
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-'24 

Hereafter, values less than 10 will be considered and repoirted 
as zero. 


The conditional probability of failure for the given event ^0% 
using Bq[* 3.56, is calculated as 


Pflo' ■> -^'-<34.8 


43^75x21.6 


r''°] 


= 6.44 X 10 




using Bi, 3.31, the probabilily of failure of the section is 
obtained as 

Vf = (O) (0-837) + (6.44 lo"^ ) (O.I63) 

-7 

= 1,08 10 

3.9 EELIABIIiaY MABTSES OP CONTIMJOUS BEJMS 

Ihe critical sections lAhere failures are likely to occur and 
cause the total collapse of the beam are to be first detennined for 
•ttie analysis of continuous beams. Pixed ends and intermediate 
supports of the beam, where masanum negative bending moments occur, 
are the one set of critical sections. The positions of critical 
sections for positive bending moments depend on the difference 
between internal resisting moment and external moment at different 
sections. In the caseof probabilistic analysis, it depends on the 
difference between the mean values of the internal and external 
moments and the stand.ard deviation of the resisting moment. Bie 
Insisting moment of the section is approximately proportional to the 
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eccentnciiy of the cable. !I!he second older equation for external 
bending moment for gL-ven load is known. The standard deviation of 
the ultimate resisting moment is approximately equal to the standard 
deviation of the strength of steel (observed later). Assuming 
normal distribution for (r-S), rt is possible to locate approximately 
the position of the critical section in each span within a few trials. 


After locating critical sections, Monte Carlo method is used 
and number of samples are generated for the resistance of each 
critical section. 33ie failure probability of each critical section 
for each loading case is calculated as explained in article 5 . 6 . 

In a redundant structure, it is known that the failure of the 
structure takes place only after sufficient number of sections have 
failed. In a continuous beam each ^an has a failure mode under a 
given load and each critical section has a certain probabilily of 
failure. Hence if there are E critical sections in the ith failure 
mode fcr given load the probability of occurrence of the failure 
mode, is given by the probability of intersection of the events 

Z^f defining the failure of each critical section. This can be 
written as 


p(Mbde i occurs under load 3 ) 

- n '^2 ^ 

Assuming events are independent. 


(^. 38 ) 



K 

k =1 ^^k 


(3.39) 
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where is the probabilily of failure of the section This 

can be evaluated as explained in article 3.8, If p^^ is the 

j 

probabilxiy of failure of the system under the load then bounds 
on p^g can be defined by (26) 


m 

1 Pfij 


(3.40) 


■where m is the number of failure modes, 
load conditions and mode resistances and p^^ 
to 1,0, bounds on p^^ are given by (9,26), 

m n 


Assuming independent 
is small compared 

3 


max 


Pfij I Pfij 


(3.41 ) 


in which p^^ is the probability of failure of the ^stem for m 
modes and n load conditions. The procedtire is illustrated with 
exampl es. 


Example 3.2 


OontinuDus PSO beams, having two equal spans and supporting 

a floor, are spaced at 5m« live load on the floor is 40Ckg/m^, 

Total dead Ipad (ino lading self wei^t of the beam) on ary- 

intermediate bean is 1500 kg/m. The beams are designed witb. 

concrete cube strength of 350 kg/cm and steel strength of 15000 kg/cm. . 

The beam has umiform cross section vifaich is shown in Eig. 3.12, 

2 

The effective span and area of steel are lOmand I0,89cm . Die 
beams are designed assuming that are cast separately euid 

composite action of slab and beam is not considered. 
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One of the inteimediate beams is only considered and the 
reliability analysis of the same at limit state of strength is 
illustrated, 

Idve load on ajqy intemediate b^m = 2000 ig/m. 

There are two possible loading coMitions for the live load 
(i) live load on both spans 
‘ (li) live load on one span only 
Ttor the case (i) the critical sections have been found to be 
(4,2mfrom the end A), C and B^(4.2mfrom the end E) of lig. 3.12 
and for the case (ii) (5ia from the end a), C and D2(5mfrom the 
end e). 

It is seen from Chapter 2 that strengths of M350 concrete and 

2 

Tram. <l> hi^ tensile steel are distributed as ir(422,8,56) kg/cm and 
w( 15680,488) kg/csm respectively. Using Ifonte Carlo method, 10000 
samples are generated for the resistance of critical sections B^, 

B^ and C of the beam* Erom the generated data parameters of the 

% 

ultimate resisting momenta of the critical sections for the events 
U and 0 are fo\ind out as illustrated in Example 3.1. Ihe results 
are given in table 3«1. 

load condition (i) : Both spans loaded: 

Magnitudes of external bending moments (due to live load and 
dead load) at the critical sections for different loading conditions 
are given in table 3«2. Erom the above table, 
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Talale 5 ‘I* Parameters of the iiltimate resisting moments of 

critical sections of the beam of Pig. 5-12 for 

FV of cr and a 
cu s 


SI. 

No. 

Section 

NR - n(m , s ) 
rum ru 

OR - EX. „ (M 

III , s rom 

, s )* 
ro"^ 

P(U) 

M 

rum 

(tm) 

s 

ni 

(tm) 

^0) 

■PH 

®ro 

1 


0.8787 

71.17 

2.2 

0.1215 

66,95 

21.8 

2 

% 

0.8718 

69 .48 

2.15 

0.1282 

65.40 

22.2 

5 

G 

0,9022 

75.66 

2.56 

0.0978 

71.06 

22.2 


*EXj jy g^^rom^^ro^ denotes Type III extremal (smallest) distribution 

with parameters M and s . 

rom ro 


Table 5»2» Magnitudes of external bending moments at critical 
sections of the beam of Pig. 5 -12 for deterministic 
load , 


Total External B.M. (tm) at 


., No. 

Live Loading on 
Spans 

B* or 

c 

D* or D 

1 2 

1 * 

AC & CE 

24.26 

43.75 

24.26 

2 

AC 

31.25 

28.15 

5.13 

3 

GE 

5.13 

28.15 

31 .25 


^Sections and coircespond to loading case 1 . 
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= 24.26 -tm ; = 43.75 tm 

where exterml moments at sections and C 

respectively. 


Probability of failure of section B^(or ) 

Using Bi. 3.35, the conditional probability of failure of the 
section B^ for the given event U is 

I’flu' 2.2 


, 24.26 - 71.17 ) 


Using Eq, 3.36, the conditional probability of failure of the 
section B^ for the given event 0 is 

Pf|o = 1 “ 1 - 2.44 X 10-10 

Die probability of failure of the section B.|, p^ ^ using Bq. 3.31, 
is given 'ty 

Pj.3 = (0.8787) (0) + (0.1213) (2.44 X to“^°) 

= 2.96 X 10^^^ 

Probability of failure of section C 


Using Bq. 3.35, the conditional probability of failure Of the 
section C for the given event U is 


,, _ ^ / 4?.75 - 75.66 ^ _ 

Pflu " 2.56 ^ ° 


Using Bq. 3.36, the conditional probability of failure of the 


section 0 for given event 0 is 
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pfro = 1 f = 2-1 « '<>■' 

using Bj, 3.31, the probabilily of failure of the section C is 


p„ = (0.9022) (O) + (0.0978) (2.1 x 10~^) 


■fC 


= 2.05 X 10 


Sections and C should fail simltaneously for the occurrence 
of failure mode in span AC, The probabilily of occurrence of the 
failure mode in span AC, using Bi. 3.39, is 


-17 


■fAC 


^fB. ^f C ^ ^ ^ ~ ^f CE 


1 


The maxunim value of the occurrence of failure mode for this load 

-17 

condition 1 is or P^qj; which is equal to 6,1 x -jo 

Defining p„ as the probabilily of failure of the beam for loading 

condition, (i), bounds on the value of p^^ , using Eq, 3.40, are 


6.1 X p. < 6.1 X + 6.1 X 10“"'’^ 

- fs^ 

6.1 X 10*'^'^:^Pfg ^12.2 X 10"’'*^ 

load Condition (ii): ^an AC loaded 
Erom table 3,2^ 


M 

“eBg = 31.25 tm ; MgQ = 28.13 tm 5 = 3. 13 tm. 

Probability of failure of section Bg 

Using Bis. 3.35 and 3.36, tiie values of conditional probability 
of failure of the seetion are 

Pfiu ' ° ' Pfio “ 
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iChe probabilily of failure of the section 

= (0.8718) (O) + (0.1282) (l.04 x io“®) 

=1.33 X io"^ 

Probability of failure of section C 

Using Bis, 3«35 and 3.36 » the values of conditional probability 
of failure of the section are 

The probability of failure of the section C is 

Pfp = (0,9022) (O) + (0.0978) (l.2 x io"®) 

-9 

=1.19 X 10 

Probability of failure of section 

Using Bis, 3.35 and 3.36, the values of conditional probability 
of failure of the section are 


fflu ' ° ' 10 ' ° 

The probability of failure of the section Dg is 
PfD " (0-9022) (O) + (0.0978) (o) = 0 
Die probability of occurrence of the failure mode in span AC, 
PfAC* ^ 


fAC PfBg ^fC 

= (1.33 xio"^) (1.19 xl0"^) = 1.58 xlO^"'® 
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Si nn larUy f "Hie probability of occurrence of "the failure node in span 
GE is evaluated and found to be equal to zero(i, e, less than 10 

It is found that the value of p^^^g is negligible compared to 

PfAC ^ Tipper bound and lower bound on the value of p„ are 

2 

same and equal to Hence 

^fs^ 1.58 X lo*"''® 

Similarly > the probability of failure of the beam, vfaea the live 
load occupies the span CE only is calculated and given in table 3.3. 

The bounds on the probability of failure of the beam, p^^, for all 
loading conditions and failure modes are calculated by using Eg,. 3.41 
and are given by 

6.10 X lo”'''^ _<p^g 12.52 X lo"''^ 

Example 3.3 

Continuous prestressed concrete beams, having three equal spans 

and supporting a floor, are spaced at 4m. live load on the floor 

is 400 kg/m . Total dead load (including self wei^t of the beam) 

on ary intermediate beam is 1000 kg/m. The beams are designed with 

concrete cube strength of 350 kg/cm and steel strength of 15000 kg^cm , 

The beams have unifbim cross section as shown, in Pig, 3.13. The 

2 

effective span and area of steel are 8m and 6,94 cm respectively, 
laie beams are designed assuming 1ii€y are cast separately and the 
couqwsite action of slab and beam are not considered. 



TaMe 3*5* Bounds on the PF of the heam of Fig* 3 *12 for deteiministic load and F7 










98 


One of the inteimediate beams only is considered and the reliabililry 
analysis of the same at limit state of strength is illustrated, 
live load on aqy inteimediate beam = 4 x 400 = 1600 kg/m. 

Gliere are five possible loading conditions for the live load 
acting on the beam as stated below; 

(i) love load on all spans 

(ii) live load on end spans AC and EG- 
(lii) live load -on middle span CE 

(iv) live load on end span AC(or EG) only and 

(v) Uve load on two adjacent spans. 

!Che critical sections are foxind to be (5.36m from the end a), 

C,D,E and (3.36m from the end g) of PLg, 3.13 for cases (i) and 
(v) and Bg (4m from the end, A), C,D,E and Pg (4m from the end g) 
for cases (ii), (iii) and (iv). 

It is known from Chapter 2 that strengths of M350' concrete and 

2 

7mm 4 high tensile steel are distributed as n( 422.8, 56) kg/cm and 

2 

Il( 15680, 488) kg/cm respectively. Using Monte Carlo technique, 
parameters of the ultimate 3?esistlng moments of critical sections 
for events U and 0 are obtained and given in table 3.4. 

Magnitudes of the eocternal bendirg moments (due to dead load 
and live load) at -Uie critical sections for different loading cases 
‘are shown in table 3»5* Pailure probabilities of the beam are 
calculated for different loading cases. 
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Tatle 5»4« Parameters of the ultimate resisting moments of 

critical sections of the heam of Pig, 3»13 

PP of a „ and a 
cu s 


SI. 

No, 

Section 

TJR - N(M , s ) 

' rum ’ ru'' 

OR - EXt^t (M 

lIIjS ^ rom 

’ ®ro^ 

I(TJ) 

M 

rum 

(tm) 

^ru 

(tm) 

P(0) 

M 

rom 

(tm) 

^ro 

1 


0.8529 

34.41 

1.06 

0,1671 

32.35 

21.3 

2 


0.8175 

35.50 

1.04 

0.1825 

51.55 

21.0 

3 

0 

0.8757 

57.60 

1.16 

0.1243 

35.35 

22.2 

4 

D 

O.6O72 

24.30 

0.75 

0.5928 

22.78 

18.7 


Table 3»5* Magnitudes of external bending moments at critical 
sections of the beam of Pig. 3 *15 deterministic 
load 


SI. 

No. 

Live Load 
on Spans 

Total External B.M. (tm) at Section 

or 

c 

D 

E 

P* or Pg 

1* 

AC,CE & EG 

15.28 

16.64 

4.16 

16.64 

13.28 

2 

AC & EG 

15.04 

11.52 

5.52 

11.52 

15.04 

5 

CE 

2.24 

11.52 

9.27 

11.52 

2.24 

4 

AC 

14.16 

13.26 

0.116 

4.66 

5.67 

5* 

AO & CE 

12.55 

18.58 

6.72 

9.^78 

1.03 

6 

EG 

5.67 

4.66 

0.116 

15.26 

14.16 

7* 

OB & EG 

1.05 

9.78 

6.72 

18.38 

12.55 


^Sections Ej and P^ correspond to loading oases 1, 5 and 7* 


I 
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Load Condition, (i): live load on all spans 

Etom table 3.5, the values of bending moments at the critical 
sections are 


'5.26 t. , 

M „ = M _ = 16.64 tm, ^ M = 4. 16 tm. 
eO eh eD 

Ihe conditional probabilities of failure of the critical sections 
C, D, E and for events U and 0 are calculated as illustrated 
in the previous examples and are given in table 3.6. The probability 
of failure of the above sections, calculated by using Bi, 3.31 are 


PfB “ ^fD " 

1 1 

p^Q = 6.8 X 10“^ ; = 6.1 X lo"''^ 


There are 3 failure modes in three span continuous beams. Sections 
aiyi C should fail for the occurrence of failure mode in the span 
AC. Its probabilily is given by 

PfAC = PfB^ PfC = = PfEG 

Sections G,D and E should fail for the occuirerce of failure mode 
in the span CE, Its probabilily of failure is 


^fCE ~ ^f C * ^fD ' Pf E 
= 2.8 X 10"^^ 

Bounds on the value of p^g , using Bq,. 3.34, are 

6.6 X 10“^®^p^g <13.2 X lO"^® 
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as 
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Sunllarly liie reliabilily analysis is done ftn* othea? Load conditioias 
and -aie bounds on the probabili-ty of failure of the beam for different 
load conditions are given in table 3,7, 

Prom the above table, the bounds on the p. for all loading 

iS 

conditions and failure modes are eq^ual to: 

1,8 10 ^ p^g ^5»0l X 10 

3,10 IROBABIIiar DISTRIBUTION OP M OP A SECTION POE EROBABIIISTIC 
VARIATIONS -OP Og AND DIKENSIONS OP SECTION 

In previous article, the probabilistic variations of atreogths 
•of materials only aere coosadered in tha resistance of the section} 
but in practice, as observed in CShapter 2, there are random discre- 
pancies between actual and nominal cross-rsectional dimensions of 
structural members, areas of steel bars and position of reinforcement. 
The differences between actual and nominal dimensions , expressed as 
percentage of the nominal dimensions, are generally not independent 
of these dimensions. Since lack of data on the variations of these 
differences with nominal dimensions for Indian conditions, 'ttie ratios 
of the mean dimensions to their respective characteristic values 
and coefficients of variation of parameters established for the 
section of Pig. 2,19 in Chapter 2, are assumed same for all sections 
in this thesis, 

I^ndom variations of Ihe variables h^fb', t^, d and D^ in 
the equations ^ven in the article 3,5, are considered in the 
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resistance of the section,. Usi33g the ratios of the mean values 
to iiieir respective oharacteixstic values of b*, t,j. and d 
and their coefficients of variation, the mean values ant! standard 
deviations of b^, b’, t^ and d of the section shorn in 5ig, 3.2 
can be oalculated, IBhese values are given in table 3.9, Using 
normal distribution (observed in Chapter 2) for a.n random varLablee 
b^, b’, t^|d, °a "ttieir parameters, Monte Carlo 

simulation method is used and the distribution of the resistance of 
the section is obtained. The histogram and cumulative distribution 
of of the section for lOOOO samples are given in Hg. 3.14. 
Normal, beta, ilSrpe I extremal and 3ype III extremal distributions 
have been tried for the generated data and none of thaa is found to 
satisfy the chi square test, Bie samples for over-reinforced and 
under-reinforced cases are arranged separately and their histograms 
and cumulative distributions are shown in Jigs, 3.15 and 3.16, It 
is observed that about 4.4 percent of the samples lie in the over- 
reinforoed case even though the beam is under-reinforced determinis- 
tically. It is fbund that the normal and lype III extimal 
(smallest) distributions satis:^ the ohi square test of the 
generated data of Jigs. 3*15 and 3.16 respectively for one percent 
level of slgrdficanoe. Hence the normal distribution for under- 
reinforced case and iiype III extremal (smallest) distribution for 
over-reinforced case are aidopted in the reliability analysis of 
the PSO beams at limit state of strength. The formulation of the 
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the reliahilily aaaly*le at limit state of atpength is same as 
explained in article 3.7. !I!he procedvire is illustrated with 
examples.. 

Example 3.4 

Ghe sin 5 >ly supported ISC beam of Example 3,1 (fig. 3, 11 ) is 
considered for probabilistio variations of strengths of materials 
and geometric properties of the section. Other particulars of 
the beam remain same as in Example 3,1, 

Using 1he coefficients of variation and the ratios of the mean 
values of the paraneters b^ b', t^, h and d to their respective 
characteristic values (found in Chapter 2), 1he mean values ard 
standard deviations of b^, b',t^and d of the beam shown in 
Eig. 3,11 are calculated as fellows: 

b. = 1.031 *< 50 = 51.55cm 

xm 

b* = 0,995 *< 15 = 14.925 <an 

m 

t^ = 1.072 X 12 = 13.936cm 

d„ = 1.083 54.7 = 59.24cm 

m 

h = 1.031 X 61 = 62.89 cm 

m 

s^^ = 0,00794 X b^= 0.409cm 

s, , = 0.00918 X iji = 0.137cm 
0 ' m 

®tt “ °*^507 = 0.7065cm 

= 0.0384 ^ = 2. 079 can 

8^ = 0.0065 = 0,409cm 
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It is known from CJiapter 2, 

2 2 
~ 4-22«S kg/cm ; s =; 56 kg/cm 

'Oum c 

0 = 15680 kg/can^ ; S => 488 kg/csn.^ 

S 6 

D„ = 7.065 mm ; = 0.0345 mm 

cxul U 

Using normal distribution for all the above random variables, 20,000 
samples are generated for applying Monte Carlo method. Hie 
number of sanqjles in events U and 0 are equal to 18600 and 1400 
respectively. Using Bq. 3.32, 

p(0) = 0.07 and l(u) * 0.93 

Bor tile under-reinfbrced case (normal) 

M = 100.51 tm ; s = 5.37 tm 
rum ru 

Bor the over-reinforoed case (l^pe III extremal smallest) 


“lom = 5 

Using Eqs. 3.35 and 3.36, smd found to be equal to 

zero and 2.0 10 respectively. The probability of failure of 


the beam is 


Pf = (0.93) (O) + (0.07) (2.0 X io"^) 


= 1.4 10 


(The method of reliability analysis of continuous BSC beams is 

same as explained in article 3.8 except that the parameters of the 

% 

resistance for events U and 0 are to be obtained takirg probabilistic 
variations of strengths of materials and gecaaetilo properties of tt^e 
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section as illustrated in the previous example. Ihe results of the 
reliability analysis of the two span beam of iig. 3-12 are given in 
tables 3*9 to 3.12 and for the three span beam of ELg. 3*13 in 
tables 3.13 to 3.16. 

3.11 COMIUfAlIOlT OP PROBABIUOY OP PAIHIEE POP HIOBABILISTIC 
VARIAPDIT OP LOAD 

Por probabilistic variation of action (load) and resistance, 
the probability of failure is given by 

00 » 

p = p(R < S) = / fJr) [ / f (s) ds ]dr (3.42) 

^ — 00 r ^ 

where f„(r) and fc,(s) are probability density fluactions of R 
and S respectively. 

The above Bq,. 3*42 can be rewritten as 

00 

p = p(r < s) = / f (r) [1 - Pg(r) ]dr (3.43) 

— 00 

inhere ^g(r) is the cumulative distribution of S. Hie area of 
the hatched portion in the iig. 3.17 gives the probability of 
failure. 

(a) Lognormal Distribution for Load and Normal Distribution for 
Resistance 

The probability distribution for the live load has been found 
to be lognormal in Chapter 2. Since the bending moment is directly 
proportional to load, the probability distribution of is also 
lognomal viiidi is given by 
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« 

Table 3.9» Mean values and standard deviations of 

dimensions of the sections of the beam of Pig. 3 •12 


SI. 

Ho. 

Item 

Mean Value 

Coefft. 

of 

Variation 

Charac- 

teristic 

Value 

(cm) 

Mean 

Value 

(cm) 

Standard 

Deviation 

(cm) 

Distri- 

bution 

Characteri- 
stic Value 

1 


1.031 

0.00794 

50 

51.55 

0.409 

Normal 

2 

b’ 

0.995 

O.OO9I8 

15 

14.93 

0.137 

Normal 

5 

*t 

1.072 

0.0507 

12 

13.94 

O.7O6 

Normal 

4 

\ 

1.083 

0.0384 

47,36 

51.29 

1.970 

Normal 

5 

\ 

1.083 

0.0384 

46,35 

50.20 

1.928 

Normal 

6 

dp 

1.083 

0.0384 

50 

54.15 

2.019 

Normal 

7 

h 

1.031 

0.0065 

56 

57.74 

0.375 

Normal 


Table 3»10. Parameters of the ultimate resisting moments of critic5al 

sections of the beam of Pig. 3*12 for PV of o , a and 

^ cn ^ s 

dimensions of section 


SI. 

Ho. 


1 

2 


Section 




TIR - N(M^,^ , s^J OE - E3 C_t- (M 9 ) 

rum ru' 111,8 rom* ro 


P(TI) 

M 

rum 

(tm) 

®ru 

(tm) 

P(0) 

M 

rom 

(tm) 

s 

ro 

0.9512 

79.14 

4.25 

0.0488 

74.70 

17.0 

0.9476 

77. 18 

4.23 

0.0524 

75.12 

16.0 

0,9613 

84.17 

4.59 

0.0387 

79.06 

17.7 


5 


C 
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Table 3*1 3 • Mean values and standard deviations of dimensions of the 
sections of the beam of Ehg, 3*13« 


SI. 

No. 

Item 

Mean Value 

Coefft. 

of 

Variation 

Charac- 

teristic 

Value 

(cm) 

Mean 

Value 

(cm) 

Standaid Distri- 
Deviation bation 

(cm) 

Character- 
istic value 

1 


1.031 

0.00794 

40.0 

41.24 

0.327 

Normal 

2 

b» 

0.995 

0.00918 

13.0 

12.94 

0.119 

Normal 

3 

t 

1.C72 , 

0.0507 

10.0 

10.72 

0.544 

Normal 

4 

% 

% 

1 .063 

0.0384 

36.05 

39.04 

1.499 

Normal 

5 

1.083 

0.0384 

35.25 

38,18 

1.466 

Normal 

6 


1.083 

0.0384 

39.0 

42.24 

•1.622 

Normal 

7 


1.083 

0.0384 

26.75 

28.97 

1.112 

Normal 

8 

h 

1.031 

0.0065 

44.0 

45.36 

0.295 

Normal 

Table 3»14« 

Parameters 
sections of 
dimensions 

of the ultimate resisting moments of 
the beam of Fig. 3*15 IV of 0 , 

of section 

crxtioal 

a and 
s 

SI. 

Section 

UR - N(M , s ) OR ■ 

^ nim ’ ra 

- EX j 

iiljS 1X321 

.S ) 
ro 

INO# 


P(TJ) 

M 

rum 

(tm) 

s P(0) 

ra ' ^ 

(tm) 

M 

rom 

(tm) 

®xo 

1 


0.9285 

38.26 

2.1 0.0715 

36.10 

16.9 

2 


0.9213 

37.28 

2.01 0.0787 

35.18 

17.0 

3 

C 

0.9503 

41 rS 0 

2.28 0.0497 

39.34 

17.4 

4 

D 

0.7558 

27.10 

1 .49 0*2462 

25.6O 

16.0 
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(M^) = ♦ [ 

e 


M 

log ( if. 


em 


( 3 . 44 ) 


me 


where 51, (M ) = owraulative distribution of M 
M e e 

e 

®*'em ^ ®me ~ parameters of following lognormal distribution, 

When the beam fails in tension (undeiwreinforced), it has been found 
in the article 3.5 that the probabilily distribution of the resistance 
of the section is normal. Hence using Hj., 5.25 for the probability 
density function of the resistance and substituting 5):i« 3.44 in 
Bq,. 3.43| the probability of failure is obtainei as 


Ph> = / 




ru 


, M - M „ 

r 1 / ru rum\2,,, 

exp [ - g ) ]x 

ru 


log (M 

n 

me 


(3.45) 


Substituting 


^ M - M 
1 / ru rum ^ 

w = — ’ V ■ ) 

^ V • 


(3.46) 


the Bq, 3.45 simplifies to 


P,=/“ e- f [, -«{ 

^ --00 /tT 


/2 w s + M 

logi SS ) 


M 


Wl 


}]dw ( 3 . 47 ) 


me 


Ihfe external bending momeiit (m ) consists of live load bending 

w 

momerct (m ) and dead load bending moment (M ): 

<1 g 
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M = M + M 
e g q 


(3.48) 


Assuming M as detennixiistic, the cumulative probabilily of M 
g e 

can be calculated from the known distribution of M . 

q 


P(M„ < ) = P(M + M < M ) 

e ru^ ^ g q 

= I(\ < - Mg) 


(3.49a) 

(3.49b) 


Hence 


e 


M - M 

Si g_) 

M ’ 

-jm. 


log(' 


(3.50) 


mq 


where M and s are parameters of live load bending moment 
qm mq 

(log-normal). As is proportional to Q, is equal to the 
dimensionaless parameter s^. Substituting the Eq, 3.50 in Bq, 3.43, 
the probabilily of failure is given by 


w 3 + M - M 

log ( 


^f - 


00 2 , ° ' M 

/ -w 1 ^ j- qm 

= y e — [1 - 


/iT 


■ } ] dw 

(3.51) 


(b) lognormal Distribution for load and !type III Ertremal 
(smallest) Distribution for Resistance 

When the beam fails in con^jression (over-reinforced), the 
probability distribution for ’the resistance of the section has been 
found to be !5rpe III extremal (smallest) in the article 3.5. 

Using Bi. 3.27 for the probability density function of the resistance, 
and following similar procedure of case (a), the probability of failure 


is obtained as 
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M w 

log 


(^t) 


ro 


Pf = 


/ e"^ [1.0 

o 


- « {- 


M 




32l. 


}] dw 


(3.52) 


Using numerLcal quadrature method, the integral equations 
3.51 and 3.52 can be evaluated. IDae above equations are in the 
standard tbrms and are evaluated in this thesis using Hennite-Gauss 
quadrature formula and Laguerre-Gecuss quadrature formula (57). 


idle parameters of Uie resistance of a section for events U 
and 0 are obtained using Monte Carlo technique. The conditional 
probability of failures Pf 1 o evaluated using Dqs. 3.51 

and 3.52, for probabilistic variation of live load. Shrther 
reliability analysis of sin^jly supported and continuous beams is 
same as explained in the articles 3,7 and 3.9. Hie procedure is 
illustrated witii examples. 

Example 3.5 

Ihe siii 5 )]y supported pSO beamof Example 3*1 (lig* 5.1l) is 
considered ftr probabilistic variations of live load and strengths 
of materials, AH particulars of the beam remain same as in 
Example 3.I. 

It is known from caaapter 2 that the live load is distributed 
as I(ir( 121.77 kg/m , 0. 368). The mode value of the live load 
bending moment is 

vs 121.77 ^5 >^ 10 ^ 


= 7610 kg m 
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Since the bending moment is proportional to the load, is distri- 
buted as (7»6l tm, 0,368). Ihe follomng results are taken from 
Ezan^jle 3.1. 


M 

g 

= 18.75 tm 



p(u) 

= 0.837 tm 

; p(o) 

= 0.163 

M 

rum 

= 90,57 tm 

ru 

= 2.84 tm 

M 

rom 

= 84.80 tan 

: s 
^ ro 

= 21.6 


Substituting the above values in Bjs. 3.51 and 3.52 and using numerical 
quadrature method (57 )» the yalues of p and p are obtained as 

i |U I I 0 

—10 -9 

6.86 ^ 10 and 6.8 x 10 respectively. Using Bi. 3.31, the 
of the section is 


p^ = (0.837) (6.86 X 10'*''°) + (O.I63) (6.8 x 10“^) 

= 1.682 X lO"^ 


Ihe reliability analysis of continuous beams for probabilistic 
variation of live load is same as e3q)lained in article 3.9 except 
that probabilistic variation of live load is to be considered and 
conditional probability of failures are evaluated using Bqs, 3.51 
and 3.52. Hie results of the reliability analysis of the two span 
beam of ilg, 3.12 in tables 3.17 to 3.19 and of the three span 
beam of Hg. 3.13 in tables 3.20 to 3.22 for probabilistic variations 
of Q , and o^ are piresented, 

Ihe method of reliability analysis of beams for probabilistic 
variations of live load, strengths of materials and geometric properties 
of the section, is illustrated vd-th an eacample. 



Table 5.17* Magnitades of external bending moments (mode' values of M ) at critical sections 
of the beam of Fig. 5.12 for probabilistic live load 
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Bounds on the probability of failure of the beam s 5*22(-20) ^ ^6.46(-20) 

4 - Sections B-, and correspond to loading case 1, 


Table 5,20. Magnitudes of external bending moments (mode values for M ) at critical sections 
of the beam of Fig. 3* 13 for probabilistic live load. 
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-(26 




5.2(-11) is read 
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Example 3*6 

The simply supported PSC beam of Example 3,1 (lig. 5*11 ) is 
considered and the reliabilily analysis of the same beam, fear 
probabilistic variations of live load, strengths of materials and 
geometric properties of the section, is illustrated. 

It is known from Ecample 5.5 that the live load is distributed 
as ^[(T.SI tm, 0.368), The value of M is equal to 18,75 tm, 

O 

The following results, obtained by using Monte Carlo method, are 
taken from Example 3.4« 


r(o) 

= 0,07 J 

P(U) = 0.93 

M 

= 100,51 tm; 

s =5.37 tm 

Tum 


ru 

M 

3X5m 

= 94.67 tm ; 

TO 


Substituting the above values in Eqs, 3.51 and 3.52, the values of 

-10 -9 

p„, and p„, are obtained as 1,16 ^ 10 and 1,08 x 10 
■^flu -^flo 

respectively. Using Bq. 3.31» the p^ of the section is 

p^ = (0,93) (l.l6 X 10"''°) + (O.O?) (1,08 X 10"^) 

= 1.84 10“^ ° 

The results of the reliabilily analysis of Mae two span beam 

of Fig, 3.12 are presented in tables 5.23 aiod 3.24 and of the three 

span beam of Pig, 3.13 in tables 3,25 and 3.26 for the probabilistic 

variations of Q, cf , o and dimensions of section . 

' cu'^ s 
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CHATTER 4 


EELIABmOS: AHAiySIS OP IRESIEESSED CONCRETE 
BEAMS AT LIMIT STATES OP CRACKING AND DEFLECTION 

4.1 GENERAL 

A structure is said to have failed if it has lost its 
resistance against load or has become unserviceable because of 

undesirable deformation or cracking etc. The limit states can 

V 

be placed in two categories: 

(a) the ultimate limit state whida coixesponds to Ihe maximmi 
load carrying capacity? 

(b) the serviceabiliV limit states which are related to the 
criteria governing normal use and durabilily. 

The causes leading to the attaiiment of serviceability limit 
states axe 

(i) excessive deformations with respect to normal use of structure 

(ii) prepiature or excessive cracking 

(iii) undesirable damage (oorrosion) 

(iv) excessive displacement without loss of equilibrium 

(v) excessive vibrations, etc. 

limit states of initial crack formation and deflection are 
considered in this thesis. Reliability analysis of sin^jly 
supported £md continuous prestressed concrete beams at limit 



132 


state of Initial crack formation is first presented for dete rmin istic 
and probabilistic loads considering probabilistic variations of: 

(i) strengths of concrete and steel 

(li) strengths of concrete and steel and geometric properties of 
the section. 

Reliability analysis at limit state of deflection is presented for 
the probabilistic variations of: 

(i) strengths of materials 

(ii) strengths of materials and live load 

(iii) strengths of materials, live load and geometric properties of 
the section. 

4,2 REUIBIHTS: AHAKSIS AT LIMIT STATE OP CRACKING 
4,2.1 Introduction 

Codes specify that prestressed concrete beams possess certain 
minimm resisting moment at the limit state of cracking which can be 
defined as the state when the maxinum tensile stress in the extreme 
fibre due to external moment reaujhes the value of modulus of rupture 
of concrete, TiTi a external moment producing first crack in a 
prestressed concrete beam is conqjuted by the elastic theory assuming 
that cracking starts when the modulus of repture of concrete is 
exceeded, Tte-p-i rrf ng M as the resisting moment of the section 

at initial crack fbimation, the beam is said to become unserviceable 
^en the external moment exceeds -yae value of 


Hence the 
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probability of cracking of the section, p^, can be written as 

p^=p(Mg>M^) (4.1) 

is a function of several random variables of strengths of 
materials and geometric properties of the section. The probability- 
distribution and parameters of are best obtained by Monte Carlo 
technique, explained in the pre-vious chapter. 

4.2.2 Equations for the Determination of M 

2X5 

The following equations, -bp determine M of a section, are 

^ rc 

necessary in the Monte Carlo me-bhod. It is assumed that Idle 

cables are fully bonded and so the trsnsfoimed area of cross 

section is -used in the calcula-feLon of M . !IIie state of stress 

pc 

and strain at limit state of cracking is shown in lig. 4.1. If 
Eg is the Young's Modulus of steel, is the tensile strain at 

initial (tracking, m is the modular ratio and is the modulus 
of rupture of concrete, the depth of neutral axis, x, is given by 

2r^+ TgCr^ + r^) 

nhera 

(b^- b') - (bjj- bO + 2ht^j(b^j-b») + 2rgd + b 'h^ 

^2 = 2 V"r 


r^ « 0.66 A^g 
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r, = E e . 

4 s ct 

rg = (m-1 ) 

denotes bi^adth of bottom flsinge and t^, the thickness of bottom 
flange. The 4.2 is valid whether the neutral axis lies in 
the web or bottom flange. 

Case 1. Neutral axis in the web: 

The foUowmg equations are used to determine M if the 
neutral axis lies in the web. 


The coopressive force in concrete, C t is 

° * * * 2(h-x) 

which acts at a distance 

5 b^X^- (b^-.b») (x-t^)^ 


(4.3a) 


(4.4a) 


from the r«utral axis. The tensile force in concrete, T^, is 
given by 

T= [b (h-x)^-(b^-b') (h-x-t. )^+ 2rg(d-x) 1-^^ (4.5a) 

° D D 2(h-x) 


which acts at a distance 

\3' 


-2-i ' 


2 r ' )(^-*-\)^+ 3rg(d-x)^ 


b^(h-x)^-(b^|-b’)(h“x-t^) + 2rg(d-x) 


] 


(4.6a) 


frcnn the neutral axis. Assuming e„_ is eqxial to one tenth of 

c© 

e and effective tensile stress in prestress steel to the maxinum 
00 
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value of 0.6 o (46), the strain in steel, e , at initial cracking 

S 3C 


is obtained, from the equation 


so 


0,66 CT 


( — ) 

ct h-x ^ 


(4.7) 


The total tensile force in steel, T^, is 
=■ 'sc ®s ^s 

whidi acts at a distance (d-x) from the neutral axis. Taking 

moments of C , T and T about neutral axis, M , is given by 
CCS ' rc' 


M = C ’ x. + T x„ + T„ (d-x) 
re cl c 2 s ' ' 


(4.9) 


Case 2. Neutral axis in bottom flange; 

The following equations are used to detemnine M if the 
neutral axis lies in the bottom flange. 


C = [ (b -b » )t . ( 2x-t . )+(b, -b ' ) (x-h+t, b 'x^ ] 

0 t t t b b 2 (h-x) 

■3 *2 P 9 (4«3b) 

2 r (V^')( 3 c-h+t^) +b'x\(b^-b’)t^( 3 x “ 3 xt.j+t^ ) 

2^ SS — 1^ — I. ■■■ III. . II— ..■■■II . I -I II l.l Pill ■■■ I— M il. M »-l II — ■■ ' ■■■ — J 

^ ^ (b^-b')(x-h+t^j)^+b«x^+(b^;-b')l:^(2x-t.^) 


T 5= fb. (h-x)^ + 2r^(d-x)] ^ — 

° ^ ^ 2(h-x) 


(4.4b ) 
(4.5b) 




Values of 


- b (h-x)^ + 3r^(d-x)^ 

I 5 ^ i (4.6b) 

b^(h-x) + 2rg(d-x) 

e , T and M are given by Bq.s, 4.7, 4*8 and 

SC S PC 


4.9 respectively 



137 


The above equations are fanctions of material properties 

^ct^ ®r» geometric properties b^, b^, t^, t^, b', 

d, h and A. of the section. As elastic theoiy is valid upto 
t»s 

initial crack fonnation (58) 
o 


^ct “ 


B 


( 4 . 10 ) 


where is the Young' Modulus of concrete which is given by (46) 


B = 18000 Aj ( 4.11) 

c cu ' 

The value of is equal to 

<^r = 2-25 ^ (4-12) 

m 

Hence o and E can be connected to the stwngth of concrete 
r c 

whose probability distribution and parameters are known. Aj.^ 
can be connected to diameter of steel whose probability distribution 
is also known. It is assumed that random variations of b^^ and 

t, are same as b . and t . . Hence probability distribution and 

U M U 

parameters of all random variables are known. 


4.2.3 Trobabilily Distribution of for Probabilistic Variations 
of Strengths of Materials 

It can be seen from Bq. 4.9 that is a function of several 

rc 

randcm vaidables. The analysis is first presented considering 
probabillstio variation of o_,^ anl wMoh are dlstritatei 

as h( 422.8,56) kg/cm and N( 15680, 488) kg^cm respectively 
(Refer Chapter 2). Using Monte Carlo method, explained in article 3.3, 
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10000 samples have been generated for M of the section shovm in 

rc 

Fig, 3,11. The histogram and cumulative distribution of the 
generated data are shown in Jig, 4.2. The skewness is veiy veiy 
!=!itia.n fliyl nomal distribution is found to satisfy the chi square 
test at one percent level of signtlfcance. Hence normal distribution 
is adopted for of section for probabilistic variations of 
and a • 

4,2,4 Computation of Probabilily of Cracking for Deterministic 
Load 

aihe probability of cracking of a section has been defined by 
Eq, 4,1. For normally distributed M , the probabilily of cracking 
is given by 

M - M 
rc 

vdiere M and s are the parameters mean and standard deviation 
rcm rc 

of The method of reliabilily analysis at limit state of 

cracking is illustrated with an example. 

Example 4,1 

SUrnply supported PSG beams, supporting a floor, are spaced at 

2 

5m, live load on the floor is 400 kg/m • Ibtal dead load 
(including self«-weight of the beam) on ary intermediate beam is 

1500 kg/m. The beams aia designed with concrete cube strength 

2 2 
of 350 kg/cm and steel strength of 15000 kg/cm . The effective 

2 

span and area of steel eire 10m and 11,83 on respectively, dhe 




52.5 5A.4 56.3 58.2 60.1 62.0 63.9 65.8, 

. . ' ■ Mrc (tm) - 

,(b) CUMULATIVE, OIStRlBUTJON 
FIG. 4.2, VARiATIOKl,QF Mre FQR PV OF cfeu AND oj 
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cross section of the beams and cable piofile ar« shown in ilg. 5*11 
which IS again given here for ready refe 3 ?enoe, 3he beams are 
designed ets sunning they are cast sepsirately and the composite action 
of slab and beam is not considered. 

One of the intermediate beams is considered and the reliabilily 
analysis of the same at limit state of cracking is illustrated, 

Eor this beam, using Monte Carlo method paraneters of M have 
already been obtained in article 4.23 and given in Fig. 4.2. Their 
values are 


M = 59,56 tm ; s = 1,64 tm 
ro ’ rc 

Die external bending moment at midspan, from Ecample 5*1> is 43.75 tm. 
Using Bi, 4.13, the probability of cracking of the beam is 


= H -- ) = 2,7 ^ 


1.64 


10 


-22 


4,2.5 Reliability Analysis of Continuous PSC Beams at limit State 
of Cracking 

Continuous beams are analysed for different possible live load 

conditions. Bor each loading case, the critical section where 

probability of failure is iuaxiiium, is determined. Bor deterministic 

load, the critical section depends upon the difference between the 

mean values of M and M and s of each section. The value 
ew ro rc 

of M of a section is approximately proportional to the eccent- 
rc 

ricity of cable at that section. The value of s^ is assumed 
approximately' e(iual to 0.9 s^Cobserved later) for probabilistic 




(b) MtDSPAN SECTION ^ 

FIG.3.!1 SIN'PW SUPPORTED BEAM-EXAMPtE 3.1’ 

ACt. DJMENStONS ARE W mm 
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variations of ® and o . The equation for M of a section is 
cu 8 ew 

known. Hence critical section can be fixed approximately. The 
probabilily of cracking of the critical section for each loading 
ease can be computed by Bq,. 4.13. If p^-^ is "the probability of 
cracking of the beam, for loading condition i, the bounds on the 
value of p^g of the beam for n independent load conditions is 
given by (26) 

n 

^ Pci 

1^1 

where p is the pmbabilily of cra<*±ng of the beam for n 
load conditions. The method is illustrated with examples. 

Earample 4.2 

Continuous PSC beams, having two equal spans and supporting 
a floor, are spaced at 5 m. live load on the floor is 400 kg/m^. 
Total dead load (including self-wei^t of the beam) on ary inter- 
mediate beam is 1500 kg/m. The beams are designed with concrete 

2 2 
cube strength of 350 kg/cm and steel strength of 15000 kg/cm . 

2 

The effective span and area of steel are 1 Cm and 10.89 cm . The 
beams have unifarm cross section whidi is shown in iig. 3. 12» The 
same figure is again given here for reacty reference. The beams 
are designed assuming that Ihqy are cast separately and the composite 
action of slab and beam is not considered. 

One of the intermediate beams only is considered and the 
reliability analysis of the same at limit state of cracking is 


illustrated. 
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Ihere are three critical live load conditions on the beam; 

(i) Live load on both spans 

(ii) live load on span AC 
(ixi) Live load on span CE 


Using the method explained in article 4.2.5, the critical sections 
of cradcing are C, Bg and Dg (lig. 3.12) for loading conditions (i), 
(ii) and (iii) respectively. The exterml moments, mean values and 
standard deviations of o and a remain same as m Example 3.2. 
Usirg the paraaeters and distributions of and cr^, 5000 

samples are generated for M of each critical section by Monte 
Carlo method. The i^esults are given below. 

Section Bg or Bg 

M = 46.79 tmj s = 1.31 tm 
raa rc 

Section G 

M = 51.60 tm; b_ = 1,42 tm 

rom rc 

load condition (i) ; live load on both spans 

Critical section is C. Irom table 3.2, the value of M__ 
at this section is 43.75 tm. Using Sj. 4.13, the probability 
of cracfcing of the beam fbr the loading condition (i) is 


Pol “ 


= 1.6 X 10 


-.8 


Load condition (ii): Live load on Span AC 

Critical section is Bg, Erom table 3*2j tide value of 
at this section is 28.125 tnu Using Bj. 4«13» the probabilily of 



145 


cracking of the beam for the loading condition (ii), found 

to be zero# 

SintLlarly, when live load is on the span CE, p^^ can be computed 
and fovuad to be equal to zero. 

Maximum value of the probability of cracking is for the loading 
condition (i). Hence, the bounds on the value of p , using Bq. 4.14, 
is 

Pci i Pee i Pol Po2 Po3 

As the values of p^g and p^^ are equal to zero, the value of p^^ 
is equal to Hence 

p = 1.6 ^ 10 

•^es 

Example 4«3 

Continuous PSC beams, having three equal spans and supporting 

2 

a floor, are spaced at 4m. live load on the floor is 400 kg/m . 

Totel dead load (including self-weight of the beam) on ansr inter- 
mediate beam is 1000 kg/m. The beams are designed with concrete 

2 2 
strength of 350 kg/cm and steel strength of 15000 kg/cm . The 

2 

effective span and area of steel are Em and 6,94 can respectively. 

The beams have vuaiform cross section which is shown in ilg, 3»13. 

The same figure is again given here for ready reference. !Ilie beams 
are designed assuming that they are cast separately and tiie coityosite 
action of slab and beam is not considered. 
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One of the intemediate beams onny is consideired and the 
reliabilily analysis of the same at limit of cracking is illustrated. 


There are seven different loading conditions of the beam for 

live load* (Eh^ are listed in table 4.1. Critical section for 

each loading case is also shown in the same table. The external 

beMing moments, mean values and standard deviations of and 

remain same as in Example 3.3. Using the parameters of 

and ^ and their distributions, 5000 samples are generated for M 
s xc 

of each critical section by Monte Carlo method. The results are 
given in table 4.1. 

loading condition (i) ; live load on all spans 


Critical section is 0. Erom table 3*5, the value of at C 

is equal to 16,64 tm. Prom table 4,1, the -value of of -the 

» rcm 

section 0 is 25.39 -tan. Using Bj, 4.13> the value of p . is 

C I 


Pci 


= H- 


16,64 - 23.39 
0.709 


) = 0 


Similarly the values of p^^, shown in table 4,1, can be 
calculated for other loading conditions. Using Eq, 4.14 and "the 
■values of p . from table 4.1, the botinds on p „ are obtaimd as 

-23 -23 -25 

2.08 xi0'^<p <0 + 0 + 0 + 0 + 2.08 x 10 0 + 2,08 x 10 

— -^cs — 

2.08 X < p < 4.16 X 10~^^ 
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Table 4,1, Parameters of M and probability of cracking of the 

beam (shown in Pig. 3.15) for each loading case for 

FV of o and 0 
cu s 


31. 

No. 

Live Load 
on Spans 

Critical 

Section 

M 

ew 

(tm) 

M 

rcm 

(tm) 

^rc 

(tm) 

^ci 

1 

AC, CE & m 

C or E 

16.64 

25.59 

0.709 

0 

2 

AC & EG 

or Pg 

15.45 

22.29 

0.646 

0 

5 

CE 

D 

9.27 

15.68 

0.508 

0 

4 

AC 

h 

14.67 

22 .29 

0.646 

0 

5 

AC & CE 

C 

18.38 

25.59 

0.709 

2.06(-25) 

6 

EG 

^2 

14.67 

22.29 

0.646 

0 

7 

CE & EG 

E 

18.58 

25.59 

0.709 

2.08(-25) 


Bounds on p of the beam : 2.06(-23) 1. p 5. 4»l6(-25) 

CS "“** c S 

* 2.08(-25) is read as 2,08 
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4.2.6 Probability Distribution of for Probab i l i stic 

Vsulations of o , a and Dimensions of Section 
cu* s 

Probabilistic variations of a „ and a only have been 

cu s 

considered in the article 4.2,3 to determine the probability 

distribution and parameters of It has been observed in 

CSaapter 2 that b^, b', t^, d, h and D^ have also random variations. 

Hence to be more realistic in analysis, the random variations of 

above parameters are also considered in the equations to determine 

M . 31ie section shown in Pig. 3.11 is considared, She mean 
rc 

values and standard deviations of o^, b’, t^, d, h and 

Dg of the section are taken from Example 3.4 and of b^^ and t^ 
are assijmed same as b^ and t^. Using those values and distributions 
of above random variables, 10000 samples are generated for M of 
the section by Monte Carlo method* Hie histogram and cumulative 
distribution of the generated data are shown in Pig. 4.3. Hormal 
distribution does not satisfy the chi square test for one percent 
level of significance. However it satisfies at 0,1 percent level. 

Ihe coefficient of skewness is very small (less than 0,0?) and the 
coefficient of kurtosis is 2.98 idiich is very close to the value of 
3 for true normal. Hence normal distribution is adopted for M 
for probabilistic variations of and dimensions of section. 

Having found out the probability distribution of the 

reliabilily analysis of beams at limit state of crac k ing is done 


as follows 



tm 
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rrrn 


rt 
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(i) Determine the critical section of the beam. 

(ii) Using normal distributions and parameters of b^, b^,t^, 

t, , b', d, h and D , generate sanqjles for M of the critical 
section by Monte Carlo method. Sample mean and standa3?d 
deviation of the genrated data -mill be the parameters of 

(iii) Calculate M at the critical section, 

' ' ew 

(iv) Compute p using Sj. 4,13, 

c 

Ihe procedure is illustrated with examples 
Example 4-4 

Qhe simply supported ISC beam of Example 4,1 (Big, 3,1 1) is 
considered. Dimensions of the section and strengths of materials 
are subjected to random variations, 

Qhe critical section is liie midspan of "Hie beam, Bae value 

of M at the critical section is 43,75 tnu Using Monte Carlo 
ew 

techniquei the values of M^^ and s^ have already been found in 
the previous article for this beam. Erom Big. 4,3f 

M_„ = 67.42 tm } s^ = 3,99 tm 
rom re 

-9 

Using Bli, 4,13, 1iie value of p^ is found to be 1,5 10 . 

ahe reliability analysis of continuous beams at limt state 
of cracking and for probabilistic variations of o^ and 

dimetisions of section remain same as explained in the article 4*2,5 
except that probabilistic variations of ^ 
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and D_ are also to be considered in deteimining the parameters of 
Ihe results of the reliabilily analysis of the two span 
beam of iig* 3.12 are given in table 4.2 and for lie three span 
beam of Pig. 3.13 in table 4.3. 

4.2,7 Oomputation of p^ for Probabilistic load 

It is known from Chapter 2 Ihat live load is distributed as 

kg/m ,0,368), The dead load is assumed as deterministic 

as in the previous chapter. The resistance of the section at 

initial crack formation is normally distributed as determined in 

articles 4.2,4 and 4.2.6, When the load and resistance are 

lognoimally and normally distributed, the probabilily of cracking 

is calculated by using Bi. 3.51, Using the above equation, ttie 

value of p for probabilistic load is given by 
c 


log(- 




s + M - M 
rc rcan g ^ 


Po“ 


= / 


e 




^ [1 

/ir 


M. 


qm 


• }’ ] dw 
(4.15) 


The reliabilily analysis of PS0 beams at limit state of cradking, 
considering probabilistic variation of Q, is illustrated witii 
examples. In all the following examples in this section, live 
load is distributed as Iiff(i21.77 kg/m , 0.368), 


Example 4.5 

The simply supported PSC beam of Example 4.1 (Pig. 3.11 ) is 
subjected to probabilistic variations of live load and strengths 
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Table 4 * 2 . Parameters of and probability of cracking of the beam 
in Pig. 5*12 for each loading case for PV of o ,0 and 
dimensions of section ^ ® 


SI. 

llo. 

Live Load 
on Spans 

Critical 

Section 

ew rom 

( tm) ( tm) 

s 

rc 

(tm) 

^*ci 

1 

AC and CE 

C 

43.75 58.51 

5.36 

7 . 5 (-^)* 

2 

AC 


28.25 52.91 

3.05 

i.98( -16) 

3 

CE 

^2 

20.25 52.91 

3.05 

1.98(-i6) 

Bounds 

on Vos 0^ 

the beam ; 7*5(-6) < p„_ < 7»5(“6) 

C S 


* 7*5('^) is read 

-6 

as 7 . 5 ><lO . 




Table 4 » 5 » Parameters of and probability of cracking 

of the 


Beam 

ou * 

in Pigf 3*15 for each loading ^ 
0g and dimensions of section 

case for TV 

of 

SI. 

Ho. 

Live Load 
on Spaas 

Critical 

Section 

ew rom 

( tm) ( tm) 

s 

rc 

(tm) 

I>ci 

1 

ACj CE & EG C 03?^ E 

16.64 28.73 

I.6O7 

2 . 5 (- 14 ) 

2 

AO & EG 

or Pg 

15.43 25.24 

1.457 

9.0(-12) 

5 

CE 

D 

9.27 17.85 

1.046 

1.1(-16) 

4 

AO 

C 

14.67 25.24 

1.457 

1 . 8 (- 15 ) 

5 

AO & CE 

% 

18.38 28.75 

1.607 

6 . 0 (-il) 

6 

EG 

E 

14.67 25.24 

1.457 

1 . 8 (- 15 ) 

7 

CE & EG 

^2 

18.38 28.73 

1.607 

6 . 0 (-ll) 

Bounds 

on Vos 

the beam ; 6 . 0 ( 

- 11 ) iPcs 112 . 95 (- 11 ). 
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of mateiials. Other particulars of the beam remain same as in 
Example 4«1 

Erom Example 4.1, the parameters of M of the midspan of the 
beam are 

= 59.56 tm ; s = 1,64 tm 

rcm rc 

Erom Example 3.7, 

M = 18,75 tm ; M ^ = 7.61 tm, 
g qm 

Substituting the values of M , s , M„ and in Bj, 4.15 and 

^ rcm' rc* g qm 

solving the same by using Gaussian Quadrature formula (explained in 
14ie previous chapter), the value of p^ is found to be 2.94 10”^, 

The reliabilily analysis of continuous beams at limit state 
of cracking for probabilistic variations of Q, and remain 

same as explained in the article 4.2.5 except that p is evaluated 
using Bq, 4.15. 

!Qie same tvfo span and three span continuous FSC beams of 
Examples 4*2 and 4,3 are considered. The values of M and M at 

O MM*' 

the cxltloal sections of the beam of Example 4.2 are available in 

table 3.17 and p£irameters of M in Example 4.2. SHie values of 

rc 

M and M at the critical sections of the beam of Example 4.3 are 
g ga 

available in table 3 *20 and parameters of in Example 4.3. 
results of the reliability analysis of the above beams at limit 
state of orackii^ are presented tn tables 4.4 and 4,5. 
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Exan^le 4*6 

Qhe siii 5 )]y supported PSC beam of Example 4» 1 is subjected to 

probabilistic variations of live load, strengths of materials and 

dimensions of the section. The live load is distributed as 
2 

Iill(l21.77 kg/m ,0.368). Other particulars of the beam remain sam^ 
as in Example 4.1. 

Using Monte Carlo method, the parameters of of the midspcr. 
of ttie beam for probabilistic variations of and dimension'' 

of the section have alrea(^ been obtained in Example 4.4. Ihe 
results axe 

M = 67.42 tm; e = 3.99 ta 
3X5IE rc 

Erom Exanqjle 3.7, 

M s= 18,75 tms M = 7.61 tm 
g qm 

Tfeing Ej* 4,15, the value of p of the beam is con^juted and found 

c 

-7 

to be 4.87 X 10 . 

Using 1iie parameters of of the critical sections of tho 
beams of Example 4.2 and 4.3 and tdis valnes of M and M at those 

o 

critical sections from tableB4.4 and 4.5, the vahies of of tho 
above beams for different loading conditions can be calculated. 
Hesults are presented in tables 4.6 and 4. 7* for two span and three 
Span continuous ISO beams (Examples 4.1 and 4.2) subjected to 
piobabilistio TaiiatiDna of Q, o^, WA dUBnoiODS of seot'jor . 
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Table 4*6 « Parameters of aad probability of cracking of the 

beam in Pig. 3.12 for each loading case for PV of Q, 

o , a and dimensions of section 

on s 


SI. 

No. 

Live Load 
on Spans 

Critical 

Section 

M 

S 

(tm) 

M 

q^m 

(tm) 

M 

rcm 

(tm) 

®xc 

(tm) 

^oi 

1 

AO & CE 

C 

18.75 

7.62 

58.31 

3.36 

7.02(-6)* 

2 

AC 


9.58 

5.67 

52.91 

5.03 

4.49(-e) 

3 

CE 

®2 

9.58 

5.67 

52.91 

5.03 

4.49(-8) 

Bounds on p^g of 

the beam 

: 7.02(-6) 

< P 

— -^CS 

± 8.0l(-6) 



* 7.02(-6) IS read as 7.02 

xio"^. 





Table 4.7» Parameters of M and probability of cracking of the beam 

in Pig. 3.15 eaoh loading case for P7 of Q, 0 , a and 

dimensions of section ° ® 

SI. 

No. 

Live Load 
on Spans 

Critical 

Section 

M 

g 

(tm) 

M 

qm 

(tm) 

M 

rcm 

(tm) 

s 

TO 

(tm) 

^^Ci 

1 

AC,GE & EG 

C or E 

6.4 

3.12 

28.75 

1.607 

8.32(-8) 

2 

AC & EG 


4.8 

3.15 

25.24 

1.457 

3.24(-7) 

3 

CE 

D 

1.6 

2.34 

17.85 

1.04£ 

1.60(-7) 

4 

AC 

% 

4.8 

2.85 

25.24 

1.457 

7.92(-8) 

5 

AC & CE 

C 

6«4 

3.65 

28.75 

I.607 

7.3l(-7) 

6 

EG 

P 

^2 

4.8 

2.85 

25.24 

1.457 

7.92(-8) 

7 

CE & EG 

E 

6.4 

5.65 

28.75 

I.607 

7.3l(--7) 


Bounds on of the heam : 7*5l(-7) 1. P„o — 21.88(“7)« 

09 OS 
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It is observed that the probability of cracking of simply 

-9 -22 

supported beam varies from 10 to 10 for deterministic load 
-6 -7 

and from 10 to 10 for probabilistic load. Ihe probabilily of 

-6 

csracking of continuous beams varies from 10 to veiy small value 

-24 -5 -7 

less than 10 for deteiministic load and from 10 to 10 for 

probabilistic load. The probability of cracking of beams are 

hi^er for probabilistic load than deterministic load. 

4.3 EELIABILITK' MAiySIS OP PSC BEAMS AT LIMIT STATE OP BEPLEOTION' 

4.3.1 Introduction 

Under working loads PSC beams do not crack usually. Since 
prestressed concrete is assumed as a homogeneous elastic body which 
ob^s quite closely the ordinary laws of flexure, the deflections 
can be computed by methods available in elementary strength of 
materials (58). Deflections due to prestress can be computed by 
considering the concrete as a free body separated from the tendons 
which are replaced by a system of forces. Total deflection is 
calculated by superposing deflections due to load and prestress. 
Codes specify that deflections must not exceed a certaiin fraction 
of span when 

(i) prestress and dead load are acting and 

(ii) prestress, dead load and live load are acting. 

Hence the beam is said to become unserviceEdile and to have reached 
serviceabilily limit state caused due to excesslTe deflection if 
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the deflection, u, caused by external load exceeds the allowable 
deflection specified by the code. If is defined as the 
probabilily of the stiucture becoming unserviceable due to excessive 
deflection constraint, then 

p^ = t( u > allowable deflection) (4. 16 a) 

Allowable deflection is specified as the fraction of span, s^ 

Tfdiere 'I' represents the fraction, Kbw B 3 .. 4.16a can be rewritten 
as 

Pd = ^( I ^ 'J' ) (4-l6b) 

Denoting 4.l6b becomes 

p^ = p(v > Tp ) (4.16c) 

Qlie deflection of a prestressed concrete beam is a function of 
several random variables namely, properties of materials, load and 
geometrio properties of the section. Hence the deflection is a 
random variable. Since deflection is a random variable, v is also 
random variable viaose distribution and parameters are best obtained 
by Monte Carlo method. 

4.5.2 Bg,uations for Determination of Deflection in PSC Beams 

It is assmed that cables are fully borded and the cable of 
the beam passes throu^ the centroid of tiie sections at ends of the 
beanu 

Ihe downward deflection, u^, of a PSC beam under dead load, live 
load and prestresa can be written as 
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u, = u + u - u 
d g q p 


(4.17) 


where u . u and u are deflections due to dead load, live load 

rip 

and prestress respeotive3y. Por parabolic cable profile, the 
equiTalent uoifOrmly diatrobutedload, Wp, due to prestress is 


given "ty 


es 


w = 
P 


(4. 18) 


where P is the total prestressing’ force in steel, y^ is the 
central ordinate of the parabolic cable profile and ^ is the 
effective span. Bbr unifomOy distributed dead load and live 
load, Bi, 4.17 can be written as 
,4 

— ( C w + 

11 P p 


u. = ^ (c„ w„ + c^w„ - c w ) 




(4.19) 


where I is the moment of inertia of the transformed area of cross 

section and c , c and c are constants which are ei’ther known 
g' 1 P 

or can be calculated for given loading and the position where the 
value of deflection is required. Taking idle value of effective 
tensile stress in prestress steel equal to 0,6 and using Bq, 4.11, 


the f i nal equation for v^ is obtained as 


/ dN Jf- ^ 

r M-~^) = " . — Cc„ + c„w - 

^ ^ 18000 j S g 1 q 

cu 


%^s Wc 


) 

(4,20) 


where v^ is the ratio of downward deflection to ’the effective 


span of the beam.. 
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The upward deflection, u^, of a PSO beam, tinder prestress and 


dead load is 


u = u - u 
u p g 


(4.21 ) 


Under dead load conditions, the bottom fibre of the beam is under 
compressive stress which majr be present for considerable period 
when Ihe live load is not acting. Under this substained stress, 
the beam has a tendency to deflect upwards due to creep, Ohe I.S, 
Code (46) gives the jEbllowing expression for this deflection 


L ds 
u 

■u- = — ST” 

cc 8h 


(4.22) 


■vdiere, 

u = deflection due to creep 
00 

= unsupported length of beam 
As = differential shortering between top and bottom 
surfaces of the beam due to creep 

Assuming creep coefficient of concrete is equal to 2,5 and stresses 
in concrete at extreme fibres of the beam are equal to allowable 
values 0.5 0 and zero, the Bq* 4,22 reduces to 

Cvl 


1 a 

- u cu 
%c “ 64h E 


(4.23) 


Therefore final equation for similar to Bj, 4.20, can be 


written as 


\ - (it) = 




18000 /a I 
cu 


“ 0 W ) 

s 


L 

^ cu 
11:5200(]h 


(4.24a) 



where is tiie latio of total upward deflection to the effective 
span of the beam. Since the moment of inertia is a function of 
geometric properties of the section. Bis, 4 . 20 and 4.24a are functions 
of strengths of materials, load and geometric properties of Uie section. 


If creep of concrete is not considered, the prestress in steel 
at transfer shovild be used instead of effective prestress in deter- 
mining the deflection due to prestress. Using the maxinum value of 
sti^ss in steel at transfer equal to 0,8 0 ^( 46 ), v^ is expressed as 


6,4 c o„ 




Wc 


18000 


/o I 
cu 


- c w ) 
g g^ 


(4.24b) 


4.3,3 Probability IDistributions of v^ ani for Probabilistic 

Variations of o and o . 

cu s 

Hie probabilistic variations of 0 ^^ and are first considered 

in the equations for determination of v^ and v^ and the reliability 

analysis at limit state of deflection is fiirst presented for 1die 

same, o and o are distributed as ir( 422. 8,56 ) kg/cau 
CU- s * 

N( 15680,488) kg/cm respectively. Hie beam of Example 3,1 (lig. 3.1l) 

is considered, Ihe critical section is midspan of the beam vdiere 

maximum deflection occurs. The values of c , c and c are same 

g’ q p 

and equal to 0,01302, Using Bqs, 4,20 and distributions and 
parameters of and Ogj 10000 samples are generated for v^ by 
Monte Carlo method, She histogram and cumulative distribution 
of the generated data are shown in Hg, 4.4* lognormal distribution 
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IS found to satisfy the generated data at one percent level of 

significance. Similarly samples have been generated for v^ using 

Bq., 4,24a and the same lognoimal distribution has been observed for 

this data also. Hence lognormal distribution is adopted for random 

variables v, and v for probabilistic variations of o and o . 

d u ^ cu s 

4.3-4 Computation of p^ 

GIhe probability of toe structure becoming xinserviceable 
due to excessive deflection constraint is given by Bq, 4.1 6c, 

I.S. Code (46) specifies values of 'I' for v, and v equal to 0.002 
and 0.0033 respectively. Therefore Bq. 4.16c for limit state of 
downward deflection becomes 


= P(Td >0.002) (4.25) 

and for limit state of upward deflection 

= P (v^ > 0.0033) (4.26) 

where, 

P^^ = the probability of the structure becoming unserviceable 
due to excessive downward deflection constraint 
p^^ = the probability of the structure becoming unserviceable 
due to excessive upward deflection constraint. 


Por lognormally distributed v^, the value of p^^ can be computed 
from the following equation 


log (• 


s 


0.002 

vd 


) 


Pad = t' 


] 


(4.27) 
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whe3?e and s_^ are the parameters of Similarly for lognormally 

distubuted the value of p, can be calculated from the equation 


log (^) 

Pau"i - 


s 


(4.28) 


vu 


where v^^ and are the parameters of v^, Bie procedure is illustrated 


wx1±i an example. 


Example 4-7 

The simply supported PSC beam of Example 4.1 (Eig. 3-11 ) is 
subjected to probabilistic variations of strergths of materials. Other 
particulars of the beam remain same as in Example 4, 1. 

a and a are distributed as H(422.8, 56) kg/cm and 

wU O 

N( 15680, 488) kg/cm i^spectively. Using Monte Carlo technique lOOOO 

samples are generated for v^ and v^ at midspan of the beam, given by 
Bjs, 4f20, 4.24a and 4.24b. The results are given below. 

= 56.56 X lo"^ ; s^ = 5.42 x 10“^ 

\aa “ ^ ’ ®um 4.24a) 

\mi "" ^ » ®um ^ 4.24b) 

using similar equations 2.28 and 2.29, the parameters of v^ and v^ 
can be computed. IheLr values are given below. 

■^dm xIO”^ ; s^ = O.C^56 

’’um =^5.11 X ^0^ . s^ = 0.11 (Bor Sb. 4.24a) 

V = 85.88 xio“^ : s = 0.091 

um ' ini 


(Btor El. 4.24b) 
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using Bg., 4«27, the value of is calculated as 


0.002 




= 1 - <|> 


log (■ 

^ 56.3 X 


r) 


10 


0.0956 


] 0 (i.e. less than 10 


If deflection due to creep of concrete is considered, the value of 
^du’ usxng Bg, 4.28, is 


log (- 


0.0033' 


•du 


= 1 - <j) [ 


45.11 X 10 


-5‘ 


0.11 


1 ilO 


If deflection due to creep of concrete is not considered, then 


log (- 


0.0033 


Pdu = ■> - <!> t 


85.88 X 10 


0.091 


•3 1 0 


Hence the probability of the beam becoming unserviceable due to 
excessive deflection constraints is zero. 


4.3.5 Reliabilily Analysis of Oon'tinuoxis PSC Beams at limit State of 
Deflection 


Ihe reliabilily analj^sis of continuous PSC beams at limit 
state of deflection is carried out in the following steps. 

(i) Determine the position of maximum deflection for each loading 
condition. 

(ii) using Monte Carlo technique, obtain parameters of v^ and v^ 
of each critical section. 

(iii) Compute p^^ and p^ for each loading condition using Bgs. 4.27 


and 4.28. 
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(iv) Using El, 4,14, calculate bounds on p which defines the 

QS 

probability of the beam becoming unserviceable for all loading 
conditions due to excessive deflection constraint. 

Ihe procedure is illustrated with examples. 

Example 4,8 

(Hie two span ISC beam of Example 4.2 (lig. 3,12) is 
considered. The beam is subjected to probabilistic variations of 
and 0g, Other particulars of the beam remain same as in Example 4.2. 

o and CT are distributed as N(422.8, 56) 'k.g/cm and 
2 

N( 15680, 488) kg/cm respectively. Critical sections where maximum 

deflection occur under different load conditions are given in table 

4,8. The values of constants c^ , c and c are also given in the 

g q. P 

same table. Using Monte Carlo method 10000 samples are generated for 

Vj and V of each critical section and the mean value and standard 
d u 

deviation of Vj and v are obtained, Dhe results for the critical 
d u 

section X^(o,422 2. from the left support) for load condition (i) are 
given below. 

Load condition (i) ; live load on both sptins 

Critical section is X^, 3he mean value and standard deviation 

of Vj and v are 
d u 

= 24,81 X 10~^ ; s^ = 2,38 x 10”^ 

V = 20,60 X 10“^ ; s = 2.22 x 10”^ (Ear Bl. 4,24a) 
jum ' vu 

= 38.46 X lo"^ } s^ = 3.49 X 10“^ 


(Eor El. 4.24b) 



TalDle 4»S» Ciitioal sections and values of c , c and c for 
two span continuous 'beam SI P 


SI. 

No. 

Live Load Critical 

on Spans Section 

Distance 

of Critical °g °q ®p 

Section from 

End A 

1 

AO (Sc OE 


0.422 £ 

0.00542 0.00542 0.00542 

2 

AC 


0.485 A 

0.00527 0.00915 0,00527 

5 

CE 


1.517 A 

0.00527 0,00915 0.00527 

Table 4 • 9 * Values of axid of the 

of a aud a 
ou 8 

beam in Pig. 5.12 for PV 


Live 

Limit State of Lownwaid < 

Critical Limit State of 

SI. 

Load 

Deflection 

Section Upward Deflection 

No. 

on Spans 


%d 

-5 

Y X 10 s p. 

um YU du 

1 

AC & OE 

24.7 0.0957 

0 

X 20.48 0.106 0 

* 58.50 0.091 0 

2 

AC 

66.95 0.0722 

0 


3 

CE 

66.95 0.0722 

0 

"3 


The value of is equal to zero. 

* For Eq, 4 » 24 'b (without considering creep of concrete). 
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— 2 38 

Coefficient of variation of v,„ = 6 , = = 0.096 

am vd 24.81 

Using similar Bis. 2.21 and 2,22, 

s , = *^log (0.096^ + l) = 0.0957 
va 

' ^dm " 2 

= 24.70 X 10^ 

Similarly parameters of v^ are calculated and given in table 4.9. 
Using Bis, 4.27 and 4.28, the values of and p^^ are calculated 
as illustrated in the previous example and found to be zero. 

Similarly the values of parameters -of v^ and v^ for other 
load conditions are obtained and shown in table 4.9. (Ehe values of 
p^^ for other load conditions are calculated in a similar manner and 
given in the same table. In every loading case, the value of p^^ is 
found to be zero and hence the probability of the beam becoming 
unserviceable due to excessive deflection constraint for all load 
conditions is zero. 

Example 4,9 

ilhe three span continuous ISC beam of Example 4.3 (K.g, 5.13) 

is considered. Ihe beam is subjected to probabilistic variations of 
0 and a , Other particulars of the bean remain same as in 
Example 4*5. 

Distributions and parameters of and remain same as in 
Example 4,3* Ilhe critical sections for lindt stetes of deflection 
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for each loading case is given in table 4.10, The values of c , c 

S ^ 

and c are also given in the same table. Using Monte Carlo method 
P 

10000 samples are generated for v^ and v^ of each critical section. 

Parameters of v, and v are calculated from the mean values and 
d u 

standard deviations of v, and v obtained flrom the Monte Carlo method. 

d u 

Parameters of v, and v for all critical sections are given in table 4,11 
d u. 

Using Bqs, 4.27 and 4.28, the values of and are computed as 
illustrated in the previous examples. Table 4.11 shows the values of 
p^^ and p^ for each loading case. It is found that values of p^ is 
zero in eveiy case and so the probability of the beam becoming unservice- 
able due to excessive deflection constraint for all loading conditions 
is zero. 

4*3.6 Probability Distribution of v^ for Probabilistic Variations 
of , Og and Q 

!Qae probabilistic variations of Q considered 

in the equation 4,20. Eor reliability analysis, the probability 

distribution of v^ is required. To study the probability distribution 

of v^, the beam of Ikample 4.1 is considered, , 0^ and Q are 

distributed as ir( 422 . 8 , 56) kg/cm^, N(15680, 488) kg/om^ and 

1^(608,85 kg/m, 0,568) respectively, Ihe maxinum deflection occurs 

at midspan of the beam. The values of c , and o are same and 

g <1 P 

equal to 0,01302. Using Bj. 4.20 and distributions and parameters 
of Q 10000 samples are generated for v^ by Monte Carlo 

method. 33ie histogram and cumulative distribution of the generated 
data are shown in Big. 4.5. distribution is hi^ly positively 
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TaMe 4*1 0» 


Critical section aTid values of c , 

6 

span continuous 'beam 


c and o for three 
q P 


SI. 

No. 

Live load 
on Spans 

Critical 

Section 

Distance 
of Critical 
Section 
from End A 

c 

e 


0 

P 

1 

AC,CE & EG 

^1 

0.444 

£ 


0,00693 

0.00693 

0.00695 

2 

AC & EG 


0.496 

i 


0.00677 

0.00990 

0.00677 

5 

CE 


1,5 a 



0.00052 

0.00677 

0.00052 

4 

AC 


0.416 

1 


0.00684 

0.00887 

0.00684 

5 

AC & CE 

"5 

0,483 

Z 


0.00685 

0,00587 

0.00685 

6 

EG 


2.584 

Z 


0,00684 

0.00887 

0.00684 

7 

CE & EG 


2.517 

z 


0.00685 

0.00587 

0.00685 

Table 4»1'1» "Values oi and i 

■du 

of the "beam in Eig. 

3.13 for 


PV 

Of and 






Live Limit State of PomOTaid 


Limit 

: State of Up- 

SI. 

No. 

Load on 

Deflection 



n„_. X-* - -1 ward Deflection 

Spans .^-5 

’’dm 

®vd 


% 

Section v x 
tun 

10 ^ s p, 

vu '^du 

1 

AC,CE & EG 

31 .62 

0,0960 


0 

^1 

29.85 0.0996 0 








* 51 .89 0.0872 0 

2 

AC&EG 

58.05 

0.0763 


0 




3 

CE 

65.55 

0«o686 


0 




4 

AC 

48.90 

0.081 


0 

^4 



5 

AC & CE 

22.29 

0.116 


0 




6 

EG 

48.90 

0.081 


0 




7 

C® & EKJ 

22.29 

0.116 


0 





The value of is equal to zero, 

*Por Bq., 4»24h (without considering creep of concrete). 
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skewed. ISie random vaiiable has taken positive and negative values. 
This IS true in realily also as prestressed concrete beam can deflect 
upwards and downwards depending upon the magnitude of live load on the 
beam. The standard distribution which has these characteristics i*e. 
hi^ positive skewness and positive and negative values for random 
variable, is QJype I (largest) extremal distribution. This distribution 
does not satis:^ chi square test because of large values of chi square 
on the lower tail even though ihe distribution fits the data well on 
the remaining portion. However as this is the most conceivable distribu- 
tion for the deflection behaviour of prestressed concrete beams, Hype I 
(largest) extremal distribution is assumed for v^ for probabilistic 

variation of o.d and Q. 

ou ' s 

ibr the Type I (largest) extremal distribution, the cumulative 

probability of v_ is given by (40) 
d 

\ (t^) = exp [ -exp { -s^ (v^ - v^)}] (4.29) 

4 

*00 < Vj < ^ 

d — 

in vftiich v^ and s^ are parameters of v^. The above parameters are 
evaluated from the following equations (48). 


^dm 


y ..^'7 

®Vd 


(4.30) 




Td 


n . 


(4.31) 


The format for tire probability of Idae structure becoming 
unserviceable due to excessive downward deflection constraint is given 
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by Bj. 4.25- Hence if is assumed to follow O^pe I extremal (largest) 
distributioti, the value of using Bg.. 4.29, is obtained from the 
following equation. 

Pdd = ■' " ^ "V ”■ 

The method of reliability analysis at limit state of deflection for 
probabilistic variation of and Q is illustrated with an example. 

Example 4.10 

Bie simply supported beam of Example 4.1 (lig. 3.1i)is analysed 

at limit state of deflection and for probabilistic variations of o . a 

cu’ s 

and Q. 


Using Monte Carlo method, tie values of v^ and s^ have already 

been obtained for this beam in the previous section, ilrom Eig, 4.5, 

the values of and s^ are -11.99 x and 13.22x10“^ respectively. 

Using Bis^ 4Tr30 and 4*31^ the parameters ot v- are calculated as 

d 


TT X 1 0^ 5 

s__. ^ * — ^ 0.0965 X 10 


^ 13.22 


"a. = 


= - 17.97 10 


Substituting the above values of and v^ in Eg. 4.32, the value of 
p^^ is obtained as 7,27 x lo"^^. 


calculated in Example 4.9. 


Ihe value of p, remains same as 
du 
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!Hae reliabilily analysis of contimous beams at limit state of 
deflection and for probabilistic variations of 0 ^, and Q is same as 
explained in article 4.3.5 except that probabilistic variation of Q is 
also to be taken into account in deteimining parameters of v^ and Ej. 4-32 
must be used t6 compute Ihe results of the reliabilily analysis 

at limit state of deflection and for probabilistic variation of 0 ^^, 
and Q for continuous ISC beams of Examples 4.2 and 4.3 are given in 
tables 4.12 and 4.13 respectively. 

4 . 3.7 Probability Distributions of v^ and v^ for Probabilistic Variations 

of 0 f CT # Q and Dimensions of Section 
cu ' s’ 

In the previous sections probabilistic variations of strengths 
of materials and live load only have been considered in determining the 
probability distributions of v^ and v^. As it is known from chapter 2 
that dimensions of section of a beam have also random variations which 
are now taken into account in liie equations of v^ and v^. Tb.e beam of 
Example 4.1 (Pig, 3.1 1 ) is considered, Bie mean values, standard devia- 
tions and distributions of o , 0 , Q,b., b', t., t, , d, b, , h and d 

are available in the Example 4.3. Using them and prediction Bjs. 4.20 
an d 4.24a, 10000 samples are generated for v^ and v^ by Monte Carlo method. 
Ihe histograms and cumulative distributions of v^ and v^ are shown in 
Pigs. 4.6 and 4.7 respectively. It is observed frcmi the Pig. 4.6, that 
the distribution of the generated data ffor v^ is hi^ly positively skewed 
and the random variable v^ takes positive and ne^tive values. Even 
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For Eq.. 2.24a (without considering creep of 

-15 

1.0(-15) i® as 1#0 ^ 10 • 
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though !iype I (largest) eactjremal distribution does not satisfy chi 
square test at one percent level of significance, the same distribution 
is assumed for v^ based on reasons given in the previous article- 
Estimation of parameters of the above distribution and computation of 
p^^ remain same as explained m the previous article. 

logno rm al distribution is found to satis:^ the chi square test 
at one percent level of significance for the generated data shown in 
Eig. 4.7 for Uie random variable v^. Hence lognonnal distribution is 
adopted for v^ for probabilistic variation of strengths of materials, 
load and geometric properties of the section. 

An example is presented to illustrate Ihe method of reliabilif 

analysis at limit states of deflection and probabilistic variations of 

o t ^ t Q dimensions of section, 
cu^ s 

Example 4.11 

!I5ie simply supported beam of Example 4.1 (Pig, 3.1 1) is 
anaiysed at linat state of deflection with probabilistic variations 
of 0 f , Q and dimensions of section. (Qie live load is distributed 
as I2J(608,85 kg/m, 0.368), Other particulars of the bean remain same 
as in Example 4,1. 

Using Monte Carlo method, the values of v, , s_, , v and s 
^ » dm^vd’umvu 

at the midspan of the beam have already been determined and sihowh in 
Pigs. 4.6 and 4.7. Ercmi the above mentioned figires 
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^dm “ X 10“^ ; = 10.0 x 10*^ 

Tyjjj = 35.3 xio”^ ; = 4-. 01 X io”*^ (Por Bi. 4.24a) 

= 66,0 xIO ^ = 6.35 ^ 10*^ (itor Hj. 4.24b) 

Assuming 5ype I (largest) extremal distribution for the parameters 
of are cailculated by using Bg.s. 4.30 and 4.31. 

= -14.91 X io“^ 5 s , =0.127 X 10^ 
dm vd 

-12 

Using Bi, 4.32, the value of p^^ is obtained as 1,8l x 10 , 

Bor lognormally distributed v^, the parameters, calculated 
using similar His, 2.21 and 2,22, are found to be 

T = 35.07 xio"^ • s =0.113 (ihr Bi. 4.24a) 

■UR VU 

V = 65.69 10"^ } s = 0.096 (R)r Hi. 4.24b) 

um ' vu 

Using Hi, 4.28, the values of p^ for both cases, i.e. considering the 
deflection due to creep of concrete and without considering the 
deflection due to creep of concrete, are found to be zero. 

Gllhe reliability analysis of continuous beams at limit state 
of deflection vdth probabilistic 'variations of Og » Q 

dimensions of section raaain same as explained in 'Hie previous 
section except that probabilistic varla'tions of dimensions of section 
are also to be taken into account in de-terming parameters of v^ and y^. 
She results of the reliabilily analysis of the two span and three span 
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continiJDus PSC beams shown in Pigs, 5.12 and 3.13 are piesented in 
tables 4.14 and 4.15 respectively. 

It is observed that the probabiliiy of PSC beams becoming 

unserviceable due to excessive deflection constraint at transfer of 

prestress is almost zero. Hae probability of PSC beams becoming 

unsezviceable due to excessive deflection constraint at design load 

is almost zero for deterministic load. IVor probabilistic load, the 

values of p,, for continuous beams are found to be smaller than the 
■^dd 

values of p^^ for simply supported beems as expected. 



Table 4 . 14 , Values of and of the beam in Pig. 5.12 for PV of 
dim«isions of section 
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CHAITJR 5 


RELIABILIOr MADTSIS OP KESTRESSED COECEETE BEMS AT 
mAHSPm OP IRESTRESS 


5.1 lETROnJGTION 

Different stages of loading are to be considered in the 
analysis and design of PSC beams. In Chapter 3» reliabilily analysis 
of PSC beams at limit state of strength at design load were presented. 

In this chapter, reliability analysis of PSC beams at limit state of 
strengthat transfer of prestress is given considering probabilistic 
variations of : 

(i) strengths 6f materials 

(xi) stirengths of materials and geometric properties of the section. 

The failure of botthm flange under compression of concrete at 

transfer of prestress occurs if tiie strain in 1±ie bottom flange reaches 

e . Ohe limit state is reached vhen the strain in bottom fibre is 
c 

equal to or the stress at the bottom f large is equal to the compre- 
ssive strength of concrete. At this limit state, (i) total compressive 
force in concrete, C , should be equal to the prestressing force at 
transfer P.j. and (ii) the distance of the centroid of the area under 
coapression from centre of gravity of steel must be equal to (m^/P.^.). 

As in the case of failure under majdLnum loading, it is assumed that 
the stress distribution can be represented by a rectangular diagraa 
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(uniform stress) over a deplii ‘a'. 5Che strain at the extreme fibre 
in compression is equal to the rupture strain 0,0035 and the uniform 
stress IS equal to 0,68 o^, Etor calculating probability of failure 
or diecking safely, the failure can be said to occur vriaen the resisting 
capacity of the section is less than 1116 prestressing force. If is 
the area under compression (shown hatched in ilg. 5«l)» then the proba- 
bilily of failure of the section at transfer of prestress, p^, is 
given by 

Pt ■= r(0.68 < P^) (5.1 ) 

is a function of P^, the above equation can be rewritten as 
p^ = p(k < l) (5.2) 

where 

« = 0-68 Pott (5.3) 

1 C represents the ratio of resistance to action. 

5.2 BQUAPIOirS IQR DETERMIFAIIOir OP ic 

5Qie failure of the bottom flange of a flanged section at 
transfer of prestress may take place with ary one ot "the following 
events s 

“ Ihe event that a < h - t^ (Pig« 5,1a) 

- the evait that a > h - t^ (l^« 5.1b) 

PoUowing equations are ^ven to determine ic . 


Since S^ 
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Event E^ (ELg. 5.1a) 

If 'a* is less than (h - t^)» the area of concrete in coaipression 
is given by 

= b' a + t^ (b^ - b') (5.4) 

where t^ is the thickness of bottom, flapge. At this limit state, 

the distance of centroid of from 1iie line of action of must 

b t 

be equal to the shift of line of action of C from the line of 

c 

action of Pj. wdiere C denotes compressive force in concrete. Using 
"C c 

the above condition, the following equation is obtained, 

Sb =r '"(f- V 

g 

where is the distance of centroid of steel from bottom. Solving 
Bqs, 5,4 and 5,5, the value of 'a' is expressed as 

a = (b 2 +/b 2 + 4b^ b^)/(2b^) (5.6) 

where 

^Pt b' 

N = 2Mg 

- g 

6 

!I!he value of 'a* obtained by Eg,. 5.6 is aubstitutecl in Hj, 5,4 to get 
the value of Knowing the value of x can be computed by using 
Si. 5.3. 
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Event Eg (lig« 5.1b) 

If is greater than (h “ ibe equation for the value of 
'a* can be similarly derived as for event and given by 

a = (og + / Cg + 4c^ c^)/( 2 c^) (5.?) 

where 

o, = P/( 2 Mg) 

°2 ' ’■t "b 

O 3 =t^ (b^ -»') tl (x^ - 0.5 tj,)] 

g 

+ (b^ -b-) (tj -h) (1 

g 

!Dhe value of Sj^ is expressed as 

S]j = b*a + t^(bjj - b») + (a - h + t.^.) (b.^. - b*) ( 5 .&) 

Knowing Sj^, the value of k can be computed by using Eb, 5.3. 

Ihere is a possibility that failure may occur at the ends of 
a ISC beam. The failure referred here is the direct strength 
failure. Secondaiy stresses (bursting and spalling) developed in 
the end block (einchorage zone) as a consequence of transfer of pres- 
tressirg force at the ends of the beam are not considered in this 
thesis, iisstmiing the centroid of steel coincides with the centroid 
of the section at the supported end of the beam, the value of for 
rectangular shape (b.j. x h) of the end section can be eoqpressed as 
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Sjj = h (5.9a) 

In case, the centroid of the cable has an eccentricity at ends of the 
beam, the value of for rectangular section is 

*= 2b. (h - d) (5.9b) 

b T! 

Knowing S^, k can be computed by using Bi, 5.3. 

!I!he stress in prestress steel is assumed as 0.8 (j , the maziEum 

o 

value specified by I.S. Code (46). Hence 

It can be seen in all oases, k is a function of t a t M , A. _ 

» s ’ cu ' g^ ts 

and dimensions of section. Monte Carlo method has been used to 
detennine probability distribution and parameters of tc . 

5.3 HlOBABIIITr DISTRIBUTION OP k POR IROBABIHSTIC ViffilATIONS OP 
STRMGTBS OP MATERIALS 

Probabilistic variations of c and are first considered 

cu s 

in governing Bis. 5.4 to 5.10 used to detemxne tc* The prestressed 
concrete section shown in ELg. 3.11 is considered. a and a_ are 

CU s 

2 2 

distributed as N(422.8, 56) kg/cm and N(15680, 488) kg/cm respectively. 
Using Monte Carlo method, 10000 samples are generated for k of the 
section using Bis. 5.4 to 5.8, The histogram and cumirlative distributior 
of the generated data are shown in KLg. 5.2. The normal distribution 
is found to satis:^ the chi square test at one percent level of signi- 
fioance. Bene© normal distribution is adopted for k fox probabilistic 
variations of o^, and o„. 

CU s 




(g) H;g'!OGRAM 



(b) CUMULATIVE DiSTRIBUTION 
FIG. 5.2 VARIATION OF K. FOR PV OF or,, AND 
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5A COMPUTATIOB OS' BROBABIHE- OP PAIlURE AT THAHSPER OP ERESTEESS 

The probabilily of failure of a section at limit state of 
strength at transfer of prestress is given by Bj, 5.2. The value of 
p^ for nozmalJy distributed k is given by 

p^. =i(k < 1) . ♦(-1^) (5.11) 

where and are the parameters mean and standard deviation of k . 

5.5 RELIABILITT ^AITOIS OP A SBmS SCJKOETED ESC BEiM AT STRMG-TH 
LIMIT STADE AT TRARSPER OP ERESTRESS 

The following st^s are involved in the reGLlability analysis of 
a singly supported ISC beam at limit state of strength at transfer of 
prestress. 

(i) Por the given beam, detenoine Ihe critical section where 
maximum compressive stress occurs. 

(ii) Using Monte Carlo method, generate samples for k of the 
section determined in the above step. 

(iii) Using Bj* 5.11, compute p^ of the section, 

!Qie procedure is illustrated vdth an example. 

Example 5.1 

The alTnply supported ISC beam of Example 3.1 is subjected to 

a total dead load (including self weight of the heaa) of 1500 kg/m. 

ahe beam is designed witii concrete strength of 350 kg/cm and steel 

2 

atraigth of 15000 kg/oBi . Ihe effective span and area of steel are 
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!I9ae method of evaluation of probability of failure at this limit stato 
for the above two cases is givsa below. 


Case 1 

Under the action of dead load and prestress, critical sections 
(otlier than end sections) of the beam ate determined as explained in 
article 3»9« Bach span can have a failure mode under the action of 
dead load and prestress and each critical section has a failure 
probability. Hence if there are K critical sections in the ith failure 
mode, the probabilily of occurrence of the failure mode is given by 
the probability of the intersection of events defining the failure 
of each critical section. The probabilily of occurrence of the mode i 
at the limit state of strength at transfer of prestress, p^, is given 
by 

p^ = p(z^nz^ (5.12) 

Assuming events are independent. 




where 'value of p^ of the section 

evaluated as explained article 4.6* If there are 
bounds on p^^ is givoa by 


(5.13) 

Ihls oan be 
m failure modes. 


m 

Ptil^tsli, Pti 

1.= I 

where is the probability of failure of the contimous beam 
at limit state of strength at transfer of pi^stress. 
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Case 2 

The complete collapse of the beam takes place when first 
yielding occurs at aqyoaBof the end sections. 3he value of 
of this section can be computed as explained in article 5.5. The 
value of p^ of the beam for this case is equal to p^ of the end 
section. 

The procedure is illustrated with examples. 

Example 5.2 

ahe two q)an continuous PSC beam of Example 3.2 is 8ub;Jected 

to a total dead load of 1500 kg/m (including self weight of the beam). 

2 

33ie beam is designed wilh conca^ete strength of 350 kg/cm and steel 

2 

strrength of 15000 kg/cm , 5he effective span and area of steel are 
2 

lOmand I0,89cm . The cross section of the beam auad cable profile 

are shown in KLg, 3.12 whLdi is again given here for ready reference. 

She end sections (A and E) are rectauagular with dimensions 50cm x 56can. 

Parameters of a and a„ remain same as in Erample 3.2. 
cu s 

Cfritioal sections of the beam for case (i) failure are C, 
and (Eig. 3.12) and section A(or e) for case (ii) failure. Bead 
load bending moments at sections C, and B^ are l8.'Rtia, 10.39tm 

and tO*39tEi respectively. Using distributions aiai paiometers of 
cr^ and Og and Bjs. 5.4 to 5.10, 10000 samples are generated for k 
of sections C, B^(or B^) and A by Monte Carlo melhod. Parameters 
of K of the above sections obtained frm Monte Carlo method are given 


below, 



i97 




TWO :DHTI?njOU5 BEAM- EXAMPLE 3.2 


Ay.l A^l ffv, wrn 
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Section G s < = 2.738 } a = O.386 

m < 

Section B^(or ) ' •'m = 2.416; = 0,532 

Section A (or E ) : ic^ = 5.9l8; = 0,806 

Ceise 1 failure 

Using Bi. 5.11, the value of of the section C, is 


PtC = 


<!>(■ 


1.0 - 2.738 >. 

0.386 ^ 


3.4 X 10 


-6 


and section is 



4(- 


1.0 - 2.416 

0.332 


■) = 1.0 


X 10 




Using Bj, 5.13, the probability of occura?ence of the failure mode in 
sp^ P^ACj* ^ 

"tAO ' ftO • ftB, = 

Hae value of p^g is equal to p^^^ . 

Bounds on the value of p^g, tising Eb.» 5.14, is given by 
3.4 X 10 l.Pts »<10^^ 


Case 2 failure 


Ehe value of p, of the section A is computed using Bi. 5.11 

"C 

and equal to 5 x 10 Bais is higher tiian the value of p^g 

obtained frcm case (i) failuxe. Hence the probability of tiae failuxo 

of the beaa at limit state of stiongth at transfer of prestress a^ 

“10 

equal to 5 x 10 . 
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Example 5*3 

Dae tharee span continuous PSC beam of Example 3.3 is subjected 

to a total dead load (includir^ self wei^t of the beam) of lOOCkg/m. 

2 2 

0^^ and Og are distaributed as n( 422.8, 56)kg/can and N’(l5680,488)kg/c‘n 

2 

a?espectively. Dae effective span and area of steel are 8*. and 6.94caa 
respectively. Dae caross section of the beam and cable parofile aa?e 
shown in KLg. 3*13 vdaich is again given here for ready refearence. Dao 
end sections A and G- are rectangular with dimensions 40caii 44cm. 

Dbe pixacedure is same as illustac^ited ip the previous example, 
Oarf.txcal sections of the besm for case (i) failaare aire C, D, E 
and (RLg. 3*13) and section A or G for case (ii) failua:^. Dae 
values of dead load bending moments at the above critical sections 
are taken ficon table 3.20. Using Monte Carlo method, 10000 samples 
are geneacated for k of sections B^,C, D and A and parameters of k 
obtained from the above method axe given in table 5.1. Die calculated 
values of p^ and the bounds on p^^ obtained as illustrated in the 
previous example are also ^ven in the same table. Erran table 5,1, 

1±ie bounds on p^^ for case (i) failure are, 

1.04 * 10"’'° £p^g £2.08 X 10"'’° 

ibr case (ii) failure, the value of p^ of the section A is 
computed ard is equal to 4.1 x 10”^°, Diis is greater Hian -fee upper 
bound value of 2.08 x 10”^° for case(i) failure. Hence the probabillly 
of failure of the beam at limit state of strength at transfer of 
prestr^s is equal to 4.1 x 10”’°, 
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Table 5.1. ReHabilily aiJalysis of the beam of ihg, 3,13 at limit state of strength 
at txansf er for PV of o and a 




1 .04(-10) is read 
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5.7 lEDBABILITr DISaHIHIHOH OP ic SDR lEOBABIIISOIC miAEOlirS 

OP 0C5U » ^ EtMENSIONS OP SECHON 

In the previous articles of this Chapter, only probabilistic 
variations of and o^ have been conside 3 ?ed in the reliabilily 
einalysis of ISC beams at limit state of strength at transfer of 
prestress. Random variations of dimensions b^, b^, t^, t^, d, b% 

h and I) of section are also now taken into account in the deteimi- 

s 

nation of probability distribution and parameters of ic . 33ie PSC 
section shown in Pig. 3.11 is considered, 33ie mean and standard 
deviations or o^. o^, b^, b^, b', h, a and a« available 

in Ibcample 3.4s> Using normal distributions and parameters of all 
the above random variables 10000 samples are generated for k of the 
section by Monte Carlo method. Die histogram of tiie generated data 
is gLv^ In Pig. 5.3. It is found that normal, lognormal, !]ype I 
extremal (largest), aype II extremal (largest) distributions do not 
satis:Qr ohi square test at one percent level. However normal an d 
lognormal distributions were suggested by Cornell, Aog and lind 
(37, 38, 39) for (H/S) in the risk based evaluation analysis of 
R.C.C. beams. 

In this chapter k r^resents the ratio of (h/s). Hence, 
as already normal distribution has been observed for ic auad used 
in article 5.4 for probabilistic variations of and Og, here 

also to be in consistent with, the same distribution is adopted for k 
for probabilistic variations of Og arsd dimensions of section. 
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Bor the assijmed normal distribution, the value of can be 

computed by Bl. 5»11» Kie procediire for Idle reliability analysis of 
a sinqiny supported PSC beam is same as explained in article 5.5. Ihis 
is illustrated with an example. 

Example 5.4 

illie singly supported PSC beam of Exanple 3.4 (ELg. 3.1 1) is 
considered vfith probabilistic variations of Og smd dimensions 

of section. Other particulars of the beam and parameters of all 
randcmi variables are same as in Example 5.4. 

'Jixe critical section is the midspan of the beam, Bor this 
section, samples have already been generated for k and shown in 
Eig. 5.3* Erom the same figure the following values are taken, 

= 1,723 and s^ = 0.296 

in ^ 

-3 

lEhe computed value of p^, using Ej, 5,11, is 7.5 x 10 . 

Ihe probability failure analysis of continuous PSC beams for 
probabilistic variations of Cg and dimensions of section is same 

as explained in article 5.7 and illustrated in Examples 5,2 and 5,3 
except that random variations of dimensions of section are also to be 
taken into account. She parameters of the random variables of the 
sections of two span (Eig. 3.12) and three span (Eig. 3*13) beams are 
available in tables 3.9 and 3. 13* Ehe iresults of the reliability 
analysis of the above beams at limit state of strength at transfer of 
presteress for probabiUslio variations of Og 8^ dimensions of 
section are presented in tables 5.2 and 5*3* 



Sable 5.2. Eeliabilily analysia of the beam of Hg. 3.12 at limit state of 
streagth at transfer for PV of o , a and dimsosione of section 
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* 8.2(-5) is read as 8.2 
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12(-4) is read as l,T2 x 10“' 



CHAPTER 6 


REKEABIIiITr BASED DESIGN OP IRESffiESSED CONCRETE 
BEAMS POE LIMIT STATE OP SOEENGTH 

6.1 INTRODUCTION 

The design of sin^jly supported and continuous PSC beams for a 
given probability of failure for limit state of strength is presented 
in this chapter. A sin^de reliability based design method is given 
for two oases considering probabilistic variations of; 

(i) Uve load and strengths of materials and 

(ii) live load, strengths of materials and geometric properties 
of section, 

ilie 3celiability based design formulation reduces to the problem 
of solving an integral equation and finding out the required mean value 
of the resistance of section for given probability of failure and known 
coefficieot of variation (ov) of the resisting moment. 

6.2 REElABIEITf BASED DESIGN OP SIMPET SUPPORTED BEAMS 
6,2.1 Introduction 

She probability of failure of a sin^ily supported beam is ©qual 
to the probability of failure of its critical section, Ohe probability 
of fa 

Pfo 


1 
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= l(l'lu) P(U) (6.1) 

= p(p|0) P(0) (6.2) 

and p^^ are probabiliti-es of failure of the eectton for under- 
reinforced and oveivreinforced cases respeetiyely. In ary design 
problem, the p^ of ttie beam orOy is specified. If Q and are 

lognormally an d normally distributed and p^ and coefficient of 

variation of are known, the only unknown in 5ig. 3*51 is M • 

3he integral equation can be solved and Idle value of be obtained. 

A section is to be proportioned fbr -ttiia mean value of 

6.2.2 Reliability Based Design Method for Probabilistic Variations 
of live load and Strengths of Materials 

In Chapter 3» the values of P^q Pf critical 

section of the different beams (Hgs, 3.11 to 3*13) for different 
loading oases have been computed. Diese values are listed in table 
6,1, Jhe values of »diich is defined as the ratio of p^ to p^, 
for each case is oaloulated aid given in the same table. 3iie average 
value of this ratio is found to be 0,3477. 'Sae coefficient of 
variation of M of each critical section of the beams in Examples 
3.1 » 3.2 and 3.3 are also given in the same table 6.1. It can be 
seen that ooeffioiaat of variation of is almost same in all 
cases and Ihe average value of coefficient of variation of is 
0,03097 vdiidh is approximately equal to Ihe coefficient of variation 
of strength of steel, X^, vdiose value is 0,0311. 
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HaviDg established the values of and eo efficient of 
variation of the following are the steps for the reliabililgr 

baised design of sinpiy supported beams fOr limit state of strength. 


(i) Bor given p^ of the beam, calculate p^ by multiplying p^ 

by Xg . 

(ii) Using the value of coefficient of vaiiatiou of as 0.031 

obtain s_, in terms of 

ru rum 

(iii) Substituting above values of p^ and e^ in Sg., 3.51, solve 
the integral equation and obtain 

(iv) If the section is given, find out the value of area of steel 

for the value of M obtained in the previous step, 

rum ^ 

It is assumed that the section is given and the area of steel 
is to be found out* Sbr flanged sections, the value of is 

given ly tvo Bqs, 3*13 and 3*15 depending on whether neutral aais is 
in the Hange or web* 

If neutral aads lies in the flange. 


M 


rum 


°sm "^ts 


V^cuu ^ 


Ibr giten section and M , the above equation Is 

rum' 

equation of !Ehe quadratic equation can be 

value of is given by 


(6.3) 


a quadratio 
solved and Idle 


id»re. 


(6.4) 
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0.75 a. 


D. = 


am 


1 a 


cum 


M 


rum 
^2 ^ d cr 


am 


If the neutral axia lies in the vreb, 


M 


rum 




+ a A. d(l - 
an Tsw ' 


0.75 Aj. 


sw an 


b‘d cr 


) 


(6.5) 


cam 


inhis is also a quadratic equation in if M and sectional 

^sw rum 

properties are known. !Dhe value of Aj.g^ is given by 

-'>V4 ) 

5 


vdaere. 


0.75 a 


^ “ b'd( 


sm 


cum 


0-T bOtJd-o.5t^)J /(a»^) 


"4 ""“rum cuin""t 

total area of steel is given by 


^ts *" \sw "^tsf 


•sshesre. 


■S:sf ** 




(6.7) 


(6.8) 


em 


Initially it be assumed "ftiat neutral axis lies in the 
web and llie value of may be caloulated. If this Is negative, 

is given by Bl. 6,2 
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3!ie method is illuBtrated by an example. 
Smsple 6*1 


A eiii5)3y supported prestressed concrete beam of effective 

spaa lOmis subjected to a dead load of 1500 kg/m (including self 

wei^t of beam) and probabilistic live load distributed as 

TjN{ 608,85 kg/m, 0.368). and are distributed as 

£f(422.8,56) kg/cm^ and N(15680,488) kg/cm^ respectively, Hie 

section of 1iie beam is given in ELg, 5.11. Deteimine 1iie area of 

—8 

steel if the p^ of 14ie beam is 10 , 

Bae value of p^ is obtained as 

Pfu “ ^ Pf = 


using coefficient of variation of = 0.031 


s = O.C^l M 
ru rum 


aibstituting p^ and s^ in B?.. 3-51, the equation beccanes 


0,3477 X 10 




log( 


rma rma g 




qm 


) 

}] 


(6.9) 


She values of M , M and s are 7,6l1aa, I8,75tm and 0.368. She 
qm^ g q 

above equation is solved in the following way, Ihe iteration is 
started with an initial value of and the ri^t hand side portion 
of the Iq, 6,9 is evaluated using numerical qtiadrature method 
(illustrated in Chapter 3 )* Depending on whether the evaluated value 
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is smallei or greater than the value on the left hanri side of the 

equation, the value of is decreased or increased coirespondingly 

in the next iteration. Iteration is stopped when the difference 

between the values on the right hand side and left barri side of the 

Bl. 6,9 is less than 0,008 percent of the left hand side value of 

the equation. Using the above procedure, the Hi. 6.9 is solved and 

the value of M is found to be 83,82 tn, 
rum 

2 

Using Bq, 6,6, the value of is obtained as 2.4cm . Since 

this is positive, the neutral axis lies in the web, Ohe value of 

2 

is computed using Bq. 6.8 and is equal to 8.34 cm . Hence the 
value of 

\st = W*'? 

Omapaiing Ihis wilh Example 3*5, it can be seen 1hat for an 
2 -9 

area of steel of 11,83 em , Ihe p^ was 1.7 x 10 , Since in iiiis 

■*8 

problem has been increased to 10 , 1die area of steel has been 

decreased as expected. 

6.2.3 R^Liabili-tgr Based Design Melhod for Probabilistic Variations 

of live Load, Strengths of Materials and Dimensions of Section 

Design procedure is same as explaxned in the previous article 
except that values of and coefficient of variation of are 
different* denotes the ratio of p^ to p^ for probabilistic 

variations of live load, strengths of materials and dimensions of 
section. Die values of and p^ calculated for different 
critical sections of beams for different loading oases, are listai 
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in table 6.2. Hie average value of is found to be 0,5847» 

Similarly the values of coefficient of variation of obtained by 

Monte Carlo method, for the same sections for probabilistic variations 

of Og, b^, b', t^, d and in Chapter 3 are also given in the 

same table. It is found that the values of coefficient of variation 

of M are almost same in all cases and the average value is 0,C543. 

Ihis is valid fbr the values of the parameters of the above random 

variables used in this Ihesis, In a general case, for given mean. 

values and standard deviations of the random variables o^, b^, b', 

t, ,d and D„, the mean value and standard deviation of M_, can be found 
t s' ru 

laj applying partial derivative meUmd (59) to Sg.. 3.15 or 3*15 as boHi 
yield almost same value. 

Using Ihe above established values of X^^ and coefficient of 

variation of M_ , the value of M is computed for given p» of the 
ru' rum ^ ■^i 

beam as e:q)lained in the previous sirticle. In this csise also the same 
equations (i.e, Bjs. 6.3 to 6,8) are valid except that in every equation 
meeua values of the randcan variables b*, t.^, and d are to be used 
along with o and a , The procedure is illustrated wilh an 
example. 

Biample 6.2 

Hte supported TSC beam of Hxaaqple 6.1 is considered. 

It is aibjected to probabilistic variations of live load, strengths 
of materials and dimensionB of section, DtstrlbutionB and parameters 



Table 6.2. Values of and CV of for different sectiooe for PV of Q 
diitffinsiotB of section 
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of Q, and remain same as in Example 6,1, !Ihe paraittetere and 

distributions of random variables b^, b’, t^ and d of this beam are 

available In Example Determine the mean value nf area of eteel 

-8 

if the p^ of idle beam is 10 . 

Using = 0.5847 (established earlier)* the value of 
Pfu = * 

UsJoig coefficient of variation of « 0JD545 (established earlier), 


s^, = 0,0545 
lu lurtt 


aibstitaung and 8^ Id Bj. 3.51 and solving the sans, th. valus 

of M is obtained as 83.66tm, Substituting the mean values of 
irum 

b^,b’,t,^ and d in Eg,. 6,5, the mean value of is computed as 

2 

1,12 am • Since litis is positive, mean value of neutral axis lies 

2 

in web. Using Eg, 6.8, the mean value of is obtained as 8.63cm . 

Hence the mean value of area of steel, is given by 


^:am “ 

dhie is less than the value of steel used in the Emaple 3.6. Ihe 
beam has been designed in this problm for higher vai^e of pj SaSi 
as expected, less area of steel is required. 


6.3 RBHABIHBT BASED DESIGN OP ODNETBOOTS BEAMS 
6.3.1 Introduction 

Continuous beanss have a number of failure modes and are 
subjected to different live load conditions. 'Caking a two span 
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continuous beam, there are three live load conditions inhich gire 
majcinmi positive and negative bending moments in the beam. Even 
though oliier load conditions i.e, partial loading of a span or spans 
are possible, ttae failure of the beam for such conditions will be 
veiy very small and insignificant compared to the values of probability 
of failtare of the beam for the three load conditions. Hence, 
considering only 3 loading oases as in Example 3.2, it is seen that 
Idle probability of failure of the beem for all load conditions and 
failuire modes vaiy from one to two times the value of lower bound (ib) 
of Pjg* I't is also noticed that IB value of corresponds to 
the probability of occurrence of a particular failure mode for the 
most critical loading condition. Hence, if it is possible to find 
out the probability of occurrence of the particular mode, for vhidi 
probability of failure is maxinum under the most critical loading 
condition, from given p^^ of the beam, Ihe beam can be designed for 
this failure mode and its probability of failure. 

6.3,2 Beliability Based Design Method for Continuous Prestressed 

Oonorete Beams 

SBaese are the following steps that are associated in Ihe 
reliability based design of continuous prestressed concrete beams 
for limit state of strei^fth. 

(i) Assume the lowet* bound value of probability of failure of 
•file beaa eqjial to the given p- of the beam (Oiis is on the 
cMinaervative side). 
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(ii) Find oat the most critical loading condition for the given 
heam ■which gives maxiirum moment at any critical section. 

Find out the critical failure mode, i,e. the prohahility of 

occurrence of which is the highest, and its critical sections 

for this loading case. The prohahility of occurrence of this 

failure mode is assumed as the value of p„ . 

IS 

(iii) ^'ix the prohahility of failure of the section which yields 
first (this is fixed as 10”^ or 10”^) in the critical 
failure mode. Hence the prohahility of failure of the other 
critical section in the same failure mode (say, there are two 
critical sections in the failure mode) is the prohahility of 
occurrence of the failure mode divided hy the prohahili'ty of 
failure of the section yielding first 

(iv) Assume effective depth of the heam at the faxst yielding 
critical section (this is fixed as maxiimm as possible) . 

Design the section, i,e. for prohahility of failure of 

the section, fixed in step (iii), based on the method explained 

in article 6,2,2, 

(v) For the calculated prohahility of failure of the second critical 

section, find out the value of of the section as explained 

in article 6.2.2, For this value of 11 _ and calculated A, „ 

rum ts 

in step (iv), find out the effective depth of the section 
using Bqs, 6,3 or 6,5 whichever is applicable. Adopt a 
parabolic cable pccofile passing through the cable posi-fcions 
fixEd at the critical sections. 
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The procedure is illustrated with examples. 

Example 6,3 

A oontimous prestressed concrete lieam having two equal spans 
is subjected to a dead load of 15OO kg/m (including self weight of 
the heam) and probabilistic live load distributed as LF(608.85 kg/m, 
0.36 b). o and o are distributed as 1^(422.8,56) kg/cm and 

CTil S 

2 

N( 15680,488) kg/om respectively. Other particulars remain same as 

in Example 3*2 (Eig. 3 •'12), The probability of failure of the "beam 
*•1 1 

is fixed as 1O , 


The LB value of the p^^ is taken as equal to the given proba- 

-11 

bility of failure of the beam i,e, 10 

The critical loading condition for the "beam is when both spans 
are fully occupied by live load. The first failure is at section C. 
The critioal failure mode occurs in span AC or CE, Considering 
failure mode in span AC, tt® two critical sections are C and 
(Refer Exsimple 5.2) , The probability of occurrence of the failure 
mode in span AC, is taken as equal to the LB value of Hence 


PfiC ■ ■"> 


FixitJg as 10"^, the pcroha'bility of failure of the sactiou is 


pfp» 

Fixing effective deptti of beam at C as 50 oms, the area of steel 

—4 2/ 

required for equal to 10^ is 7 •23 om (o'btained as illustratai 
in Emisple 6,1 ). 
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-7 

Using «* 10 (calculated earlier), the value of M 
1 

from Eq, 3,51 is comiuted as illustrated in Exan^le 6,1 and is 
equal to 41 .48tm, Using Eq, 6,3» 


c 0.75 K ^ 

41.84 « 1o5 = A. d [1 --4-Jt2_25 ] 

ts sm Dj-d a 

t cum 


Substituting the values of A , a , h and a , the value of d 
^ ts' sm' t cum' 

is obtained as 44 •92cm, Using this d at section the section 

is ohecJced whether neutral axis lies in the flange. It has been 
checked ard found that neutral axis lies in the flange, A parabolic 
cable profile is adopted. 


Example 6,4 

A continuous prestressed concrete beam having three equal 

spans is subjected to a dead load of 1000 kg/m (including self 

weight of the beam) and probabilistic live load distributed as 

LH(487»08 kg/m, O.368), Mean and standard deviations of and 

and other pairticulars of the beam remain same as in Example 5«3 

(Fig, 5«15), The probability of failure of the beam for all load 

-11 

conditions and failure mode is fixed as 10 , 


The lower bound value of the p^^ of the beam is taken as 
equal to 10 . 

The critical loading condition is when live load occupieB 
two adjacent spans say AC and CE. For this loading condition the 
critical ihdlure mode occurs in span AC, The critical sections in 
this failure mode are and C( Fig, 3»'53)» The probability of 



221 


ooourrenoe of this failure mode in span AC is taken as equal to LB 

>-11 


value of i,e. p, 


fAC 


10 


“4 

Fixing p^ = 10^ (Section C yields first), the probahxlity 


of failure of the sections, is 

1 


"fB 


10 


-7 


1 


Fixing effective depth of beam at section C as 39 cm, the area of 

1 ^ 2 

steel required for the value of p^^ equal to 10 is 5*95 cm . 


Using the above value of A^^ and p^ , the value of M 


1 


rum 


is computed, using Eq. 3 *51 and method explained in article 6.2.2, 

as 23.42tin. The value of effective depth at B. for this M , 

1 rum' 

using Eq, 6,3> is obtained as 40.75om. This is greater than tie 
maximun possible effective depth i.e, 39 cm . This can be solved in 
two wayss 

(i) By inoreaaing the probability of failure of section C and 
thereby reducing the value of p^ , This will require more 

fBi 

area of steel for the sane d„ and less d„ . 

C B^ 

( ii) By usixig the maximum value of dg , required A^^ for the value 
of « 23 . 42 tm can be computed. This will result more 
area of steel and hence the probability of failure of section 

_ -4 

0 will be less than 10 which will satisfy the design criteria. 


By using the second method, the value of has been obtained 


as 4»11om for dg 39cm and 25.42tm. A parabolic cable profile 


rum 


is adopted 
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Sinilaxly the same heama of Exan^les 6.5 and 6,4j subjected to 
protahilistic variations of live load| strengths of materials and 
dimensions of section, oan he designed for a given prohahility of 
failure except that individual sections are to he designed based on 
the procedure explained in article 6.2.3» 

The design method is simple. The solution is possible if 

the coefficient of variation of the resisting moment of the section 

for under-reinforced case is known. In the design fornfiilation the 

mean value of the ultimate resisting moment of the section for 

under-reinforced case is first determined for given reliability of 

the beam. Assuming section is given, the method of finding area 

of steel for calculated M has been illustrated. If a section is 

rum 

to be designed, it can be done using the ultimate strength design 

of a PSC section fcr required mean value of the resisting moment 

using mean values of o and o . 

^ ou s 



CHA.PEER 7 


SHII -PROBABILISTIC LDCET STATE DESIGN OP PHESTRESSED 
CONCRETE BEAMS 


7 . 1 Introdviction 

The philosophy in structural design is to provide a stmcture 
which can serve the required purpose under certain external forces 
with an assured degree of reliahility. 

Structural design methods can "broadly "be classified into the 
following: 

1 , Detecministio design method - In this method, design parameters, 
load and resistances are considered as deterministic and the 
safety is ensured "by limiting action less than resistance. 

2. Pro'ba'tdlistic design method (60) - In this method "basic strength 
and load parameters are treated as random and the safety 
conditions are eacpressed hy specifying a priori the probability 
that the structure will preserve its functional capacity ov®c 

a given period of time, 

5. Semi-probabilistic design method - In this method loads and 
resistances are fixed from statistical analysis and assigned 
pro"bability of failure. The safety is enaired "by limiting 
the action less than resistance. 

In the send-Hprobabilistio limit state design melhod, safety is 
assured "by the selection of vaiioos coefficients (for the strengths 
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of mterials, loads etc.) dependent on the limit state considered. 

The -rolues of these coefficients, wherever is possible, are proposed 
in this ohaiJter based on the statistical analysis of materials and 
reliability analysis of KG beams from Chapters 2 to 6, An example 
is presented illustrating the semL-probabdlistio design method. 

7.2 CHABACTERISTIC STRENGTHS OP MTERIALS 

The characteristic strengths of materials are, by definition, 
those which have an agreed probability of not being obtained. It is 
given ly 

( 7 . 1 ) 

where 

= characteristic strength 
= mean strength 

s^ » standard deviation of strength 

The CEB-PIP Committee recommends the value of k as t.64 (9^ probability) 

for strengths of concrete and steel. In Chapter 2 the value of k 

for o has been establidied as 1,6 for the field data. This is 
ou 

oloa© to the value suggested by CEB-PIP Committee. Henoe the same 
vabue, 1,64 (safer side), is recommended for k for strength of concrete, 
Henoe 

o , = o - 1.64 s (7 •2) 

ctjk cum ^ c \i ■/ 

Regarding the value of k for the statistical analysia of 

field data on o in Chapter 2, shows wide range of values for k 
s 

varying from O.C® to 12, IJumbex of projects being less, more data 
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are necessary on a to arri-re at the value of k for a for Indian 

s s 

conditions. Till then, the same value of I.64 recommended by 
CEB-PIP Committee may be used. Hence 


o, = o - 1 ,64 s (7 .5) 

^ am ^ 8 ^ 

7.5 DESIGN STRENGTHS OP MATERIALS 

•j 

Reduction factors — are applied for various materials, to the 
characteristic strengths, to take into account the uncertain factors 
and those factors -which cannot be easily evaluated and measured. 

For a given material the value of Y depends on the behaviovir of 
material and seriousness of risk involved in attaining limit state. 

The calculated reduced s-trength, R^, is called as the design strength; 


R, 




(7.4) 


CEB-PIP ConBnittee recommends y = 1 for o aid. at serviceability 

cu s 

limit states and Y = 1.15 and 7 = 1 .4 to 1.6 (depending on 

q^uality control) for at ultimate limit states. The same values 
are used in this thesis. 


7.4 CHARACTERISTIO LOADS 

For loads where a possible increase could be detrimental 
the ohairaoteristio load, is defined liy 

where 

*. the value of the most unfavourable loading wildi a fifty 
percent probability of its not being exceeded during the 
expected life of the structure. 
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k =■ coefficient depending on the agreed probahility of loading 
greater than 

6 s> coefficient of variation of the maxLnnxm loading. 

w 

If on the other hand a reduction in the value of load will endanger 
the stahility of the structure, the characteristic load, is 
defined as 

W.' =. w’ (1 - k 5 ) (T,S) 

1C m '' W'd; 

where 

«» the value of the most unfavourahle loading with a fifty 
percent prohahility of its not falling helow in the 
expected life of the structure, 

W is obtained from statistical analysis of loads on structures, 
m 

2 

The load survey (Chapter 2) shows a mean live load of 150*55 kg/m 

in office rooms at IIT, Kanjur. The I.S, Code (55) specifies 250 
2 

to 400 kg/m for office "buildings. Depending on this the value of k 
also varies. Local government and state and central government 
office structures in India are subjected to higher loads ats compared 
to IIT, Kanpur "buildings, Henoe it will not be possible to fix the 
value of k unless extensive load data on similar office "buildings are 
available. When adequate statiatioal data on loads are not available, 
the minifflam loads specified "by codes are taken as characteristic values. 

7 .5 DESIGN LOADS 

Design loads, axe o"btained "by multiplying the oharact eristic lo 
if Isad factors corresponding to particular limit state. 



227 


W, = F W, 
a j 


There F is load factor. 


(7.7) 


P. Srinivasa Rao and C.S. Kriahnamoorthy (61) discussed the 
criteria for fixing up the load factors for dead and live loads. 

Two partial load factors for live load, one for the vaidation of 
live load and the other for covering uncertainties over the exact 
distrihution of live load, errors in analysis and for the accepted 
prohability of failure, were introduced. The load factor for live 
load was expressed as the product of the above two partial load 
factors for live load. They concluded that even though the strength 
factors adopted are different in various codes, the load factors are 
reconiEiended such that moment carrying capacities of PSC singly 
reinforoed beams under working loads are found to be nearly the same. 

Load factors differ from code to code to a limited extent. The 
actual margin of safety depends not only on load factors and material 
reduction factors but also onihe quality conirol specifications. 

Even, if the load factors are large and if also the coefficient of 
variation of material strength is largajthen the margin of safety 
is not proportional to load factors. A rational appxoeush to load 
factors is to fix a probability of failure of the section and then 
airrive at the load factors. Probability of failure of a section 
can be stated as 

^ T W I 


(7.8) 
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As per tiltimate strength design, the value of is 


4 j=M ssFM + FM 
r u g g q. q 


M - F(M+M) = FM 
r g q' o ew 


(7.9) 


where F is the combined load factor. 

0 

Assuming normal distribution for resistance, the pjcobability of 
failure can be expressed as 

M - 

. / ew rm\ 

♦ ( = Pf 


M - M 




or M « M + k s 
rm ew r 


(7.10) 


Ualrig s * 5 M * the abo-ve equation is rewritten as 
^ r r rm ^ 


M (l - k6 ) = M 
rm ^ x' ew 


(7.11) 


Substitution of Eq. 7*9 in 7.11 gives 


5’ “ (T:r^") 

0 ^ 1-4c 0 ^ 

* r 


(7.12) 


Pig. 7.1 ^ows the variation, of P^ and p^ for different values of 6. 

Separate load factors for dead load and live load are obtained 
as follows. Using Eq. 7.11» 


(F M+P M)(l-k6) - M+M 
^ g g a ^ ^ S % 

M 


P^ - [1 + V (1^^+ kd^ Pp] /(I - k«^) 


(7.13) 
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where v « « Hence load factor for live load depends on 

(i) the ratio of dead load hcnding inonu3nt to live load bending 
moment, (ii) load factor for dead load, (iii) the value of Ic which 
con^esponds to the level of probability of failure fixed and 
(iv) coefficient of variation of resisting moment. 

Pigs, 7 *2 to 7 A show the variation of P and p It can be seen 

Si 

from the above figures that 


(i) 


(ii) 

(iii) 


for constant values of v and P , P increases if ( a) p» 

g q. \ 

decreases and(b) 5^ increases for constant value of 

for constant values of 6 and P , P increases if (a) p„ 

deoireases and (b) ^ increases for constant 

for constant values of 5 and v , P increases if (a) p- 

X 9, 1 

decreases and(b) P^ decreases for constant . 


Prom the reliability analysis of PSC beams (Chapter 5) at limit 

state of strength it is observed that coefficient ot variation of 

ultinate resisting moment of a PSC section varies from 3 to 8 percent. 

If the random vaadations of geomefcrio properties of a section are not 

considered, is about 3,1 percent for under~reinfbiDed oase and 

about 5 #5 percent for over-reinforced case . Generally beams are 

designed as uidex-reinforced. However there is a probability of 

the sai^ beams becoDdng over-reinforced. Hence is assumed as 

5 percent fear fixing load factor. In building floors, the value 

of V varies from 0,5 to 1,5* “faking v «= 0.75 and fixing p^ » 

and P •» 1,2, the value of P , using Eq,. 7*15, is found to be 1,4* 

g 1 
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where ve> (M Hence load factor for live load depends on 
(i) the ratio of load bonding moment to live load bending 
momenti (ii) load factor for dead load, (iii) the value of k which 
corresponds to the level of probability of failure fixed and 
(iv) coefficient of variation of resisting moment. 

Pigs, 7*2 to 7*4 show the variation of P and p™. It can be seen 
from the above figures that 


(i) 


(ii) 

(iii) 


for constant values of v and P , P increases if (a) p„ 

g q. \ 

decreases and(b) 5^ increases for constant value of 

for constant values of d and P , P„ increases if (a) p, 

r g' q. \ ^ 

decreases and (h) v increases for constant 

for constant values of 6 and v , P insreases if (a) p, 

X ^ I 

decreases and(b) P^ deoreanes for constant p^ . 


Prom the reliability analysis of PSC beams (Chapter 5) at lisdt 
state of strength it is observed that coefficient of variation of 
ultimate resisting moment of a PSC section vsiries from 3 to 8 percent. 
If the random variations of geometric properties of a section axe not 
consider^, is about 3»1 percent for under-rein&iDed case and 
about 5 *5 percent for over-reinforced case. Generally beams are 
designed as under-reinforced. However there is a probability of 
the saai!© beams becoming over-reinforced. Hence 5^ is assumed as 
5 peroent for fixing load factor. In building floors, the value 
of V varies from 0.3 to 1.5. Taking v « 0.75 and fixing * 10 
and P " 1.2, the value of P , using Eq. 7»13t is found to be 1.4* 


/ 
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For prestressed concrete "baildings P i? 1 .2 and P - 1 «4 are 

© ^ 

recommended. 


At limit state of cracking, the coefficient of variation of the 
resisting moment of a PSC section is ot®erved to be about 2,8 percent 
(Chapter 4) for probabilistic variation of c and a • To be on the 

OU 8 

safer side, it has been taken as 3 percent for fixing the load factor. 

Taking P = 1.1 as per I.S. Code of practice on prestressed concrete 
S 

(46) and fixing = 10 , the value of P^, using Eq. 7*15» is found 

to be 1.216 which is less than 1.25 specified by the I.S. Code. 


7.6 CIMEACTERISTIC MD DESIGN HfflSTRESS 

The characteristic value of the prestressing force, 
defin€Jd by the relation 

V (7.14) 

where is the mean value of the prestressing force and 5^ is 
the coefficient of variation of P arrising out of different sources 
of exrcr. 

The design value of the prestressing force, Pd, is givoa by 

(7.15) 

In these formlae, the positive sign (and as a corollary > I) 

or the negative sign (and < 1) is chosen according to whether 

an inoircase or a decrease in prestress is detrimental. CBB-PIP 

Coamittee su^ests Y » 0,9 for limit state of strength under 

F 

maxitaim loading ooisiition and Y^« 1 for limit state of strength 
under •glrrimum loading condition. 
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7*7 CHECK ON SAKETY 


For the limit state undGr consideration the effects of design 

loads mast he, at worst, equal to the absolute value permitted hy the 

design strengths of the materials. The external bending moment 

caused by the ultimate load is a function of dead load (DL), live 

load (LL), load factors and geometry of the structure. The 

strength of the section is a function of strengths of materials, 

material reduction coefficients and geometric properties of. the 

section. Hence the characteristic formula is written in the following 

manner for limit state of strength. 

a a 

V^g \ ^ LL) 1 ^ts) 

c s 

where Y and Y are material reduction factors for concrete and 
o s 

steel respectively. For statically determinate structures, 

M (F X DL), M (F X 111), etc. are equal to F M (DL), F M (liL),eto. 

Henoe Eq, 7 •16 becomes 


U 

F M + F M < M (— ^ V , d, t. , A. ^ ) 

g g q q- r^ Y, ’ Y * t* t» ts 


(7.17) 


a o 
IS £ 

r ’ T 

0 s 

Similar equations can be written for any o-Bier limit state. Knowing 
characteristic strengths of materials, mterial reduction factors, 
load factors etc., the semi-probabilistio design and analysis of 
PSC beams can be done as usual. This is illustmted with an. 


example 
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ExBJmple 7 ‘I 

A simply supported PSC "beam of effective span lOm is sutijeoted 
to a dead load of 15OO kg/m (including self weight of the heam) and 
a live load of 2000 kg/m. 

Given : 

h » 50 5 h « 60 cmj 

^2 2 
a 422 .Skg/ cm ? ^0 * kg/cm 

cum ' 2 2 

a s 156GOkg/cm ; s » 4G8kg/cm 
sm s 

P=>1.2 5Fa1,4 

g q. 

The heam Is designed based on limit state of strength at mfod-Eaim 
loading ( ultimate) at transfer of prestress. 

Lindt state of strength at maxinum loading: 

Top flange of the heam and area of steel are designed based on 
this limit state. 

ultimate load = P xLL+PxiJj 

s q. 

• 1.2 X 1500+ 1.4 X 2000 
*» 460Ckg/m 

ultimate external bending moment is 

2 

» 4600 57.5 tm 

o 

Using Etj, 7.2, the oharaoteristio straagth of concrete, 

%k “ - 1.64 X 56 - 551 kg/on? 

Similarly, using Bq^. 7*5» characteristic strength of steel 
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i» I568O - 1,64 ^ 488 » 1488O kg/om^ 

The design value odf strength of concrete, is given hy 

%uk 

'^cud * y 
0 

using Y =• 1.5 
0 

a , = « 220. 7kg/ cm^ 

cud 1,5 ' ^ 

Similarly, using Y « I.I5, the design value of strength of steel, 

s 

o^, is oomiuted 

‘^sd * 'kg/ cm 

The effective depth of heam xs fixed as 54ciQS leaving 6cm3 from the 
hottom of the team to the centroid of steel. 


Asaxming lever arm at failure as 0,95d i.e. 0.95 ^ 34 - 51 *5 o™e 
and the failure stress in steel is equal to its ultimate strength, 
the ultinnte resisting moment of the section is given "by 


M « T X lever arm 
r 8 


T 


s 


X 51.5 


From Eq, 7»15» 


Heiu>e 


F M + F M < M 
Eg 9 q, “ r 

57.5x10^151.3 T 


T > ^ tonnes 

8— 51.5 


To equal the above value of T^, the cross seotioreil area of ocaKsrete 
in oompeession mist he 
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113.2 xio^ 
“ 220.7 


513 cm 


For neutral axis to lie inside the flange, apjnroxiinate requiiBd 
thioknoss of flange at top is 

t. =» " 10.3 cm 

t 50 

H om. is adopted. 

It is now necessary to chedk the assumptions concerning , 
failure stress in steel and the position of neutral axis. 

The heam is designed as under-reinforced. Following I.S, speci- 
fications (46), the failure stress in steel is taken sia equal to its 
ultimate strength for under-reinforced case. Assuming the neutral 
axis lies in the weh (Fig. 7*5)» 


30 - 12)11 X 220,1 

15000 


7.07 cm 


Resisting moment of overhanging portion of the flange, is 

M^-7.07 X 13000 (54-0.5 x1l) 

• 46tm 

Moment to he resisted hy weh, is given hy 


« 57,5 - 46 « 11.5 tm 
Equation horizontal forces (Fig, 7»5)» 


0 « 1)* a «» 0 A A. 

0kA sd tsw 


(T.18) 





?l 
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|4r, 
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(.;) r,F-bi(>N O'- TOP FiANGfF - LIMIT 5 TATC OF STRtNCTI- 
AT MAXIMUM LOADING 



it! OEGIGM OF BOTTOM FLANGE- 
LIMIT 5 TA 7 E OF 5 TRENOTH 
AT TRANSFER OF PRE - 
STRESS 



FiG.7,5 EXAMPLE 7.1 


(c) SECTION 
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Equating’ moments and using a reduction factor of 0,9 » 


0-5 ^tsw - f) - \ 


(7.19) 


substituting Eq, 7 *18 in Eq. 7 . 19 , 

A o 

^ ^ ^ « =i A _tswad-v 


21)' do w 

cud 


A X 13000 

0.9 X 13000 « 54 (1 - s^nbwtJao:? ’ ' ■"•5 


Solving tlK) above equation, the -value of is found to be 2om , 


\m °sd 2 X >15000 

b'do^d " 12 X 54 x 220.7 


0,182 < 0.24 


Hence the beam is under-reinforced and the failure stress in steel 
is equal to its ultimate strength. The value of area of steel is 


A. » 7 .07 -t 2 » 9 .07 om^ 

uS 

Limit state of strength at transfer of prestress; 


The bottom flange of the l)eam is designed tased on this limit 
state . 

As an increase In prestress is detrimental to the beam at the 
limit state of -transfer of prestress, the characteiistio strength 
of steel is gi-Tsn by 

^sk “ ‘^sm ^ 

Taking the value of stress in steel at transfer as 0,8u^, the 
oharcKsteristio pre stressing force at transfer, P^, is 
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"dk ^ts 

» 0.8 X 16480 X 9.07 = 119500 kg 
The design value of initial prostijessing force is given ly 


td 




using Yg = 1» the value of is 119500 kg* Decrease in dead load 

is detrioental to the "beam at this limit state. Design dead load, 

M , , is taken as 0,9 times charactenstio dead load. Hence 
gd 

10^ 

^gd ” ^ ^500 X = l6.R75tni 


Compression area req,viired 


^oud 


545cm 


Assuming (h-t, ) > a (Fig, 7 *515) and equating the distance of centroid 
^ M 

of compression area from the centroid of steel to , the following 

nd 

equation is obtained: 


■b'a(f - x„) - ( yb') - 0.^) 

P P 

/td ^ ^td 

oud 

M 

or Va(0,5a - 3^)-('b^-b' )t^(3!:^-0.5\) « . 

cud 

Adopting t^ =» t^ = Horn , 

12a(0.5a-^) - (b. -12)12(6-5.5) « 

220,7 

Segleoting ( b^- 12) 12(6-5 *5) » being small, the above equation is 
solved and the value of ' a’ is found to be 42 ,25 om, fhis is less 
ttan (h-t^) which is equal to 49 ob. Bence tbs assun^tion a< (h-t^) 


is satisfied. 
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Equating the resisting capacity of the section to the prestressing 
force, 


or 


tVa+ (Vl=')t^]. P 


td 




td 


cud 


12 X 42.25 + (hjj- 12)11 « 543 


or ^ “ 15-30II1 

But on practical considerations, the hottom width of flange is 
adopted as 20cm . The designed section is shown in Pig, 7*5o» 


The semi-prohahilistic limit state design of a PSC heam. has 
"been illustrated for strength limit state and strength limit state 
at transfer of prestress. Similarly other limit states can also 
be considered and the design can be done. The semi -probabilistic 
limit state design method is suited for practical design purposes. 



CH/lPTER 8 


DISCUSSIONS AND CONCLUSIONS 


8.1 GENERAL 

Ihe iBsults of the reliability analysis and design of PSC 
beams, presented in the previous chapters, are discussed herein. 

Prom the outcome of the results, conclusions are enumerated. 

Finally suggestions for further woiic are indicated. 

8.2 RESULTS OF THE STATISTICAL ANALYSIS OP STRENGTH OP COHCKETB 

The strength of concrete varies as a random phenomenon suboeotod 
to normal or lognormal distribution. The frequency distributions of 
the strengths of M550 and M2 50 concretes axe normal with five percent 
level of significanoe. In the case of M200 concrete, where good 
supervision was provided while concreting, the distribution is normal 
with one percent level of significanoe. The distribution of the 
half of ihe groups of MI50 concrete have satisfied the normal 
distribution with one percent level of significanoe while the other 
half g3?ou.ps have not. However, all the groups together have 

satisfied lognormal distribution with five percent level of signi- 
ficance. The skewness in Ml 9) concrete is due to the fact that 
it Is baaed on the nominal mix proportions with limited supervision. 
The concrete samples of MI50 belong to many diffeirent projects 
and the quality of concrete is not of reliable nature in eome of 
the projaots. Considering all these aspects it can be oonoladed 
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that the strength of concrete follows nomal distribution with atleast 
one percent level of significance it miniraum control on the quality 
of the concrete is ensured. 


Ninely one to ninety seven percent of Ml50, M200 and M 350 concrete 

camples satisfy the I. S. Code specifications. Only ei^iypercent of 
H250 concrete samples satisfy the Indian Code sjjecifications. The mean 
value of the strength of M150 samples has been found to be much higfier 
than the characteristic value; therefore most of the samples satisfied 
the code specifications even though variance of the strength is large. 

One of the most inqportant questions to be answered in the seai- 
probabilistic design is, what are the eaqjected mean strength arai 
variance of a concrete in the field, for a specified characteristic 
value, Ihe -three values can be connected thou^ a linear relation as 
recoimaended ty the Offi-flP Committee as 


a , = o 
cut cum 


(1 -U) 


in ishich the value of k assigns the probability of failure. The -values 
of k for different groups of concrete vary from 1.21 to 2.01, The least- 
square fit of the field data of several groups of concrete indicates 
that the prob£i)ility of failure of -the concrete is 5.5 percent. The 
■value of k corresponding to the above probabili-fy of failure is 1,6. 
Kjis is close to the value of 1 ,64 (probabiUL-ty of failure is 5 percent) 
rewaamended 1y COT-KCP (SMmittee (45). 
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8.3 R33SIILTS OP THE STATISTICAI. MANSIS OP STRMGTH OP STEEL 

The strength of steel is foxmd to be a random variable foir.o’«i.ng 
nomnal diBtilbution, The frequency distribution of strength of hiqh 
tensile steel wires is nownal wi1ii five percent level of significano - 
The proof strength of cold twisted deformed steel has also normal 
distribution with one percent significance level. The coefficient cf 
variation of strenglh of 7 nm 4 HTS wire and CTD steel are 3-1 4-7 

percent respectively while it has been observed to be lees 1 han oo*' 
percent in the case of 5 nm EPS wire. As the value of the ratio of 
mean strer^h of 5 mm 41 HTS wire to its characteristic strength is high 
and the coefficient of variation is small, the probability of failure 
of the strength of 5 mm <)> HTS is found to be less than The 

probability of failure of 7 ma HIS wrire is found to be abo^rfc 0.03. 

Hie probability of failure of CTD steel (diameter varying frcm 25 mm t*^ 
32 nm), based on yield strength, is 0.25. 

8.4 RESULTS OP HIE STATISTICAL ANALfSIS OP HOOK LOAD IH OPPICE R0«S 

The frequency distribution of the floor load in office rocmn 
is fCHind to be lognormal with five percent significeince level, Hio 
equivalent load intensity varl^ froo 56 kg/m to 293 kg/m . The 
coefficient of variation of load is observed to be about 0.4. The 
probability of load greater than tte specified load depends on the 
characteristic load. Eor characteristic load of 400 kg/k , the 
probability of lead greater than "Uie above value is fovekl “to be 

6 X . 
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2 

Ittie mean value of tiie floor load is 130.77 kg/m . Even tbou^ 
this observed mean value of the floor load is small compared to 1116 
minimum loads specified by the I.S. Code fbr office buildings, the 
coefficient of va3?iation of the observed load is hi^. Both parameters 
mean and coefficient of variation of the load are important as both 
have influence on the probabilily of failure of the structure. More 
extensive load data in similar office buildings is neoessaiy to arrive 
at the characteristic values of the loads with certain confidence 
level so that Ihey can be introduced into the code. 

8.5 RBSTITS OP !KE STATISTICAl MAIMSIS OP GBOMEHtlC PROPERTIES OF A 
SEOEEOH 

It has been observed that dimensions of a section are 
random variables with small coefficient of variation. All parameters 
of the section are normally distributed with atleast one percent level 
of significance, DSie coefficients of variation of widlh. of flange, 
breadth of web and total depth of the section are negligible and 
significant variation has been observed in thickness of flange and 
effective depth of tbe beaa. 

8.6 REailTS OP ms REEEABIIITr AMITSIS OP ISC BEAMS AT UMIT 
STATE OP STBEHOTH 

It has been observed that even if a PSC section is iinder- 
reinforced based on deteamiinlstie analysis there is a probability of 
the seotLon becoming ovei>-reinforced. If probabilistic variations of 
and Og are considered, the probability of the most critical sections 
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of the beame in Examples 5.1, 3.2 and 3.3 becoming ove3>-reinforced 
is about 0.16, 0*10, and 0,124 respectively. Dae corresponding 

*^tsTv 

values of the ratios ) of the above sections are 0.234» 

cuk 

0,208 and 0,224 respectively. Hence, the probability of under- 
reinforced sections becoming over-reinforced increases with the 
increase in area of steel as expected. If the probabilistic 
variatiom of a ,* cr„ aod dimensions of section are comidered, 
the probability of the same sections becoming over-reinforced is 
about 0,07 , 0.04 and 0,05 respectively whida are less than the values 

if probabilistic variations of and only are considered, Hiis 

t "^tsiv^sk \ 

is because that mean veilues of the ratios { —. - n , ) are smaller as 

b 'do , 
cuk 

the mean values of ^ ' and d are hi^er. 

51ie probability of failure of a sia 5 )ly supported PSC beam 

(designed by woriiug stress method as per I,S, Code) for deteiministic 

-7 

load and probabilistic variations of and is i.08 x 10 . 

It is found that even thou^ the probability of the section becoming 

over-reinforced is aaall, the probability of failtrce of the section 

depends on tto seotion becoming over-reinforced ani failure taking 

place. If probabilistic variations of and eJri dimensions 

of section are considered, the probability of failure of the saaare 
-7 

beam is 1,4 x iO , It can be concluded that ttie probability of 
failure of slr^ly supported PSC beams (designed on vnrking stress method) 
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-7 

at limit state of strength at design load is of the order of 10 
for deterministic load. 

Bor probabilistic load, the probabililgr of failure of the 

-Q 

siB^idy supported PSO beam is 1,68 x 10 for probabilistic variations 

-10 

of and Cg aiid 1,84 x 10 for paxibabilistic variations of 
®cu' '^s dimensions of section. Ihese failure probabilities 

are smaller t h a n the values for deterministic load. OMs is because 
the mean value of the load is far less than lixe characteristic load 
when compared with the deterministic load even tbou^ coefficient 
of va3?iation of the load is lao^ge. Here it is fbund that there are 
contributions to p^ values by under-reinforced and over-reinforced 
cases. Hence it can be said that the pirobabilily of the section 
beccaaing over-reinforced is also to be considered Tidiile confuting 
failure probabilities of the beam at limit state of strength. It 

can be concluded that the probability of failure ofsingjty supported 

-9 

PSC beam for probabilistic load is about 10 . 

Die probability of failure of continuous PSC beams shovm in 

3»12 and 3*15 for different cases are given in table 8,1, 

It can be said fKm the above table that probability of failure of 

coatiimuyas PSC.beaas for deterministio load and fhr probabilistic 

-17 

variatiorts of a and is of the order 10 and for probabiliBtio 

variations of o and a„ and dimensions of section is of the order 10 

cu s 


« 



SCabl© 8*1. Probability of failure of contimous pa3 beams at l imi t state of strength 
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5br probabilistic load^ the probabill.ty of failuare of oontinuoua beams for 

-20 

probabilistic Tariations of 0 ^^^ and 0 ^ is of the order 10 and for 

probabilistic variations of and dimensiom of sections is 

-22 

of the order 10 , The probability of failure of continuous 

beams is smaller than probability of failure of simply supported 
beams as expected. 

8.7 R£S[JL3!S OF HE REHAEUITSf ANAJZSIS OF FSC BEAKS AT UMIT 
STATE OP CRACKIMG 

Fbr deterministic load, the probability of cracking of the 

-22 

singly supported beam is 2.7 x 10 when probabilistio variations 

of 0 ^^ and CTg are considered. When probabilistic variations of 

a , a and dimensions of section are considered, it is observed 
cu* a ’ 

that there is a sli^t increase in M but large inoreeise in s_„ 

rcm rc 

which results an increase in the value of p vidda is equal to 
-9 

1.5 x 10 . For probabilistic load, the value of p of tiie BiEQ>ly 

w 

supported beam, for probabilistic variation of ani Cg > is 

2.94 X 10 ”^ and for probabilistic variation of Og sod dirasnsions 

-7 

of section, is 4,87 x 10 . 

The values of p for continuous beams aare given in table 

8*2^ It 1b obsewed that for detemLnistic load the probability 

of cracking of continuous BSC beams vari^ frt» 10 to very small 
-24 

value less than 10 and &r parobabilistic load, it varies ficm 



Table 6.2* Probabilily of cracking of oontLmoiae PSC beams 
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-5 -7 

10 tolO . It is also observed that for piobabilistic load, the 
probabilily of crack i ng of isc beams is greater than the probabi- 
lity of failure of ‘Hae same beams at limit o£ etrength. 

8.8 RESUITB OP !IHE REHABIUTy AHAITSIS OP PSD BEMS AT LIMIT 
STATS OP DEPtEOTION 

It is observed that the probability of the deflection exceeding 
the limit value at transfer of prestress is almost zero fbr simply 
supported and contdraaous PSC beams. This shows a large margin of 
safety. 

Ebr sin^ily supported and continuous PSC beams, the probability 

of the deflection exceeding the limit value at d^ign load is almost 
zero fbr deteuninistic load and probabilistic variations of 9^^ and 

0g, Prom table 8.3, it is observed that for probabilistic load, 

the value of p,, for simply supported PSC beams varies from 
dd 

7.27 ^ to 1.81 X 10 ard for continuous PSC beams from 

10~^^ to 10~^^, The values of for continuous PSC beams are 
smaller tdian the values of for sin5)ly supported PSC beams as 
expected* It is also seen that probability of failure of PSC 
beaaas at limit state of deflection is smaller than the values at 
limit state of strength and cracking, 

8.9 R^LTS op THE RELIABIIiaS IMETSIS OP iKC BEfflilS AT LIMIT 
STATE OP SORMGTH AT TRAHSEffl OP PRKTRESS 

The probability of failure of the sinply supported beam adt 

Halt state of streigth at transfer of pi^strese is 5.52 * 10 
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for probabilistic variations of and Og. When probabilistic 

variations of 0 ^,,. o„ and dimensions of section are also consi- 

dered, it is observed that mean value of k representing the ratio 

of resistance to action, slightly increases and coefficient of 

variation of < increases largely. Daese result increase in the 

-5 

value of to 7*5 ^ 10 for the same beam, 

Ihe values of p^g for continuous PSC beams are given in 
table 8 , 4 . Erom tiie same table it can be said that vdien probabi- 
listic varlatiore of o_, and a„ are considered, the value of p. 
of iiie beam is of 121 e order 10 and ^ifaen probabilistic variations 

of o , Cg and dimensions of section are taken into account, it is 

_8 

of the order 10 , Por contLmxous PSC bems -Oae two iiossible 

failures ares 

(i) the first yielding taking place at a section other than the 
end sections and the beam collapsing after Hie occurrence of 
a failure mode 

(ii) Hie first yielding taking place at end sections and the besaa 
reducing into overhang b«aa which aitomatically loads into 
collapse. 

Both failures are to be investigated to detefrndne the probability 
of failure of continuous PSO beam at limit state of sti’ength at 


transfer of prestress. 
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8. 10 EFFECT OF ON AREA OF STEEL FOR PBDBABIUSnC LOAD 

The effect of failU3?e probability on area of steel and mean 
value of -tile ultimate resisting moment of the section has been 
studied in the case of probabilistic loads. The relation between 
and and p^ and are shown in Figs. 8,1 and 8.2, This 
has been obtained fbr the simply supported PSC beam shown in 
Fig. 5. 11 considerdLng probabilistic variations of and , 

It is seen 1diat an increase in p^ reduces the value of required 
resisting moment of -ttie cross section and area of steel as anticipated. 
Same thii^ is observed when probabilistic variations of a^, 
and geametidc proi>erties of the cross section are tafcen into account. 

8.11 EFFECT OP FAIILRE PRDBABILIffi ON LOAD FACTOES 

Eelation between probability of failure and load factors have 
been studied In Chapter 7 and are given in Figs. 7.1 to 7.4. 

It is seen that the value of combined load factor increases 
with decrease in the value of probability of failure for constant 
coefficient of variation of resistance. It is also obaesrred that 
the value of F^ increases with increase in 6^ for cons'tent p^.. 

It is observed from FLg. 7.2 to 7.4 liiat 

(i) for constant values of v and Fg, F^^ Incsreasee if (a) 
deor^ises and (h) 5^, increases for constant value of p^ 

(ii) for constant values of 8 and. F , F increases if (a) p^ 
decreases and (b) w increases fbr constant p^ 
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(iii) for constant values of 5 and v » 5* increases if 

a. 

(a) decreases and (b) decreases for constant p^. 

Coefficient of va3:lation of ulti m ate resisting nument of PSO 

sections varies fram 3 to 6 percent whraa probabilistic variations 

of a-od Og are considered and from 5.5 "to 7»6 percent ishen 

probabilistic variations of and gecanetrlc properties of 

.6 

tbe section are considered. If safety having probability of 10 
is fixed, the value of P has been found to be 1,4 (chapter 7) 

M 

for 6 = 0,05 and P_ =* 1 .2, Ihia value is close to the value 

I" e 

reconmended by Bussian codes, Sgme values of P »1,2 and P » 1,4 

& Vl 

are recommended for PSD building floors, aihese factors are based 
on the failure of the beam at limit state of strength, 

!Che coefficient of variation of the resisting moment of a P6C 

section at limit state of cracMng is observed to be about 2,8 percent, 

PafcLng its value of 3 percent and fixing as 10 , Ihe value of 

p has been found to be 1,22 for P = 1.1 and the ratio of dead load 
q g 

to live lc«d equal to 0,75. The value of P *= 1.22 is slightly 
l®3s than 1,25 recommended by I.S, Code (46) for limi t state of 
initial oracic. The values of Pg » 1.1 to 1,2 and Pq= 1,2 to 1,3 are 
recommended for Umlt state of initial oracdt, 

8.12 C»HcnPSIOIS AJ© HS«)WEKDATIOIS 

1, Sbe strength of concrete follows noiroal distribution with 
atLeast one percent significance level if nrinlnami control 
in Ihe quality of wmrete is ensured. 
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2» !Ihe most probable value of the failure of concrete cubes is 
5 .5 percent. !I!he value of 1.64 is recommended for liie 

characteristio coefficient which ensures a preassigned proba- 
bility (5 percent) of the test results less than the doaraot- 
eristio strength, 

3, Ihe strength of hi^ tensile steel follows normal distribution 
with 5 percent significance level, 

4, The probabilily of fail\ire of 7mm 4 high tensile steel wire 

is 8,0 percent and of 5mm 4 hi^ tensile steel wire is less 
-24 

than 10 based on the characteristic values specified by Ihe 
I.S. Code (46). 

5, Sloor load in office rooms follows lognormal distribution 
with five percent significance level, She mean value of the 
load is very much less tiian Ihe characteristic load given by 
the I.S, Code (53 )j but its coefficient of variation is high* 

6, Geometric properties of ISO sections follow normal distribution 
with at least one percent level of significance, 

7, Even thou^ a section is under-reinforced based on determincLstie 
analysis, there is a probability of the same secfeLon becoming 
over^reinfbrced. This must be considered in computing 
probability of failure of PSC beams at limit state of str^agth 
at design load. 

At limit state of strength, 'ttie probability of failure of 

-7 

(a) sin^jly supported PSC beams is of the order 10 for 


8. 
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■*9 

detenairdstio load and 10 for probabUiatic load 

(b) continuous ISC beams varies frma 10”^^ to for 

—20 —22 

detemLnistio load and from 10 to 10 for probabilistic 

load. 

9« At li m it state of cracking, tiie probabilily of cracking of 

(a) siii 5 )ly supported ISC beams varies from 10*^ to lo""^^ for 

-6 -7 

deteiminiatio load and from 10 to 10 for probabilistic 

load 

(b) continuous ISC beams varies from 10”^ to very small value 

-24 -5 -7 

less than 10 for detenainistic load and from 10 to 10 

for probabilistic load . 

It is observed that for probabilistic load, the probabiHly of 
cracking of prestressed concrete beams is hi^er than the value at 
limit state of strength, Ibis is because of the fact that for 
almost same combined standard deviation of (e-S) at both limit 
states, the difference between mean values of the resistance of the 
beaaa and action is amall at cradcirg (that is overlap of distribu- 
tions is more) oon5>ared to the large difference between mean values 
of resistance and action at limit state of strength. Sbr deterministic 
load, the same oonolu^on is ea^jeoted provided Ihst the resistance 
of the section has same probability distribution and coefficient 
of variation at both limit state- crac^ng and strength, aowevesr, 
the relative piobahilily of feilnre at limit states of strength and 
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cracking depends on mac^- parameters such as 

(i) coefficient of Taan-ation of resisting mcaaent and external 

moment (ii) lype of distribution functions for resisting moment 

and external moment and (iii) the level of inteieeotion of both 

distributiorffi. In general, the probability of cracking of preedEreaned 

concrete beams is bl^er than the probability of failure at limit 

state of strength, 

10. At limit states of deflection at transfer of prestress, the 

probability of deflection in prestressed concrete beams 

, -24 V 

exceeding the limit value is zero (less than 10 ), 

At limit state of deflection at design load, the probability 
of the deflection exceeding the limit value 

(a) for sin^jly supported and continuous prestressed concrete beams 
As zero (less than 10 fbr deterministic load 

(b) for singly supported prestressed concrete beams varies frcan 

—10 ~12 

10 to 10 and confeimous prestressed concrete beams frcm 

10”^^ to 10 fbr probabilistic load. 

11. At limit state of strength at transfer of prestress, the 

probabilily of failure of singly supportei prestressed concrete 

beams vartes from 10*"^ to 10*^ and confeinuoas prestressed concrete 

-8 -10 
beams from 10 to 10 . 


12. JProbabilily of failure at tiransfer condition of prestressed 

concrete beams is fairer than ttkat at limit states of strength 
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arid deflection under woiking load condition. The probability 
of cracking of prestressed concrete beams is generally hi^er 
than the probability of failure at limit states of strength 
and deflection under woiteing load condition. 

8.13 SOGOESTIONS TOR H3RTHEE TORE 

More data need, be collected on the strength of steel so Ihat 
the value of characteristic coefficient, ?*ach ensures the pressigned 
probability of the test results less than the characteristic strength, 

can be fixed with certain confidence level and recommended to I.S. 

Code, Extensive survey of floor loada in s im ilar office buildings 

nood be conducted to arrivo at rational voluos of charcctoristic 
loads which can be includod in the code. 

Other limit states such as failure of anchorage zone (end bloclcs), 
buckling, torsion eto« are also to be oonsictered in the re lia b ili ty 
analysis and design of PSC beams. 

As it IS known that e^, 6^.^., <3^ and aM mutually iadependent, 
it is siiggested that data may be collects on '5^, and mutual 

dependence of the above variables may be established. Hence further 
research on the reliability analysis of prestressed concrete beams mj be 
carried out taking into account of mutual dependence of the abo c v-; L‘ Ic" 
and random variations of prestressing operations. 
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AlPEHDIX 


Table 

1. Sjuivaleat unlfomly distributed load In office rooms 
in IIT, Kanpur 

SI. 

H’o. 

Room 

No, 

Office of 

Area 

2 

m 

Total 

load 

kg 

Bjuivalent 

TTOL 

kg/m2 

1 

101 

U,G. Acad«nic 

147-8 

19321 

131 

2 

102 

Mministrat ive 

147,8 

18140 

123 

3 

105 

Accounts 

147.8 

17038 

115 

4 

104 

Accounts 

57.8 

aCH9 

139 

5 

105 

Aocoimts 

65.1 

8218 

126 

6 

208 

Registrar 

46.5 

5456 

117 

7 

205 

lostgraduate 

28 

3656 

131 

8 

210 

Postgraduate 

46.5 

6576 

141.5 

9 

211 

Faculty 

administration 

46.5 

4756 

102 

10 

204 

Audit 

28 

2548 

91 

11 

276 

Dean of 
faculty 

46.5 

6276 

135 

12 

212 

Institute works 
aepartment 

69.7 

13368 

192 

13 

220 

Estate 

28 

2329 

83.3 

14 

219 

Import section 

46.5 

7976 

168 

15 

218 

Director 

46.5 

7334 

158 

16 

214 

Di3reotor 

28 

1569 

, 56 

17 

215 

Director 

28 

3243 

115.5 

18 

223 

Dean of 1 & D 

28 

2435 

87 

19 

275 

Placemmt 

46.5 

4696 

101 
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SI. 

Fo. 

Roam 

No. 

Office of 

Area 

2 

m 

!Dotal load 
kg 

Bjuivalent 

UDL 

kg/m^ 

20 

274 

Standing 

Oommittee 

46.5 

8456 

182 

21 

263 

Elect. dept. 

23.2 

3388 

146 

22 

313 

Civil dept. 

46.5 

5415 

116.5 

23 

332 

Aero. dept. 

46.5 

6143 

132 

24 

564 

Mech, dept. 

23.2 

4888 

210 

25 

381 

Haysice dept. 

23.2 

5896 

254 

26 

412 

Met, dept. 

28 

1868 

67 

27 

432 

(Jhemical dept. 

46.5 

5582 

120 

28 

462 

CSaeoiical dept. 

28 

1942 

69.5 

29 

461 

Chemioal dept. 

14 

2228 

159 

30 

513 

Maths, dept. 

46.5 

4189 

90 

51 

567 

Mucation 

development 

centre 

23.25 

2548 

109-5 

32 

613 

Humanities dept. 

46.5 

5062 

109 

33 

WL 121 

Dean of students 

46.5 

4896 

105 

34 

WI 122 

Asso. Dean of 
students 

23-25 

2022 

87 

35 

273 

Purchase 

46.5 

7364 

158-5 

36 

221 

Director 

14 

4101 

295 

37 

201 

Mason 

23-25 

14® 

64 

38 

OS 101 

Coaputer centre 

20 

3329 

166.5 



